Identification and Quantification of Polymorphism in the Pharmaceutical Compound Diclofenac Acid by Terahertz Spectroscopy and Solid-State Density Functional Theory

Using terahertz spectroscopy to quantify bioactive flavonoids in Moxa Wool as predictor of rheumatoid arthritis treatment outcomes

BACKGROUND

Moxibustion, a traditional Chinese medicine practice, employs Moxa Wool, derived from Artemisia argyi. Flavonoids, the key pharmacological constituents in Moxa Wool, are known for their anti-inflammatory and analgesic properties. The purity of Moxa Wool, particularly its flavonoid content, directly influences the efficacy of moxibustion treatments. However, quantifying these bioactive flavonoids accurately and non-destructively has been a challenge.

PURPOSE

This study introduces terahertz spectroscopy as a non-destructive optical detection method for qualitative detection and quantitative analysis of flavonoids in Moxa Wool. By establishing a mathematical model between spectral signals and clinical efficacy, a reliable correlation between flavonoid concentration and the therapeutic effect of moxibustion can be established, providing a potential predictive model for the treatment outcomes of rheumatoid arthritis.

STUDY DESIGN

We adopted terahertz spectroscopy technology and combined it with terahertz metamaterial biosensors to achieve rapid, efficient, and non-destructive testing of the quality of Moxa Wool. This method reduces the detection time from hours to minutes while lowering the sample detection limit, overcoming the limitations of traditional detection methods in pharmacological research.

METHODS

Through terahertz metamaterial biosensors, rapid detection of the purity of Moxa Wool has been achieved. A combination of molecular simulation and terahertz spectroscopy was used to quantitatively analyze the flavonoid content in different purities of Moxa Wool. To ensure accuracy, the quantitative results of flavonoids obtained by terahertz spectroscopy were validated using high-performance liquid chromatography (HPLC). In addition, moxibustion treatment was performed on rats with rheumatoid arthritis using Moxa Wool, and medical indicator information was recorded. A mathematical analysis model was established to evaluate the correlation between flavonoid content and analgesic and anti-inflammatory effects.

RESULTS

Terahertz spectroscopy analysis shows that there is a direct correlation between the flavonoid content in moxibustion and the absorption peak intensity. The maximum R2 in the model analysis is 0.98, indicating a high accuracy in predicting the purity of Moxa Wool. These results were also validated by HPLC. In a rat model, the purity of 30:1 Moxa Wool samples showed a 50 % decrease in TNF-α, IL-1β, and IL-6 levels during treatment compared to low-purity samples, significantly reducing inflammation markers and pain symptoms. Meanwhile, The PLS prediction model established a correlation between terahertz-detected flavonoid levels and treatment outcomes (PWL and IL-1β). The maximum R2 in the model is 0.91, indicating a high correlation between flavonoid levels and the anti-inflammatory and analgesic effects of moxibustion treatment.

CONCLUSION

This study not only demonstrates the effectiveness of terahertz spectroscopy in the pharmacological quantification of bioactive compounds but also establishes a novel predictive model for the efficacy of moxibustion in rheumatoid arthritis treatment. It underscores the potential of integrating traditional medicine insights with advanced technology to enhance therapeutic strategies in pharmacology.

Single-Shot Direct Transmission Terahertz Imaging Based on Intense Broadband Terahertz Radiation

There are numerous applications of terahertz (THz) imaging in many fields. However, current THz imaging is generally based on scanning technique due to the limited intensity of the THz sources. Thus, it takes a long time to obtain a frame image of the target and cannot meet the requirement of fast THz imaging. Here, we demonstrate a single-shot direct THz imaging strategy based on a broadband intense THz source with a frequency range of 0.1~23 THz and a THz camera with a frequency response range of 1~7 THz. This THz source was generated from the laser-plasma interaction, with its central frequency at ~12 THz. The frame rate of this imaging system was 8.5 frames per second. The imaging resolution reached 146.2 μm. With this imaging system, a single-shot THz image for a target object with a size of more than 7 cm was routinely obtained, showing a potential application for fast THz imaging. Furthermore, we proposed and tested an image enhancement algorithm based on an improved dark channel prior (DCP) theory and multi-scale retinex (MSR) theory to optimize the image brightness, contrast, entropy and peak signal-to-noise ratio (PSNR).

Quantification of resistant starch content in rice after hydrothermal treatments using terahertz spectroscopy

Since hydrothermal treatments can enhance resistant starch (RS) content in rice and provide health benefits when consumed, a less laborious and non-destructive method to determine RS content is needed. Terahertz (THz) spectroscopy is hypothesized as a suitable method to quantify RS content in rice after hydrothermal treatment with its sensitivity for the intermolecular forces increase in the formation of RS. In this study, we first used the traditional in vitro hydrolysis method to determine the content of RS in rice. Then, the potential of starch absorbance peaks to quantify RS content after three commonly used hydrothermal methods, soaking, mild heat-moisture treatment, and parboiling, was investigated. The second derivative intensities of the peak at 9.0, 10.5, 12.1, and 13.1 THz were confirmed as being correlated with RS content and showed the high accuracy to predict RS content in samples (R2 > 0.96). Our results indicate the RS content of hydrothermally treated rice can be accurately quantified using these peaks.

Harmonic and anharmonic studies on THz spectra of two vanillin polymorphs

Polymorphism commonly exists in organic molecular crystals. The fingerprint features in low-frequency vibrational range are important information reflecting different intermolecular interactions of polymorphs. Interpreting these features is very helpful to understand vibrational property of polymorphs and reveal the thermodynamic stability. In this work, the low-frequency vibrations of form I and II of vanillin are investigated using terahertz time-domain spectroscopy. Static DFT calculation and ab initio molecular dynamics (AIMD) are employed to interpret their low-frequency vibrations of both forms in harmonic and anharmonic ways, respectively. Their low-frequency vibration characteristics in harmonic calculations are discussed, and anharmonic mode couplings between OH bond stretch and the stretching and bending motion of hydrogen bonds are uncovered. Moreover, the thermodynamic energies including electronic potential energy and vibrational/kinetic energy arising from nuclear motions are calculated. The result reveals that the stability order of the two forms is mainly dependent on their electric potential energy difference.

Frequency-domain terahertz optoacoustics for non-contact quantitative detection of gas, liquid, and solid samples

Terahertz optoacoustics (THz-OA) combines the advantages of abundant molecular characteristic absorptions in a terahertz band and the low attenuation through ultrasonic detection. Frequency-domain THz-OA, benefiting from the compact and the low cost of a continuous-wave THz source, has been used in gas detection and sensing. However, liquid and solid detections are hard to achieve due to the sensitivity limitation of existing technologies. Here we present a high-sensitivity frequency-domain THz-OA system with customized optoacoustic cells to accomplish non-contact quantitative detection of gas, liquid, and solid samples. The relationships between signal amplitudes and sample concentration, volume and temperature are discussed separately, revealing a potential application of this technology.

In Situ Cocrystallization via Spray Drying with Polymer as a Strategy to Prevent Cocrystal Dissociation

The aim of the present study was to investigate how different polymers affect the dissociation of cocrystals prepared by co-spray-drying active pharmaceutical ingredient (API), coformer, and polymer. Diclofenac acid-l-proline cocrystal (DPCC) was selected in this study as a model cocrystal due to its previously reported poor physical stability in a high-humidity environment. Polymers investigated include polyvinylpyrrolidone (PVP), poly(1-vinylpyrrolidone-co-vinyl acetate) (PVPVA), hydroxypropyl methyl cellulose, hydroxypropylmethylcellulose acetate succinate, ethyl cellulose, and Eudragit L-100. Terahertz Raman spectroscopy (THz Raman) and powder X-ray diffraction (PXRD) were used to monitor the cocrystal dissociation rate in a high-humidity environment. A Raman probe was used in situ to monitor the extent of the dissociation of DPCC and DPCC in crystalline solid dispersions (CSDs) with polymer when exposed to pH 6.8 phosphate buffer and water. The solubility of DPCC and solid dispersions of DPCC in pH 6.8 phosphate buffer and water was also measured. The dissociation of DPCC was water-mediated, and more than 60% of DPCC dissociated in 18 h at 40 °C and 95% RH. Interestingly, the physical stability of the cocrystal was effectively improved by producing CSDs with polymers. The inclusion of just 1 wt % polymer in a CSD with DPCC protected the cocrystal from dissociation over 18 h under the same conditions. Furthermore, the CSD with PVPVA was still partially stable, and the CSD with PVP was stable (undissociated) after 7 days. The superior stability of DPCC in CSDs with PVP and PVPVA was also demonstrated when systems were exposed to water or pH 6.8 phosphate buffer and resulted in higher dynamic solubility of the CSDs compared to DPCC alone. The improvement in physical stability of the cocrystal in CSDs was thought to be due to an efficient mixing between polymer and cocrystal at the molecular level provided by spray drying and in situ gelling of polymer. It is hypothesized that polymer chains could undergo gelling in situ and form a physical barrier, preventing cocrystal interaction with water, which contributes to slowing down the water-mediated dissociation.

Determining changes in crystallinity of rice starch after heat-moisture treatment using terahertz spectroscopy

To investigate the potential of Fourier-transform terahertz (FT-THz) spectroscopy to follow crystalline structure changes in rice starch after heat-moisture treatment (HMT), we measured the crystallinity by X-ray diffraction (XRD) spectra and found its correlation with THz spectra. According to A-type crystal structure and Vh-type crystalline structure of amylose-lipid complex (ALC) in rice starch, crystallinity is divided into A-type and Vh-type. The intensity of second derivative spectra peak at 9.0 THz was highly correlated with both A-type and Vh-type crystallinity. Additionally, other three peaks at 10.5 THz, 12.2 THz, and 13.1 THz were also sensitive to Vh-type crystalline structure. These results indicate that after HMT, the crystallinity of ALC (Vh-type) and A-type starch can be quantified using THz peaks.

High-Mobility Naphthalene Diimide Derivatives Revealed by Raman-Based In Silico Screening

Charge transport in crystalline organic semiconductors (OSCs) is considerably hindered by low-frequency vibrations introducing dynamic disorder in the charge transfer integrals. Recently, we have shown that the contributions of various vibrational modes to the dynamic disorder correlate with their Raman intensities and suggested a Raman-based approach for estimation of the dynamic disorder and search for potentially high-mobility OSCs. In the present paper, we showcase this approach by revealing the highest-mobility OSC(s) in two series of crystalline naphthalene diimide derivatives bearing alkyl or cycloalkyl substituents. In contrast to our previous studies, Raman spectra are not measured, but are instead calculated using periodic DFT. As a result, an OSC with a potentially high charge mobility is revealed in each of the two series, and further mobility calculations corroborate this choice. Namely, for the naphthalene diimide derivatives with butyl and cyclopentyl substituents, the estimated room-temperature isotropic electron mobilities are as high as 6 and 15 cm2 V-1 s-1, respectively, in the latter case even exceeding 20 cm2 V-1 s-1 in a two-dimensional plane. Thus, our results highlight the potential of using the calculated Raman spectra to search for high-mobility crystalline OSCs and reveal two promising OSCs, which were previously overlooked.

Pharmaceutical Cocrystals of Ethenzamide: Molecular Structure Analysis Based on Vibrational Spectra and DFT Calculations

Pharmaceutical cocrystals can offer another advanced strategy for drug preparation and development and can facilitate improvements to the physicochemical properties of active pharmaceutical ingredients (APIs) without altering their chemical structures and corresponding pharmacological activities. Therefore, cocrystals show a great deal of potential in the development and research of drugs. In this work, pharmaceutical cocrystals of ethenzamide (ETZ) with 2,6-dihydroxybenzoic acid (26DHBA), 2,4-dihydroxybenzoic acid (24DHBA) and gallic acid (GA) were synthesized by the solvent evaporation method. In order to gain a deeper understanding of the structural changes after ETZ cocrystallization, terahertz time domain spectroscopy (THz-TDS) and Raman spectroscopy were used to characterize the single starting samples, corresponding physical mixtures and the cocrystals. In addition, the possible molecular structures of ETZ-GA, ETZ-26DHBA and ETZ-24DHBA cocrystals were optimized by density functional theory (DFT). The results of THz and Raman spectra with the DFT simulations for the three cocrystals revealed that the ETZ-GA cocrystal formed an O−H∙∙∙O hydrogen bond between the -OH of GA and oxygen of the amide group of the ETZ molecule, and it was also found that ETZ formed a dimer through a supramolecular amide–amide homosynthon; meanwhile, the ETZ-26DHBA cocrystal was formed by a powerful supramolecular acid–amide heterosynthon, and the ETZ-24DHBA cocrystal formed the O−H∙∙∙O hydrogen bond between the 4-hydroxy group of 24DHBA and oxygen of the amide group of the ETZ molecule. It could be seen that in the molecular structure analysis of the three cocrystals, the position and number of hydroxyl groups in the coformers play an essential role in guiding the formation of specific supramolecular synthons.

Classification of Amino Acids Using Hybrid Terahertz Spectrum and an Efficient Channel Attention Convolutional Neural Network

Terahertz (THz) spectroscopy is the de facto method to study the vibration modes and rotational energy levels of molecules and is a widely used molecular sensor for non-destructive inspection. Here, based on the THz spectra of 20 amino acids, a method that extracts high-dimensional features from a hybrid spectrum combined with absorption rate and refractive index is proposed. A convolutional neural network (CNN) calibrated by efficient channel attention (ECA) is designed to learn from the high-dimensional features and make classifications. The proposed method achieves an accuracy of 99.9% and 99.2% on two testing datasets, which are 12.5% and 23% higher than the method solely classifying the absorption spectrum. The proposed method also realizes a processing speed of 3782.46 frames per second (fps), which is the highest among all the methods in comparison. Due to the compact size, high accuracy, and high speed, the proposed method is viable for future applications in THz chemical sensors.

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.

Diclofenac Ion Hydration: Experimental and Theoretical Search for Anion Pairs

Self-assembly of organic ions in aqueous solutions is a hot topic at the present time, and substances that are well-soluble in water are usually studied. In this work, aqueous solutions of sodium diclofenac are investigated, which, like most medicinal compounds, is poorly soluble in water. Classical MD modeling of an aqueous solution of diclofenac sodium showed equilibrium between the hydrated anion and the hydrated dimer of the diclofenac anion. The assignment and interpretation of the bands in the UV, NIR, and IR spectra are based on DFT calculations in the discrete-continuum approximation. It has been shown that the combined use of spectroscopic methods in various frequency ranges with classical MD simulations and DFT calculations provides valuable information on the association processes of medical compounds in aqueous solutions. Additionally, such a combined application of experimental and calculation methods allowed us to put forward a hypothesis about the mechanism of the effect of diclofenac sodium in high dilutions on a solution of diclofenac sodium.

The fingerprints of nifedipine/isonicotinamide cocrystal polymorph studied by terahertz time-domain spectroscopy

Cocrystal is constructed to improve physicochemical properties of active pharmaceutical ingredient and prevent polymorphism via intermolecular interactions. However, recent examples on cocrystal polymorphs display significantly different properties. Even though some analytical techniques have been used to characterize the cocrystal polymorphic system, it remains unclear how intermolecular interactions drive and stabilize the structure. In this work, we study the cocrystal polymorphs of nifedipine (NFD) and isonicotinamide (INA) using terahertz (THz) spectroscopy. Form I and form II of NFD-INA cocrystals show spectral fingerprints in THz region. Temperature-dependent THz spectra display distinguished frequency shifts of each fingerprint. Combined with solid-state density functional theory (DFT) calculations, the experimental fingerprints and their distinct responses to temperature are elucidated by specific collective vibrational modes. The vibrations of hydrogen bonding between dihydropyridine ring of NFD and INA are generally distributed below 1.5 THz, which play important roles in stabilizing cocrystal and preventing the oxidation of NFD. The rotations of methyl group in NFD are widely distributed in the range of 1.5-4.0 THz, which helps the steric recognition. The results demonstrate that THz spectroscopy is a sensitive tool to discriminate cocrystal polymorphs. It has the potential to be used as a non-invasive technique for pharmaceutical screening.

Comparison of Proton Acceptor and Proton Donor Properties of H2O and H2O2 in Organic Crystals of Drug-like Compounds: Peroxosolvates vs. Crystallohydrates

Two new peroxosolvates of drug-like compounds were synthesized and studied by a combination of X-ray crystallographic, Raman spectroscopic methods, and periodic DFT computations. The enthalpies of H-bonds formed by hydrogen peroxide (H₂O₂) as a donor and an acceptor of protons were compared with the enthalpies of analogous H-bonds formed by water (H₂O) in isomorphic (isostructural) hydrates. The enthalpies of H-bonds formed by H₂O₂ as a proton donor turned out to be higher than the values of the corresponding H-bonds formed by H₂O. In the case of H₂O₂ as a proton acceptor in H-bonds, the ratio appeared reversed. The neutral O∙∙∙H-O/O∙∙∙H-N bonds formed by the lone electron pair of the oxygen atom of water were the strongest H-bonds in the considered crystals. In the paper, it was found out that the low-frequency Raman spectra of isomorphous crystalline hydrate and peroxosolvate of N-(5-Nitro-2-furfurylidene)-1-aminohydantoin are similar. As for the isostructural hydrate and peroxosolvate of the salt of protonated 2-amino-nicotinic acid and maleic acid monoanion, the Raman spectra are different.

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.

Terahertz and Raman Spectroscopic Investigation of Monohydrate Cocrystal of Antitubercular Isoniazid with Protocatechuic Acid

Pharmaceutical cocrystal provides an alternative modification strategy for the formulation development of drugs owning to their potential ability to improve the physicochemical properties of active pharmaceutical ingredients (APIs) efficiently by changing inter-molecular interactions between raw materials. Isoniazid (INH) is an indispensable main drug for the treatment of tuberculosis, but its tablet formulation is unstable and prone to degradation. In the present study, the monohydrate cocrystal of INH and protocatechuic acid (PA) was prepared by solvent evaporation using PA as cocrystal former to optimize the properties of INH. The parent materials and corresponding 1:1 molar ratio INH-PA monohydrate cocrystal have been characterized by the terahertz time-domain (THz-TDS) and Raman spectroscopy. The THz absorption spectra displayed that there were obvious differences between the peaks of experimental cocrystal and the parent materials, and the same situation was found in Raman vibrational spectra. In addition, density functional theory (DFT) was applied to simulating and optimizing the structure of INH-PA monohydrate cocrystal and supplied corresponding vibrational modes. Our results provided a unique method to characterize the formation of INH-PA monohydrate cocrystal at the molecular-level and a lot of information about cocrystal structure and intra-molecular and/or inter-molecular hydrogen bond interactions in the emerging pharmaceutical cocrystal fields.

Ru(II)/diclofenac-based complexes: DNA, BSA interaction and their anticancer evaluation against lung and breast tumor cells

Ruthenium(ii) diclofenac-based complexes of the general formula [Ru(dicl)(P-P)(bpy)]PF6 [dicl = diclofenac, bpy = 2,2'-bipyridine, and P-P = 1,4'-bis(diphenylphosphino)butane (dppb) (1), 1,2'-bis(diphenylphosphino)ethane (dppe) (2), 1,3'-bis(diphenylphosphino)propane (dppp) (3) and 1,1'-bis(diphenylphosphino)ferrocene (dppf) (4)] are synthesized. The complexes (1-4) are characterized by elemental analyses, infrared, NMR, and UV-vis spectroscopy and (3) and (4) are characterized by single crystal X-ray diffraction. The DNA binding of complexes (1-4), studied by circular dichroism (CD) and Hoechst 33 258 staining assay, indicates their binding with the minor grooves. The complexes interact with BSA with binding constants (Kb) in the range of 2.5 × 103-5.5 × 104 M-1. The complexes exhibit high cytotoxicity against the tumor cell lines A549, MDA-MB-231, and MCF-7 with IC50 values ranging from 0.56 to 15.28 μM. The complexes are more selective for the hormone-dependent MCF-7 breast tumor cell line and complex (1) is the most potent one. The study demonstrates the anticancer activity of ruthenium(ii)/diclofenac-based complexes.

Quantitation of Commercially Available API Solid Forms by Application of the NMR-qSRC Approach: An Optimization Strategy Based on In Silico Simulations

Physical forms of active pharmaceutical ingredients (APIs) play a crucial role in drug discovery since 85% of API molecules exhibit polymorphism and sometimes complicated phase behavior, often resulting in important differences in the respective biochemical and physical properties. Characterization and quantitation of the different forms are becoming more and more essential in the pharmaceutical industry: once these characteristics are known, it is easier to choose the best solid form for development, formulation, manufacturing, and storage. Time domain-nuclear magnetic resonance (TD-NMR) has recently been used to develop a quantitation protocol for solid mixtures, named qSRC, based on the linear combination of T₁ saturation recovery curves (SRCs) collected on a bench-top instrument. Despite its potentials and ease of use, a limited number of application cases have been reported in the literature since its development and many aspects remain to be clarified for the technique to be adopted as a robust routinely industrial analytical tool. In the present work, the reliability of the qSRC approach has been studied by focusing on the role played by key experimental variables, including mixture composition, signal-to-noise ratio, and T₁ differences. In silico simulations were carried out for a wide range of theoretical cases to predict the expected level of accuracy obtainable for a given sample-parameter acquisition set and to clearly define the range of applicability of the method. Results of the simulation are presented alongside a comparison with three real-case studies of commercially available APIs: piroxicam, naproxen sodium, and benzocaine.

Determination of starch crystallinity with the Fourier-transform terahertz spectrometer

We measured the terahertz (THz) spectra of native, amorphous, and dried starches derived from corn and potato using the Fourier-transform (FT) system and compared these spectra to the X-ray diffraction (XRD) patterns. Both native corn and potato starches had seven absorption peaks in the terahertz regions, but five peaks were observed in the amorphous states. While spectral changes slightly occurred in corn starch even after drying, increase and decrease in the terahertz peak intensities were obtained in potato starch during drying. Similar changes in both starches during amorphization and drying were obtained in the X-ray diffraction patterns, and the correlations were found between terahertz peaks and the X-ray signals. Since the intensity of the peak at 9.0 THz was correlated with crystallinity obtained using an X-ray diffraction (r² = 0.98), our data indicate that the Fourier-transform terahertz spectrometer can be a new analytical device to measure the starch crystallinity.

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.

Crystal structure of a new hydrate form of the NSAID sodium diclofenac

The crystal structure is reported of sodium 2-[2-(2,6-di-chloro-anilino)phen-yl]acetate 3.5-hydrate or tetra-μ-aqua-κ8 O:O-deca-aqua-bis-{μ3-2-[2-(2,6-di-chloro-anilino)phen-yl]acetato-κ3 O:O:O}tetra-sodium(I) bis-{2-[2-(2,6-di-chloro-anil-ino)phen-yl]acetate}, {[Na4(C14H10Cl2NO2)2(H2O)14](C14H10Cl2NO2)2} n , which re-presents a new hydrate form of the NSAID sodium diclofenac (SD). The triclinic unit cell contains one ionic compound with formula Na4(C14H10Cl2NO2)4(H2O)14, in which two symmetry-related carboxyl-ate anions C14H10Cl2NO2 - are bonded to a centrosymmetric [Na4]4+ core cationic cluster, while the others are only hydrogen bonded to the cationic cluster. The conformation for the anions is similar to that found in other diclofenac compounds, and the [Na4(Ocarbox)2(H2O)14]4+ cluster displays an unprecedented geometry, which can be described as an incomplete dicubane cluster formed by face-sharing incomplete cubes. A complex framework of O-H⋯O hydrogen bonds stabilizes the crystal structure. The herein reported crystal structure for SD·3.5H2O in space group P is different from those previously reported for other hydrates, namely SD·4.75H2O (P21) and SD·5H2O (P21/m).

Simulation and Assignment of the Terahertz Vibrational Spectra of Enalapril Maleate Cocrystal Polymorphs

The identification of crystalline drug polymorphs using terahertz vibrational spectroscopy is a powerful approach for the nondestructive and noninvasive characterization of solid-state pharmaceuticals. However, a complete understanding of the terahertz spectra of molecular solids is challenging to obtain because of the complex nature of the low-frequency vibrational motions found in the sub-3 THz (sub-100 cm^-1) range. Unambiguous assignments of the observed spectral features can be achieved through quantum mechanical solid-state simulations of crystal structures and lattice vibrations utilizing the periodic boundary condition approach. The terahertz spectra of two polymorphs of enalapril maleate are presented here to demonstrate that even large pharmaceuticals can be successfully modeled using solid-state density functional theory, including cocrystalline solids comprised of multiple distinct species. These simulations enable spectral assignments to be made, but also provide insights into the conformational and cohesion energies that contribute to the polymorph stabilities. The results reveal that the Form II polymorph of enalapril maleate is the more stable of the two under ambient conditions, and that this stability is driven by a greater intermolecular cohesion energy as compared to Form I.

Terahertz spectroscopic characterizations and DFT calculations of carbamazepine cocrystals with nicotinamide, saccharin and fumaric acid

Carbamazepine cocrystals with nicotinamide, saccharin and fumaric acid were synthesized and characterized by time-domain terahertz spectroscopy. Lattice vibrations of cocrystals with their individual constituents were investigated by means of the dispersion-corrected density functional theory with and without cell parameter constraints. The simulated THz spectra successfully reproduce the features of all the crystals in their experimental spectra. A better agreement between experimental and theoretical THz spectra is achieved when the cell parameter constraints are applied in geometry optimization. Some intensive modes of neat carbamazepine and cocrystals were discussed in terms of the motions of hydrogen bonds. The effect of lattice vibration on these cocrystallizations was further examined to gain insights into the thermodynamics. It is found that lattice vibration is favorable for all these cocrystal formations.

Combined X-ray Crystallographic, IR/Raman Spectroscopic, and Periodic DFT Investigations of New Multicomponent Crystalline Forms of Anthelmintic Drugs: A Case Study of Carbendazim Maleate

Synthesis of multicomponent solid forms is an important method of modifying and fine-tuning the most critical physicochemical properties of drug compounds. The design of new multicomponent pharmaceutical materials requires reliable information about the supramolecular arrangement of molecules and detailed description of the intermolecular interactions in the crystal structure. It implies the use of a combination of different experimental and theoretical investigation methods. Organic salts present new challenges for those who develop theoretical approaches describing the structure, spectral properties, and lattice energy Elatt. These crystals consist of closed-shell organic ions interacting through relatively strong hydrogen bonds, which leads to Elatt > 200 kJ/mol. Some technical problems that a user of periodic (solid-state) density functional theory (DFT) programs encounters when calculating the properties of these crystals still remain unsolved, for example, the influence of cell parameter optimization on the Elatt value, wave numbers, relative intensity of Raman-active vibrations in the low-frequency region, etc. In this work, various properties of a new two-component carbendazim maleate crystal were experimentally investigated, and the applicability of different DFT functionals and empirical Grimme corrections to the description of the obtained structural and spectroscopic properties was tested. Based on this, practical recommendations were developed for further theoretical studies of multicomponent organic pharmaceutical crystals.

Structural insights into anhydrous and monohydrated forms of 2,4,6-trihydroxybenzoic acid based on Raman and terahertz spectroscopic characterization

In order to characterize molecular structures of 2,4,6-trihydroxybenzoic acid (PCA) by means of vibrational spectroscopic techniques, we report investigation of PCA monohydrated form and its anhydrous polymorphic one by using terahertz and Raman spectral characterization. The experimental THz spectra show that the monohydrated PCA only has two absorption bands at 0.69 and 1.65 THz respectively in the frequency region from 0.2 to 1.8 THz, meanwhile the anhydrous form has a few significantly different absorption bands at 0.75, 1.01, 1.46 and 1.64 THz, respectively. Furthermore, Raman spectra characterized such differences of vibrational modes shown within 200-1800 cm^-1 region about the monohydrated and anhydrous forms of PCA. In view of various possible theoretical structural forms that may exist in anhydrous PCA and its monohydrated one, density functional theory calculations were performed to simulate optimized structures and vibrational mode of above two PCA polymorphic forms. Theoretical results and experimental THz/Raman spectra of anhydrous PCA show that the dimer synthon via the carboxylic group ••• carboxyl group and its ortho-phenolic hydroxyl group inter-molecular hydrogen bonding interaction establishing the theoretical form I (AH-I) is more consistent with experimental observation than other theoretical forms (AH-II and AH-III). Meanwhile, the theoretical monohydrated form I (MH-I), which is formed by the linkage of carboxyl group and its ortho-phenolic hydroxyl group with water molecule, is also much more agreement with experimental spectral observations of PCA monohydrate than other monohydrated forms (MH-II and MH-III). Our study demonstrates effectively qualitative analysis of both micro-molecular structures and dehydrated transitions between anhydrous and hydrated polymorphic forms of PCA, thus providing rich information on the corresponding structural changes of anhydrous and hydrated PCAs due to various inter-molecular and intra-molecular interactions based on their finger-print vibrational spectra combined with theoretical simulations.

Structure and spectroscopic characterization of pharmaceutical co-crystal formation between acetazolamide and 4-hydroxybenzoic acid

Co-crystals have great potential for drug research and development because the formation of co-crystal is accompanied by changes inter-molecular interactions between starting materials that enable to improve both physical and chemical properties of active pharmaceutical ingredients. In order to provide a more profound insight into the structural changes of specific drugs upon co-crystallization, spectroscopic characterization of solid-state acetazolamide (ACZ), 4-hydroxybenzoic acid (4HBA) and their co-crystal prepared by mechanical grinding approach has been performed with spectral techniques including terahertz time-domain spectroscopy (THz-TDS) and Raman spectroscopy. Experimental THz spectra show that the ACZ-4HBA co-crystal has a few significantly different absorption peaks in 0.82, 1.16, 1.28 and 1.64 THz respectively compared with parent materials in the frequency region from 0.2 to 1.8 THz. Likewise, such differences between the co-crystal and starting compounds could also be characterized by Raman vibrational spectra. Moreover, density functional theory (DFT) calculations were performed to simulate optimized structures and vibrational modes of three kind of possible co-crystal theoretical forms (form I, II and III) between ACZ and 4HBA. Theoretical results and THz/Raman vibrational spectra of ACZ-4HBA co-crystal show that the 4HBA links to the thiadiazole acetamide fragment of ACZ via the double-bridged heterodimeric synthon C(N)NH⋯HOOC inter-molecular hydrogen bonding interaction establishing the theoretical form I, which is more consistent with experimental observations than other two possible theoretical co-crystal forms. These results provide rich information and unique method for characterizing the composition of co-crystal structures and also inter-molecular interactions shown within pharmaceutical co-crystallization process at the molecular level.

Structural investigation of a 2:1 co-crystal between diflunisal and isonicotinamide based on terahertz and Raman spectroscopy

In order to characterize molecular structure changes of drugs upon co-crystallization by means of spectroscopic techniques, vibrational spectra of solid-state diflunisal (DIF), isonicotinamide (ISO) and their 2:1 co-crystal have been investigated by using terahertz time-domain spectroscopy (THz-TDS) and Raman spectroscopy. A 2:1 DIF-ISO co-crystal between DIF and ISO has been synthesized by slow solution crystallization from ethanol. The experimental THz spectroscopy shows that the co-crystal has a few significantly different absorption peaks compared with raw parent materials within the frequency region from 0.2 to 1.6 THz. Likewise, some differences of vibrational spectra between the co-crystal and starting compounds could also be characterized by Raman spectral results. Density functional theory (DFT) was used to simulate optimized structures and vibrational modes of two kind of possible co-crystal theoretical forms (form I and II) between DIF and ISO. Theoretical co-crystal form I is shown with 2:1 theoretical binary-adduct formed by carboxylic acid-amide and carboxylic acid-pyridine under inter-molecular hydrogen bonding. Theoretical co-crystal form II has a similar structure as form I, meanwhile the only difference is that O63 atom simultaneously forms hydrogen bond with H33 and H64. Also the hydroxyl -OH and carboxyl group -COOH establish molecular heterocycle under intra-molecular hydrogen bonds in both forms. The theoretical results show that both THz and Raman spectra of co-crystal form II between DIF and ISO is more consistent with the experimental observations than those of co-crystal form I. These results provide us with a wealth of information and unique method for characterizing the composition of co-crystal structures and also inter-molecular hydrogen bonding interactions shown within pharmaceutical co-crystallization at the molecular level.

Effects of Solvent Stabilization on Pharmaceutical Crystallization: Investigating Conformational Polymorphism of Probucol Using Combined Solid-State Density Functional Theory, Molecular Dynamics, and Terahertz Spectroscopy

Solid-state density functional theory (DFT), molecular dynamics (MD), and terahertz (THz) spectroscopy were used to study the formation of enantiotropically related conformational Form I and Form II polymorphs of the pharmaceutical compound, probucol. DFT calculations were performed on the crystal systems to compare relative lattice energies and the solvent stabilization of the metastable Form II structure. The thermodynamics of solvent inclusion in the Form II·MeOH crystal system were determined from MD simulations, as was the favored conformation of molecular probucol in methanol and ethanol solutions. The findings from both solid-state DFT and MD calculations suggest that the preferred molecular orientations of the probucol molecule in solution and the probable inclusion of methanol in the crystal lattice during the crystallization process lead to the solvent selectivity of the probucol polymorph formation. The additional stabilization energy provided by the crystallization solvent facilitates the nucleation and growth of the Form II polymorph under conditions that favor this metastable crystal form over the thermodynamically stable Form I, despite the higher energy molecular and crystalline configurations of probucol Form II. We demonstrate the influence of solvent on the formation of pharmaceutical polymorphs and provide a molecular-level view of complex interactions leading to polymorphism using a combination of computational methods and THz spectral data.

Optical Properties of Meloxicam in the Far-Infrared Spectral Region

One of the most commonly used nonsteroidal anti-inflammatory active pharmaceutical ingredient called Meloxicam has been characterized spectroscopically both by Terahertz (THz) time domain spectroscopy (THz-TDS) and by Fourier Transform Infrared (FTIR) spectroscopy in far-IR regions of electromagnetic spectrum; 0.2 THz to 20 THz. While many relatively sharp features are observed in the far-IR range between 2 THz to 20 THz as expected for being an organic substance, very distinct and relatively strong absorption bands are also observed at 1.00, 1.66, 2.07 and 2.57 THz in the THz range. These well separated, defined, and fairly strong spectral features can be used for discrimination and quantification of Meloxicam in drug analysis. Frequency dependent refractive index of the drug was determined in a range of 0.2 THz and 2.7 THz, where an almost constant index was observed with an average index of 1.75. Powder XRD, and solid-state Density Functional Theory (SS-DFT) calculations were utilized to determine the crystalline form of the Meloxicam sample in its enolic crystalline form. Single molecule DFT calculations were also performed in all four possible structures of Meloxicam. In addition, the capability of THz waves transmission through common packaging materials is demonstrated for possibility of future on-site analysis. The results suggest that drug analysis will be possible to perform not only at every stage of manufacturing without destruction but also directly at the shelf of a market after development of portable THz technologies.

Relationship between electron-phonon interaction and low-frequency Raman anisotropy in high-mobility organic semiconductors

Recent theoretical studies have shown that charge transport in high-mobility organic semiconductors is limited by low-frequency vibrations because of strong non-local electron-phonon interaction. Here we investigate two high-electron-mobility organic semiconductors with similar molecular structures but considerably different crystal packings, TCNQ and F2-TCNQ, and reveal the relationship between the experimental low-frequency Raman spectra and the calculated contributions of various vibrational modes to the electron-phonon interaction. We suggest that the combination of Raman spectroscopy with solid-state DFT is a powerful tool for probing electron-phonon interaction and focused search for high-mobility organic semiconductors.

A terahertz study of taurine: Dispersion correction and mode couplings

The low-frequency characteristics of polycrystalline taurine were studied experimentally by terahertz (THz) absorption spectroscopy and theoretically by ab initio density-functional simulations. Full optimizations with semi-empirical dispersion correction were performed in spectral computations and vibrational mode assignments. For comparison, partial optimizations with pure density functional theory were conducted in parallel. Results indicate that adding long-range dispersion correction to the standard DFT better reproduces the measured THz spectra than the popular partial optimizations. The main origins of the observed absorption features were also identified. Moreover, a coupled-oscillators model was proposed to explain the experimental observation of the unusual spectralblue-shift with the increase of temperature. Such coupled-oscillators model not only provides insights into the temperature dynamics of non-bonded interactions but also offers an opportunity to better understand the physical mechanisms behind the unusual THz spectral behaviors in taurine. Particularly, the simulation approach and novel coupled-oscillators model presented in this work are applicable to analyze the THz spectra of other molecular systems.

Solvation of diclofenac in water from atomistic molecular dynamics simulations - interplay between solute-solute and solute-solvent interactions

The solubility-permeability relationship of active pharmaceutical ingredients determines the efficacy of their usage. Diclofenac (DCL), which is a widely used nonsteroidal anti-inflammatory drug, is characterized by extremely good membrane permeability, but low water solubility limiting drug effectiveness. The present research focuses on the fundamental explanation of this limitation using the combination of ab initio and classical molecular dynamics simulations of different ionic forms of DCL in water, namely, ionized, un-ionized and the mixture of them both. The analysis of diclofenac solvation in an aqueous environment is used to understand the origin of drug precipitation, especially in gastric pH. The used computational approach reveals the formation of micelle-like self-associated aggregates of diclofenac in water as the result of intermolecular π-π interactions and C-Hπ hydrogen bonds. The DCL aggregation in water is shown to depend mostly on drug concentration, protonation and temperature of the aqueous environment. The detected self-association properties of the drug in water are likely to be of great importance during the development of new drug formulations and fabrication of drug adsorbents for wastewater.

The effect of the flexibility of hydrogen bonding network on low-frequency motions of amino acids. Evidence from Terahertz spectroscopy and DFT calculations

Low-frequency modes of L-Asp and L-Asn were studied in the range from 0.1 to 3.0THz using time-domain Terahertz spectroscopy and density functional theory calculation. The results show that PBE-D2 shows more success than BLYP-D2 in prediction of THz absorption spectra. To compare their low-frequency modes, we adopted "vibrational character ID strips" proposed by Schmuttenmaer and coworkers [Journal of Physical Chemistry B, 117, 10444(2013)]. We found that the most intense THz absorption peaks of two compounds both involve severe distortion of their hydrogen bonding networks. Due to less rigid hydrogen bonding network in L-Asp, the side chain (carboxyl group) of L-Asp exhibits larger motions than that (carboxamide group) of L-Asn in low-frequency modes.

Spectroscopic investigation on structure and pH dependent Cocrystal formation between gamma-aminobutyric acid and benzoic acid

Vibrational spectroscopic methods, including terahertz absorption and Raman scattering spectroscopy, were utilized for the characterization and analysis of gamma-aminobutyric acid (GABA), benzoic acid (BA), and the corresponding GABA-BA cocrystal formation under various pH values of aqueous solution. Vibrational spectroscopic results demonstrated that the solvent GABA-BA cocrystal, similar as grinding counterpart, possessed unique characteristic features compared with that of starting parent compounds. The change of vibrational modes for GABA-BA cocrystal comparing with starting components indicates there is strong inter-molecular interaction between GABA and BA molecules during its cocrystallization process. Formation of GABA-BA cocrystal under slow solvent evaporation is impacted by the pH value of aqueous solution. Vibrational spectra indicate that the GABA-BA cocrystal could be stably formed with the solvent condition of 2.00≤pH≤7.00. In contrast, such cocrystallization did not occur and the cocrystal would dissociate into its parent components when the pH value of solvent is lower than 2.00. This study provides experimental benchmark to discriminate and identify the structure of cocrystal and also pH-dependent cocrystallization effect with vibrational spectroscopic techniques in solid-state pharmaceutical fields.

Solid-state cocrystal formation between acyclovir and fumaric acid: Terahertz and Raman vibrational spectroscopic studies

The vibrational spectra of solid-state acyclovir, fumaric acid and their cocrystal have been investigated by using terahertz time-domain spectroscopy (THz-TDS) and Raman spectroscopy at room temperature. In experimental THz spectra, the cocrystal has absorption peaks in 0.65, 0.94 and 1.10THz respectively, while the raw materials are absolutely different in this region. Raman spectra also show similar results about differences between the cocrystal and raw materials. Density functional theory (DFT) was performed to simulate vibrational modes of different theoretical forms between acyclovir and fumaric acid. The calculation of theoretical THz spectra shows that O8C7N1H27 and the carboxyl group COOH establish a dimer theoretical cocrystal form by the hydrogen bonding effect, which makes contributions to the formation of absorption peaks in 0.70, 1.01 and 1.34THz, and agrees well with experimental observations. The theoretical Raman result also indicates that this dimer form matches with experimental results. The characteristic bands of the cocrystal between acyclovir and fumaric acid are also assigned based on the simulation results from the DFT calculation.

Multivariate quantification of mebendazole polymorphs by terahertz time domain spectroscopy (THZ-TDS)

An analytical method based on terahertz-time domain spectroscopy (THz-TDS) and PLS regression models to quantify mebendazole polymorphs in pharmaceutical raw material is presented.

Examination of l-Glutamic Acid Polymorphs by Solid-State Density Functional Theory and Terahertz Spectroscopy

The ability of l-glutamic acid to crystallize in two different forms has long been the subject of study due to its commercial importance. While a solvent-mediated phase transformation between the α and β polymorphs is the prevailing theory, recent reports indicate a thermal solid-solid transformation between the two may be possible. However, determining accurate thermodynamic stabilities of these crystals has been challenging. Here new low-temperature single-crystal X-ray diffraction data coupled to solid-state density functional theory simulations have enabled a detailed description to be achieved for the energetic parameters governing the stabilization of the two l-glutamic acid solids. The temperature-dependent Gibbs free-energy curves show that α-glutamic acid is the preferred form at low temperatures (<222 K) and the β form is most stable at ambient temperatures. Terahertz time-domain spectroscopy was utilized to evaluate the quality of the intermolecular force modeling as well as to provide characteristic low-frequency spectral data that can be used for quantification of polymorph mixtures or crystal growth monitoring.

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.

Uncovering the Terahertz Spectrum of Copper Sulfate Pentahydrate

Terahertz vibrational spectroscopy has evolved into a powerful tool for the detection and characterization of transition metal sulfate compounds, specifically for its ability to differentiate between various hydrated forms with high specificity. Copper(II) sulfate is one such system where multiple crystalline hydrates have had their terahertz spectra fully assigned, and the unique spectral fingerprints of the forms allows for characterization of multicomponent systems with relative ease. Yet the most commonly occurring form, copper(II) sulfate pentahydrate (CuSO4·5H2O), has proven elusive due to the presence of a broad absorption across much of the terahertz region, making the unambiguous identification of its spectral signature difficult. Here, it is shown that the sub-100 cm(-1) spectrum of CuSO4·5H2O is obscured by absorption from adsorbed water and that controlled drying reveals sharp underlying features. The crystalline composition of the samples was monitored in parallel by X-ray diffraction as a function of drying time, supporting the spectroscopic results. Finally, the terahertz spectrum of CuSO4·5H2O was fully assigned using solid-state density functional theory simulations, helping attribute the additional absorptions that appear after excessive drying to formation of CuSO4·3H2O.

Component spectra extraction from terahertz measurements of unknown mixtures

The aim of this work is to extract component spectra from unknown mixtures in the terahertz region. To that end, a method, hard modeling factor analysis (HMFA), was applied to resolve terahertz spectral matrices collected from the unknown mixtures. This method does not require any expertise of the user and allows the consideration of nonlinear effects such as peak variations or peak shifts. It describes the spectra using a peak-based nonlinear mathematic model and builds the component spectra automatically by recombination of the resolved peaks through correlation analysis. Meanwhile, modifications on the method were made to take the features of terahertz spectra into account and to deal with the artificial baseline problem that troubles the extraction process of some terahertz spectra. In order to validate the proposed method, simulated wideband terahertz spectra of binary and ternary systems and experimental terahertz absorption spectra of amino acids mixtures were tested. In each test, not only the number of pure components could be correctly predicted but also the identified pure spectra had a good similarity with the true spectra. Moreover, the proposed method associated the molecular motions with the component extraction, making the identification process more physically meaningful and interpretable compared to other methods. The results indicate that the HMFA method with the modifications can be a practical tool for identifying component terahertz spectra in completely unknown mixtures. This work reports the solution to this kind of problem in the terahertz region for the first time, to the best of the authors' knowledge, and represents a significant advance toward exploring physical or chemical mechanisms of unknown complex systems by terahertz spectroscopy.