Next Article in Journal
A Fully Human Monoclonal Antibody Targeting cKIT Is a Potent Inhibitor of Pathological Choroidal Neovascularization in Mice
Next Article in Special Issue
Novel Polymorphic Cocrystals of the Non-Steroidal Anti-Inflammatory Drug Niflumic Acid: Expanding the Pharmaceutical Landscape
Previous Article in Journal
Probucol Pharmacological and Bio-Nanotechnological Effects on Surgically Transplanted Graft Due to Powerful Anti-Inflammatory, Anti-Fibrotic and Potential Bile Acid Modulatory Actions
Previous Article in Special Issue
Obtaining Cocrystals by Reaction Crystallization Method: Pharmaceutical Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Centre for THz Research, China Jiliang University, Hangzhou 310018, China
2
Department of Chemistry, Zhejiang Sci-Tech University, Hangzhou 310018, China
*
Author to whom correspondence should be addressed.
Submission received: 7 July 2021 / Revised: 15 August 2021 / Accepted: 16 August 2021 / Published: 20 August 2021
(This article belongs to the Special Issue Cocrystal Applications in Drug Delivery (Volume II))

Abstract

:
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.

Graphical Abstract

1. Introduction

Tuberculosis (TB) is a global emergency case announced by the World Health Organization (WHO) and among the top 10 causes of death [1,2]. The International Union against TB and Lung Diseases has recommended the use of fixed-dose combinations (FDCs) of anti-TB drugs to combat monotherapy and drug resistance [3,4]. Isoniazid (INH, as shown in Figure 1), the principal component in the so-called ‘triple therapy’ and different FDCs [5], tends to degrade in the FDCs involving INH and rifampicin (RIF). Moreover, oxidative stress will occur in the patients treated with anti-TB drugs such as INH and RIF, which tend to cause hepatotoxicity [6]. Therefore, in order to improve its pharmaceutical performance and enhance the stability of tablet formulations, it is urgent to explore new solid forms [7].
Cocrystal is generally defined as the crystalline single-phase materials [8] formed by two or more solid-state neutral compounds in a certain stoichiometric ratio under the influence of inter-molecular interactions [9]. Compared with cocrystal, salt formation is limited to acid–base reaction. However, cocrystallization provides different avenues, where any pharmaceutical substance may form cocrystals without having ionizable groups [10]. Additionally, on the premise of almost no effect on the chemical structure, intrinsic biological activity and pharmacological activity of APIs [11], cocrystallization strategy can significantly optimize the physicochemical properties and clinical performance of specific APIs [12,13]. Protocatechuic acid (PA, whose molecular structure is shown in Figure 2), as a typically natural antioxidant and antibacterial agent, has been used in pharmacology and food industry because of its antioxidant activity by eliminating free radicals via donating H-atom [14]. Additionally, containing the special hydroxyl-benzoic structures, as the naturally occurring phenolic acids [15], PA possesses a great deal of hydrogen bonding sites including the carboxylic acid group and the hydroxyl groups used as both hydrogen bond donor and hydrogen bond acceptor. Therefore, PA is an ideal co-former that is prone to form steady inter-molecular hydrogen bonds with other molecules.
In the past few years, many new solid forms involving INH have been studied, especially the cocrystallization of INH and different co-forms. Pawel Grobelny et al. [16] firstly reported the cocrystal of the anti-tuberculosis drugs INH and pyrazinamide at a 1:1 molar ratio by using X-ray diffraction (XRD), differential scanning calorimetry (DSC) and thermal gravimetric analysis (TGA) techniques, which included two symmetrical independent COOH…Npyridine hydrogen bonds. Liu Fang et al. [17] synthesized and characterized the cocrystal of INH and syringic acid using the technology of XRD and DSC, presenting a new opportunity for effectively reducing the hepatotoxicity of INH. In order to solve the problem of oxidative stress and low stability in tablet formulations involving INH, Mashhadi et al. prepared and characterized the INH-antioxidant hydroxybenzoic acid [18] cocrystals and INH-PA cocrystal [19] by the technology of infra-red (IR) spectroscopy, XRD and DSC. In addition, there are certain studies on the cocrystals regarding polyphenols [20] and aromatic carboxylic acids [21] as co-formers.
Most of the studies of cocrystals relating to INH above are implemented by a few traditional ways such as differential scanning calorimetry (DSC), X-ray diffraction (XRD), which are time-consuming in some cases or not detailed enough about the structures of molecules and the influence of intra-molecular and/or inter-molecular interactions within the cocrystals [22]. Thus, the vibrational spectroscopy including Fourier transformation infrared (FTIR), Raman and the emerging terahertz (THz) spectroscopy with high sensitivity and selectivity [23,24], providing a wealth of information about the structure of molecules, intra-molecular and/or inter-molecular interactions [25], is considered to be beneficial to solid-state characterization. Especially, THz spectroscopy can provide the finger-print spectra referring to special molecular structure [26,27] and much information about the weak collective molecular vibrations and various crystalline structural changes in the low frequency region, especially below 3.0 THz [28,29].
In the present work, the monohydrate cocrystal between INH and PA was synthesized by slow solvent evaporation and THz and Raman spectral techniques were used to observe the fingerprint information from INH, PA, their physical mixture and cocrystal. The terahertz spectrum is composed of a mixture of hydrogens, weak covalent bonds, ionic bonds and van der Waals force in the solid states [30]. These complex and highly entangled molecular interactions make the terahertz spectrum results difficult to explain, so it needs relatively complex calculation and simulation [31]. Density functional theory (DFT), as a very useful means of quantum mechanics calculation, is able to explain and assign in detail the complex molecular motions obtained from THz experiments and understand the basic origin of material properties [32]. Combined with DFT calculations, terahertz spectroscopy can be used to understand the vibrational motions in ordered crystal structures [1]. Comparing the experimental results with the theoretical results, the structural changes of molecules and the interactions between and within molecules could be better comprehended in the cocrystal formation process.

2. Materials and Methods

2.1. Material Preparation and the Synthesis of Cocrystal

Anhydrous solid-state INH (purity 98%) and sample PA (purity 99%) were acquired by J&K Chemical Company (Shanghai, China). Methanol was purchased from Hangzhou Gaojing Fine Chemical Company (Hangzhou, China). All chemicals were used as received without further purification.
The physical mixture of parent materials was obtained by mixing the two compounds at a molar ratio of 1:1 in a glass vial by using a vortex mixer. The cocrystal was obtained by slow solvent crystallization, which was commonly used in the preparation of cocrystal. INH (0.137 g, one mmol) and PA (0.154 g, one mmol) were respectively dissolved in a mixture of 20 mL methanol and de-ionized water (1:1 volume ratio; details about the generation and characterization of de-ionized water used could be found in materials of the supporting information), and then mixed together. The mixed solution was heated to 70 °C for 10 min and then held in slow evaporation for 15 days. The cocrystal was obtained by isolating the orange prism like crystals at room temperature. The separated cocrystal weighing 0.25 g was taken out and gently ground to make its average particle size several microns, so as to minimize the scattering effect of sample particles in terahertz spectrum measurement. Then the cocrystal substance was put into the tablet press to impose 4 MPa pressure for about two minutes, and the physical mixture performed the same operation. The sample prepared according to the above steps was a circular plate with a diameter of 13 mm and a thickness of 2 mm, which was put into a sealed bag for subsequent THz detection. For the detection of Raman spectrum, only a small amount of powder sample was needed without further sample preparation.

2.2. THz/Raman Device and Detection

THz spectra of INH, PA, their physical mixture and cocrystal were detected by THz-TDS Z2 system (Zomega Co., New York, NY, USA), which was applied to investigate the THz transmission spectra of target substances based on photoconductive switches. The optical pulse sequence was from Ti: Sapphire oscillator ultrafast laser pulse system (Spectra Physics, Owen, CA, USA) with the frequency 75 MHz, the pulse width 100 fs and the center wavelength 780 nm. Due to the strong absorption of THz radiation by water, nitrogen gas with high purity was continuously flushed into the detector during the experiment. Three measurements for each set of sample and reference were carried out and took the average value as the ultimate spectral data. Since the outputs of the experiment were time domain spectra of the electric field intensity of THz, they need to be converted into corresponding frequency domain spectra by fast Fourier transform (FFT) process, and the THz absorption spectra were acquired by dividing the sample frequency response by that of the reference. The Raman spectra of the above substances were detected by the Fourier transform Raman (FT-Raman) spectrometer (Spectra Physics, Owen, CA, USA) after over 500 scans at 2 cm−1 spectral resolution over the wavenumber range between 150 and 3500 cm−1. A diode pumped (wavelength at 1064 nm) solid-state laser was selected as the near-infrared light source with 150 mW operating power and the total analysis time for each sample was of the order of 5 min intervals.

2.3. DFT Theoretical Calculations

The crystal engineering method of supramolecular heterosynthesis between the API and the coformer is an effective strategy for the rational selection of suitable coformers and the design of new pharmaceutical cocrystals. It is important to understand the hydrogen bond interactions, which lead to strong heterosynthons in the cocrystal compared to homosynthons in individual components [33]. Therefore, it is necessary to theoretically predict and analyze the terahertz spectra of crystal materials. On the one hand, the less accurate semi-empirical method that allows full crystal calculations of the terahertz vibrations, not limited to one construction unit can reduce the calculation time [34]. In addition, when the correction factors are used to compensate the systematic deviations, ab initio and DFT methods are considered to be very accurate in predicting molecular vibrations. In this paper, Gaussian 09 software was used to simulate and optimize the theoretical structures of INH, PA and their possible cocrystal forms, and B3LYP was selected as a credible DFT geometry majorization method in the computational process [35,36]. Based on the structural characteristics of INH and PA, three possible cocrystal forms were simulated. In form I, the hydrogen bond was formed by pyridine N–carboxylic acid synthons of two raw materials (as shown in Figure 3a). In form II and form III, hydrogen bonds were formed between raw materials and single water molecule, composing different ring structures (as shown in Figure 3b,c). Although Mashhadi et al. [19] have synthesized and characterized the cocrystal of INH and PA by means of XRD and DSC, there is no structural report on the analysis of the cocrystal of the above substances based on THz and Raman spectra. To further characterize the hydrogen bond motifs shown within the INH-PA hydrate cocrystal structure, the results of THz and Raman spectra were combined with that of DFT simulation, which is of great significance to better comprehend the structural characteristics and vibrational modes of INH-PA hydrate cocrystal. In the simulated THz and Raman spectra, a full-width half-maximum (FWHM) value of 4.0 cm−1 was used to convolute the Lorentzian line shapes into the calculated vibrational modes.

3. Results and Discussion

3.1. THz Spectral Characterization and Analysis

The low frequency vibrational spectra of molecules provide many useful values for molecular vibration, rotation modes, and some weak inter-molecular interactions such as hydrogen bond interactions, van der Waals forces and so on. THz spectroscopy can be used to obtain such low frequency vibrational information of crystalline compounds, which makes it a vigorous analytical tool for studying the structure and inter-molecular interaction of various molecules [37]. The THz absorption spectrum of INH, PA, their physical mixture and corresponding INH-PA monohydrate cocrystal in the range of 0.2–1.6 THz is shown in Figure 4. It could be clearly seen that INH has two distinct characteristic peaks at 1.18 and 1.46 THz, which makes it possess good spectral resolution. Additionally, PA engendered three relatively flat absorption peaks at 1.20, 1.40 and 1.49 THz. In addition, the physical mixture mainly has three obvious absorption peaks, which are located at 1.18, 1.40 and 1.46 THz, respectively. There is no new characteristic peak, but a simple superposition of respective absorption spectra of the two original substances in the THz absorption spectrum of the physical mixture. Nevertheless, the THz absorption spectrum of INH-PA monohydrate cocrystal shows a few characteristic peaks at 0.88, 1.22, 1.46 and 1.53 THz differing from both two individual raw substances and their simple physical mixture. Therefore, above-mentioned characteristic peaks indicate that INH and PA have successfully synthesized new substance by slow evaporation under the influence of inter-molecular and/or intra-molecular interactions. This also verified the ability of THz spectroscopy to supply obvious fingerprint information for the molecular structure of solid-state crystal. Moreover, the cocrystal was characterized by the H-NMR spectrum as shown in Figure S1, which proved that the compound did contain water but did not contain other solvents.
Three possible theoretical cocrystal forms were simulated in this paper. In form I, the hydrogen bond was formed by pyridine N–carboxylic acid synthons of two raw materials (as shown in Figure 3a). In form II and form III, hydrogen bonds were formed between raw materials and single water molecule (as shown in Figure 3b,c). The comparisons of THz absorption spectra between experimental and theoretical results are shown in Figure 5. There are some obvious characteristic peaks at 0.4, 0.85 and 1.50 THz for theoretical form I and 0.46, 0.65, 1.10, 1.14 and 1.45 THz for theoretical form II, respectively. Similarly, theoretical form III reveals characteristic peaks with high peak intensity and resolution at 0.48, 0.86, 1.16 and 1.53 THz, respectively. However, the theoretical form I lacks the characteristic peak at 1.20 THz and the theoretical form II does not appear the corresponding characteristic peak at 0.88 THz, compared with the experimental absorption spectra. Surprisingly, the absorption peaks of theoretical form III at 0.86, 1.16 and 1.53 THz are in good agreement with the experimental characteristic peaks at 0.88, 1.20 and 1.53 THz. Therefore, the THz spectra of the theoretical cocrystal form III is more consistent with that of experimental results, which also proves the formation of cocrystal. Although, the absorption peak of theoretical cocrystal form III is unable to match the experiment at 0.48 and 1.46 THz and a certain red shift exists at 0.86 and 1.16 THz, which may result from the fact that the theoretical calculation is carried out at absolute zero temperature while the actual THz spectral detection is measured at room temperature [38]. Additionally, the theoretical simulation only calculates one construction unit of the whole cocrystal molecular network, while the experimental cocrystal spectrum is used to characterize the whole solid-state cocrystal molecular network. In addition, different computing settings will have a certain impact on consistency between theory and experiment as well.
GaussView is one kind of software specially designed for Gaussian, whose main purpose is to construct Gaussian input file and to display Gaussian output calculation results in the form of a graph. Through the GaussView software, the various vibrational modes of theoretical cocrystal form III between INH and PA (shown in Figure 6) could be conveniently observed and the vibrational modes at different characteristic peaks are assigned, as shown in Table 1. Corresponding to the absorption peak of experimental cocrystal at 0.88 THz, the characteristic peak of theoretical cocrystal form III at 0.86 THz is mainly caused by the out of plane rocking vibrations of INH and PA molecules. The out of plane rocking vibration of INH and in-plane rocking vibration of PA make the theoretical cocrystal form III format an obvious characteristic peak at 1.16 THz, which is consistent with the absorption peak of the experimental cocrystal at 1.20 THz. The crystal vibration of the absorption peaks at 1.53 THz of the experimental and theoretical cocrystal is mainly owing to the out of plane rocking vibration from the ring of INH and the torsional vibrations from O11-C12-N13-H17 of INH and H30-O29-C27-O28 of PA. Moreover, the above three vibrational modes are accompanied by the vibration of water molecules (H36-O35-H37), which is consistent with the theoretical simulation of INH, PA and water molecule forming O35-H36…O33, O35-H37…O11 and O33-H34…N14 inter-molecular hydrogen bonds. Through the analysis of vibrational modes, it can be proved that the INH-PA monohydrate cocrystal is formed chiefly due to the generation of inter-molecular hydrogen bonds, which plays a key role in changing the molecular structure and makes the vibrational modes different from those of the original parent materials. Therefore, THz spectroscopy can be applied to characterize the hydrogen bond effect sensitively, owing to its absorption mechanism affected by intra-molecular and inter-molecular vibrational excitation [39]. In addition, the powder diffraction experiments of INH, PA and their hydrate cocrystal were performed, and the powder diffraction pattern of INH-PA cocrystal was simulated according to the crystal structure data from the Cambridge Crystallographic Data Centre (CCDC). The experimental X-ray diffraction (XRD) pattern is almost consistent with the calculated pattern, which further confirms the formation of INH-PA monohydrate cocrystal (PXRD pattern is displayed as shown in Figure S2 in Supporting Information).

3.2. Raman Spectral Characterization and Analysis

Raman spectroscopy has become one of the preferred methods to study the solid-state polymorphism of drugs (polycrystalline, solvent, hydrate, salt, eutectic, etc.), not only because of its sensitivity to structural changes, but also its rapidity, non-destruction, no sample preparation and capability of characterization in the water environment [40]. Different from THz spectroscopy based on absorption mechanism, Raman spectroscopy mainly relies on light scattering to obtain an abundance of structure information from target substances, such as the vibrational transitions in molecules, electronic polarizabilities, and vibrational energy transitions of the molecules [41,42]. Figure 7 shows the Raman vibrational spectra of INH, PA, their physical mixture and INH-PA monohydrate cocrystal acquired by experiments. As it can be seen, the characteristic peaks of the physical mixture are mainly attributed to the simple linear superposition of the two raw materials, so it does not involve intra-molecular and/or inter-molecular interactions such as hydrogen bonding just during the simple mixing process, which is in good agreement with the above THz spectral results. Compared with the Raman spectra of INH-PA physical mixture, some corresponding characteristic peaks were not discovered at 364 and 798 cm−1 in the Raman spectra of INH-PA monohydrate cocrystal. Additionally, there is a new characteristic peak at 1654 cm−1 in Raman spectra of INH-PA monohydrate cocrystal, which are not found in the original parent materials and physical mixture. For the convenience of observation, the above-mentioned two kinds of characteristic peaks are covered with cyan shadows in Figure 7. In addition, corresponding to the characteristic peaks of physical mixture at 680, 755, 1002 and 1309 cm−1, the relevant characteristic peaks of cocrystal have a certain blue shift and/or red shift, respectively, located at 699, 773, 1023 and 1319 cm−1. Such characteristic peaks are shadowed with light green color in Figure 7. According to the analysis of the above spectral diversity, it can be inferred that a new substance between INH and PA is produced, and the characteristic peaks of the substance which reflect the molecular structure are completely different from those of the original materials in some specific locations.
The Raman spectra of the experimental results compared with the theoretical calculation are shown in Figure 8 (the theoretical cocrystal consists of three possible forms) and some differences among the three theoretical spectra are observed. Corresponding to the characteristic peak at 1184 cm−1 in the experimental Raman spectrum, theoretical form II and theoretical form III are in good agreement, while theoretical form I lacks the relevant characteristic peak. Comparing the theoretical form II with the theoretical form III, the characteristic peak of the theoretical form III at 1442 cm−1 is well matched with that of experiment, but the theoretical form II is short of such characteristic peak. Additionally, located in the high frequency region of 1610 cm−1, the number of characteristic peaks of form III is more consistent with the experimental spectrum. Therefore, the Raman spectrum from theoretical form III is in better agreement with the experimental results, which is consistent with the above THz spectrum results and proves the formation of cocrystal. In addition, the distinct vibrational modes of INH-PA monohydrate cocrystal are assigned in Table 2. Combined with the theoretical simulation results, it can be obtained that the characteristic peak at 798 cm−1 is mainly attributed to the scissoring vibration of O33-C22-C21-H26, deformation vibration of PA ring and the out of plane rocking vibration of H36-O35. Another characteristic peak with high resolution at 1318 cm−1 is mainly formed by the in-plane rocking vibration of H17-N13-C12, the out of plane rocking vibration of H15-N14-H16 and stretching vibration of N14-N13. In addition, two obvious characteristic peaks located in the high frequency region are accompanied by the vibration of the water molecule during the crystal vibration. For instance, the characteristic peak at 1632 cm−1 can be considered as the result of deformation vibration of INH ring and the scissoring vibration of water (H36-O35-H37). Additionally, the characteristic peak at 1655 cm−1 can be attributed to the deformation vibration of PA ring and the scissoring vibration of H34-O33-C22 and water (H36-O35-H37). Among them, O35-H36... O33 and O35-H37... O11 are the two important hydrogen bonds forming the cocrystal structure. By analyzing the vibrational modes of the four Raman characteristic peaks with high resolution, it can be seen that most of the vibrations of INH-PA hydrate cocrystal are influenced by inter-molecular interactions, especially the effect of inter-molecular hydrogen bonds between INH, PA and water molecule. Therefore, the hydrogen bonds play an important role in the formation of INH-PA monohydrate cocrystal.
It is well known that the change of bond length between different atoms is the most direct reason for the characteristic vibration mode information exhibited in the Raman spectra of various compounds [41,42]. Therefore, to further understand the specific micro-molecular structure of INH and PA monohydrate cocrystal, the research on the typical bond length change of cocrystal structure and the relationship with corresponding characteristic peaks were also added. Typical bond lengths of INH-PA monohydrate cocrystal are shown in Figure 9 and some typical bond length variations upon the formation of cocrystal are also given in Table 3. The bond lengths of O11-C12, N13-N14, N14-H15 and N14-H16 from INH changed by 0.251, 0.071, 0.016 and 0.022 Å, respectively, while H34-O33 and O33-C22 from PA changed by 0.007 and 0.061 Å, respectively. Moreover, the bond lengths of H37-O35 and O35-H36 in water molecule also increased. The change of the length of above-mentioned chemical bond is closely related to the two inter-molecular hydrogen bonds (O35-H37…O11 and O35-H36…O33) forming INH-PA monohydrate cocrystal. In addition, it can be seen from Table 3 that the chemical bonds with the largest bond length change (including O11-C12 of INH and O33-C22 of PA) form inter-molecular hydrogen bonds with water molecule, which greatly contribute to the vibrational modes at 798, 1318, 1632 and 1655 cm−1 in the above Raman spectra, which also confirms the important influence of inter-molecular hydrogen bonding on the formation of cocrystal. Combining the results of chemical bond length variations with THz and Raman spectra, it can help us better understand the offset, emergence and vanishing of characteristic peaks in the spectrum of INH-PA monohydrate cocrystal and provides us with abundant information about the molecular structure and inter-molecular interactions (particularly the hydrogen bonds between raw substances).

4. Conclusions

In this paper, the vibrational spectra of solid-state INH, PA, physical mixtures and INH-PA monohydrate cocrystal were obtained by THz and Raman spectra. It can be perceived that there are some differences between the vibrational spectra of INH-PA monohydrate cocrystal and those of original materials. Three possible theoretical cocrystal forms were simulated by DFT so as to compare with the experimental spectra, then the vibrational modes at different absorption peaks were assigned according to the above simulation results. According to the above analysis results, it can be concluded that the monohydrate cocrystal between INH and PA is formed on the score of non-covalent bond interactions between molecules, particularly the influence from inter-molecular hydrogen bonds. Through the comparison of three theoretical cocrystal forms, the theoretical cocrystal form III (INH, PA and water form O35-H37…O11 and O35-H36…O33 inter-molecular hydrogen bonding) is more consistent with the experimental results. In addition, the study on typical hydrogen bond length also confirmed the important role of hydrogen bonds in the formation of cocrystal. Raman and terahertz spectroscopic techniques are complementary to the investigation into vibrational modes and micro-molecular structure of INH-PA monohydrate cocrystal. The study shows the superiority of vibrational spectroscopy adequately certifying the molecular structure and inter-molecular hydrogen bond interactions in the study of specific cocrystals, which provides theoretical and experimental basis for the structural research of emerging pharmaceutical cocrystal fields at the micro-molecular level.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/pharmaceutics13081303/s1, Figure S1: H-NMR spectrum of INH-PA monohydrate cocrystal in the region of 0.5–8.5 ppm. The signal at 2.51 ppm is DMSO from the solvent., Figure S2: PXRD comparison of INH (a), PA (b), INH-PA monohydrate cocrystal from experiment (c) and INH-PA monohydrate cocrystal from theoretical simulation (d) in the 2θ range 5–50°, Figure S3: The THz spectra of INH (a), PA (b), their physical mixture (c) and INH-PA monohydrate cocrystal (d) in 0.6–1.6 THz spectral region. Experimental Details and Characterization Techniques.

Author Contributions

Data curation, Y.B., J.F., Z.Z. and J.L.; formal analysis, Y.B., J.F. and J.L.; funding acquisition, J.X., Z.H. and Y.D.; investigation, Z.Z., J.X., Z.H. and Y.D.; methodology, Y.B., J.X., J.L. and Z.H.; project administration, Y.D.; software, Y.B., Z.Z. and J.X.; supervision, Z.H. and Y.D.; visualization, J.F. and Z.Z.; writing—original draft, Y.B., J.F. and Z.Z.; writing—review and editing, Z.H. and Y.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partly funded both Natural Science Foundation of Zhejiang Province (NSFZJ, Grant No. LY19B050003) and National Natural Science Foundation of China (NSFC, Grant Nos. 21973082, 61875179).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Neu, J.; Nikonow, H.; Schmuttenmaer, C.A. Terahertz Spectroscopy and Density Functional Theory Calculations of dl -Norleucine and dl -Methionine. J. Phys. Chem. A 2018, 122, 5978–5982. [Google Scholar] [CrossRef]
  2. Babu, N.J.; Nangia, A. Solubility advantage of amorphous drugs and pharmaceutical cocrystals. Cryst. Growth Des. 2011, 11, 2662–2679. [Google Scholar] [CrossRef]
  3. Carvalho, S.A.; da Silva, E.F.; de Souza, M.V.N.; Lourenço, M.C.S.; Vicente, F.R. Synthesis and antimycobacterial evaluation of new trans-cinnamic acid hydrazide derivatives. Bioorganic Med. Chem. Lett. 2008, 18, 538–541. [Google Scholar] [CrossRef] [PubMed]
  4. Bhutani, H.; Singh, S.; Jindal, K.C. Drug-drug interaction studies on first-line anti-tuberculosis drugs. Pharm. Dev. Technol. 2005, 10, 517–524. [Google Scholar] [CrossRef] [PubMed]
  5. Iseman, M.D. Tuberculosis therapy: Past, present and future. Eur. Respir. J. Suppl. 2002, 20, 87–94. [Google Scholar] [CrossRef] [Green Version]
  6. Singh, M.; Sasi, P.; Gupta, V.H.; Rai, G.; Amarapurkar, D.N.; Wangikar, P.P. Protective effect of curcumin, silymarin and N -acetylcysteine on antitubercular drug-induced hepatotoxicity assessed in an in vitro model. Hum. Exp. Toxicol. 2012, 31, 788–797. [Google Scholar] [CrossRef] [PubMed]
  7. Wang, Y.; Xue, J.; Wang, Q.; Jin, S.; Zhang, Z.; Hong, Z.; Yong, D. Structural investigation of a 2:1 co-crystal between diflunisal and isonicotinamide based on terahertz and Raman spectroscopy. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2019, 216, 98–104. [Google Scholar] [CrossRef]
  8. Wang, Y.; Xue, J.; Qin, J.; Liu, J.; Yong, D. Structure and spectroscopic characterization of pharmaceutical co-crystal formation between acetazolamide and 4-hydroxybenzoic acid. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2019, 219, 419–426. [Google Scholar] [CrossRef]
  9. Vishweshwar, P.; McMahon, J.A.; Bis, J.A.; Zaworotko, M.J. Pharmaceutical Co-Crystals. J. Pharm. Sci. 2006, 95, 499–516. [Google Scholar] [CrossRef]
  10. Zhang, Z.; Cai, Q.; Xue, J.; Qin, J.; Liu, J.; Du, Y. Co-crystal formation of antibiotic nitrofurantoin drug and melamine co-former based on a vibrational spectroscopic study. Pharmaceutics 2019, 11, 56. [Google Scholar] [CrossRef] [Green Version]
  11. Korotkova, E.I.; Kratochvll, B. Pharmaceutical Cocrystals. Procedia Chem. 2014, 10, 473–476. [Google Scholar] [CrossRef] [Green Version]
  12. Zhixiong, H.; Xuesong, G.; Huanzhong, C. Study on biosensors based on terahertz time—Domain spectroscopy. Electron. Test. 2019, 20, 24–26. [Google Scholar] [CrossRef]
  13. Liu, F.; Wang, L.; Yu, M.; Li, Y.; Wu, Z.; Yan, C. A new cocrystal of isoniazid-quercetin with hepatoprotective effect: The design, structure, and in vitro/in vivo performance evaluation. Eur. J. Pharm. Sci. 2020, 144, 105216. [Google Scholar] [CrossRef]
  14. Chen, S.; Wang, X.; Chen, D. Antioxidant Activity and Mechanism of Protocatechuic Acid in vitro. Funct. Foods Health Dis. 2011, 1, 232–244. [Google Scholar] [CrossRef]
  15. Wei, J.; Liu, K.; Du, C. A Novel Mannich Derivative of Protocatechuic Acid: Synthesis, Crystal Structure and Antioxidant Activity. Proc. Natl. Acad. Sci. India Sect. A Phys. Sci. 2017, 87, 181–188. [Google Scholar] [CrossRef]
  16. Grobelny, P.; Mukhrrjee, A.; Deiraju, G.R. Drug-drug co-crystals: Temperature-dependent proton mobility in the molecular complex of isoniazid with 4-aminosalicylic acid. CrystEngComm 2011, 13, 4358–4364. [Google Scholar] [CrossRef]
  17. Liu, F.; Jiang, F.B.; Li, Y.T.; Liu, R.M.; Wu, Z.Y.; Yan, C.W. Cocrystallization with syringic acid presents a new opportunity for effectively reducing the hepatotoxicity of isoniazid. Drug Dev. Ind. Pharm. 2020, 46, 988–995. [Google Scholar] [CrossRef]
  18. Mashhadi, S.M.A.; Yunus, U.; Bhatti, M.H.; Tahir, M.N. Isoniazid cocrystals with anti-oxidant hydroxy benzoic acids. J. Mol. Struct. 2014, 1076, 446–452. [Google Scholar] [CrossRef]
  19. Mashhadi, S.M.A.; Yunus, U.; Bhatti, M.H.; Ahmed, I.; Tahir, M.N. Synthesis, characterization, solubility and stability studies of hydrate cocrystal of antitubercular Isoniazid with antioxidant and anti-bacterial Protocatechuic acid. J. Mol. Struct. 2016, 1117, 17–21. [Google Scholar] [CrossRef]
  20. González-González, J.S.; Martínez-Santiago, A.M.M.; Martínez-Martínez, F.J. Cocrystals of Isoniazid with Polyphenols: Mechanochemical Synthesis and Molecular Structure. Crystals 2020, 10, 569. [Google Scholar] [CrossRef]
  21. Diniz, L.F.; Souza, M.S.; Carvalho, P.S.; da Silva, C.C.P.; D’Vries, R.F. Novel Isoniazid cocrystals with aromatic carboxylic acids: Crystal engineering, spectroscopy and thermochemical investigations. J. Mol. Struct. 2018, 1153, 58–68. [Google Scholar] [CrossRef]
  22. Zhang, Z.; Fang, J.; Bo, Y.; Xue, J.; Liu, J.; Hong, Z.; Du, Y. Terahertz and Raman Spectroscopic Investigation of Anti-tuberculosis Drug-Drug Cocrystallization Involving 4-aminosalicylic Acid and Pyrazinamide. J. Mol. Struct. 2021, 1227. [Google Scholar] [CrossRef]
  23. Ayala, A.P.; Caetano, M.W.C.; Honorato, S.B.; Mendes Filho, J.; Siesler, H.W.; Faudone, S.N.; Cuffini, S.L.; Martins, F.T.; Da Silva, C.C.P.; Ellena, J. Conformational polymorphism of the antidiabetic drug chlorpropamide. J. Raman Spectrosc. 2012, 43, 263–272. [Google Scholar] [CrossRef]
  24. King, M.D.; Buchanan, W.D.; Korter, T.M. Identification and Quantification of Polymorphism in the Pharmaceutical Compound Diclofenac Acid by Terahertz Spectroscopy and Solid-State Density Functional Theory. Anal. Chem. 2011, 83, 3786–3792. [Google Scholar] [CrossRef] [PubMed]
  25. Zhang, G.C.; Lin, H.L.; Lin, S.Y. Thermal analysis and FTIR spectral curve-fitting investigation of formation mechanism and stability of indomethacin-saccharin cocrystals via solid-state grinding process. J. Pharm. Biomed. Anal. 2012, 66, 162–169. [Google Scholar] [CrossRef]
  26. Du, Y.; Fang, H.X.; Zhang, Q.; Zhang, H.L.; Hong, Z. Spectroscopic investigation on cocrystal formation between adenine and fumaric acid based on infrared and Raman techniques. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2016, 153, 580–585. [Google Scholar] [CrossRef]
  27. Du, Y.; Xia, Y.; Zhang, H.; Hong, Z. Using terahertz time-domain spectroscopical technique to monitor cocrystal formation between piracetam and 2,5-dihydroxybenzoic acid. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2013, 111, 192–195. [Google Scholar] [CrossRef]
  28. Du, Y.; Wang, Y.; Xue, J.; Liu, J.; Qin, J.; Hong, Z. Structural insights into anhydrous and monohydrated forms of 2,4,6-trihydroxybenzoic acid based on Raman and terahertz spectroscopic characterization. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2020, 224, 117436. [Google Scholar] [CrossRef] [PubMed]
  29. Du, Y.; Zhang, H.; Xue, J.; Tang, W.; Fang, H.; Zhang, Q.; Li, Y.; Hong, Z. Vibrational spectroscopic study of polymorphism and polymorphic transformation of the anti-viral drug lamivudine. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2015, 137, 1158–1163. [Google Scholar] [CrossRef]
  30. Bawuah, P.; Zeitler, J.A. Advances in terahertz time-domain spectroscopy of pharmaceutical solids: A review. TrAC Trends Anal. Chem. 2021, 139, 116272. [Google Scholar] [CrossRef]
  31. Ruggiero, M.T. Invited Review: Modern Methods for Accurately Simulating the Terahertz Spectra of Solids. J. Infrared Millim. Terahertz Waves 2020, 41, 491–528. [Google Scholar] [CrossRef] [Green Version]
  32. Li, Q.; Bond, A.D.; Zeitler, J.A. Studying mechanical properties and phase transitions of aspirin polymorphs with terahertz spectroscopy and ab initio simulations. Int. Conf. Infrared Millim. Terahertz Waves IRMMW-THz 2019. [Google Scholar] [CrossRef]
  33. Aljohani, M.; Pallipurath, A.R.; McArdle, P.; Erxleben, A. A comprehensive cocrystal screening study of chlorothiazide. Cryst. Growth Des. 2017, 17, 5223–5232. [Google Scholar] [CrossRef]
  34. Chamorro-Posada, P. A Study of the Terahertz Spectra of Crystalline Materials (Polyethylene, Poly(Vinylidene Fluoride) form II, and α-D-Glucose) Using NDDO Semiempirical Methods. J. Appl. Spectrosc. 2018, 85, 552–559. [Google Scholar] [CrossRef]
  35. Shahab, S.; Almodarresiyeh, H.A.; Filippovich, L. Geometry optimization and excited state properties of the new symmetric (E)-Stilbene derivative for application in thermostable polarizing PVA-films: A combined experimental and DFT approach. J. Mol. Struct. 2016, 1119, 423–430. [Google Scholar] [CrossRef]
  36. Aree, T.; Jongrungruangchok, S. Enhancement of antioxidant activity of green tea epicatechins in β-cyclodextrin cavity: Single-crystal X-ray analysis, DFT calculation and DPPH assay. Carbohydr. Polym. 2016, 151, 1139–1151. [Google Scholar] [CrossRef] [PubMed]
  37. Dang, Y.; Li, S.; Liu, H.; Li, S.; Zhang, J.; Zhao, H. Low-frequency vibrations of indole derivatives by terahertz time-domain spectroscopy. J. Electron. Sci. Technol. 2016, 14, 329–336. [Google Scholar] [CrossRef]
  38. Hakey, P.M.; Hudson, M.R.; Allis, D.G. Examination of phencyclidine hydrochloride via cryogenic terahertz spectroscopy, solid-state density functional theory, and x-ray diffraction. J. Phys. Chem. A 2009, 113, 13013–13022. [Google Scholar] [CrossRef] [PubMed]
  39. Cai, Q.; Xue, J.; Wang, Q. Solid-state cocrystal formation between acyclovir and fumaric acid: Terahertz and Raman vibrational spectroscopic studies. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2017, 186, 29–36. [Google Scholar] [CrossRef]
  40. Strachan, C.J.; Rades, T.; Gordon, K.C. Raman spectroscopy for quantitative analysis of pharmaceutical solids. J. Pharm. Pharmacol. 2007, 59, 179–192. [Google Scholar] [CrossRef]
  41. Du, Y.; Cai, Q.; Xue, J. Structural investigation of the cocrystal formed between 5-fluorocytosine and fumaric acid based on vibrational spectroscopic technique. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2017, 178, 251–257. [Google Scholar] [CrossRef] [PubMed]
  42. Du, Y.; Zhang, H.; Xue, J.; Fang, H.; Zhang, Q.; Xia, Y.; Li, Y.; Hong, Z. Raman and terahertz spectroscopical investigation of cocrystal formation process of piracetam and 3-hydroxybenzoic acid. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2015, 139, 488–494. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Molecular structure of INH.
Figure 1. Molecular structure of INH.
Pharmaceutics 13 01303 g001
Figure 2. Molecular structure of PA.
Figure 2. Molecular structure of PA.
Pharmaceutics 13 01303 g002
Figure 3. Schematic diagram of molecular structure for possible theoretical cocrystal (a) form I, (b) form II, (c) form III.
Figure 3. Schematic diagram of molecular structure for possible theoretical cocrystal (a) form I, (b) form II, (c) form III.
Pharmaceutics 13 01303 g003aPharmaceutics 13 01303 g003b
Figure 4. The THz spectra of (a) INH, (b) PA, (c) physical mixture and (d) INH-PA monohydrate cocrystal in 0.2–1.6 THz spectral region.
Figure 4. The THz spectra of (a) INH, (b) PA, (c) physical mixture and (d) INH-PA monohydrate cocrystal in 0.2–1.6 THz spectral region.
Pharmaceutics 13 01303 g004
Figure 5. Comparisons of THz absorption spectra of (a) theoretical cocrystal form I, (b) theoretical cocrystal form II, (c) theoretical form cocrystal III and (d) experiment result.
Figure 5. Comparisons of THz absorption spectra of (a) theoretical cocrystal form I, (b) theoretical cocrystal form II, (c) theoretical form cocrystal III and (d) experiment result.
Pharmaceutics 13 01303 g005
Figure 6. Vibrational mode description of INH-PA monohydrate cocrystal at positions of (a) 0.88, (b) 1.20 and (c) 1.53 THz.
Figure 6. Vibrational mode description of INH-PA monohydrate cocrystal at positions of (a) 0.88, (b) 1.20 and (c) 1.53 THz.
Pharmaceutics 13 01303 g006
Figure 7. Experimental Raman spectra of (a) INH, (b) PA, (c) physical mixture and (d) INH-PA monohydrate cocrystal.
Figure 7. Experimental Raman spectra of (a) INH, (b) PA, (c) physical mixture and (d) INH-PA monohydrate cocrystal.
Pharmaceutics 13 01303 g007
Figure 8. Comparisons of Raman spectra of (a) theoretical cocrystal form I, (b) theoretical cocrystal form II, (c) theoretical form cocrystal III and (d) experiment result.
Figure 8. Comparisons of Raman spectra of (a) theoretical cocrystal form I, (b) theoretical cocrystal form II, (c) theoretical form cocrystal III and (d) experiment result.
Pharmaceutics 13 01303 g008
Figure 9. Typical bond length/Å of INH-PA monohydrate cocrystal.
Figure 9. Typical bond length/Å of INH-PA monohydrate cocrystal.
Pharmaceutics 13 01303 g009
Table 1. Vibrational mode assignment of the monohydrate cocrystal formed between INH and PA shown in THz spectra.
Table 1. Vibrational mode assignment of the monohydrate cocrystal formed between INH and PA shown in THz spectra.
Experimental Results f/THzTheoretical Results f/THzMode Assignment
0.880.86Out of plane rocking vibration of INH and PA
1.201.16Out of plane rocking vibration of INH and in-plane rocking vibration of PA
1.46-
1.531.53Out of plane rocking vibration of the ring of INH, torsional vibration of O11-C12-N13-H17 and H30-O29-C27-O28
Table 2. Vibrational mode assignment for characteristic peaks of the monohydrate cocrystal formed between INH and PA shown in Raman spectra.
Table 2. Vibrational mode assignment for characteristic peaks of the monohydrate cocrystal formed between INH and PA shown in Raman spectra.
ModeTheoretical Wavenumber/cm−1Experimental Wavenumber/cm−1Mode Assignment
ν1212225 (w)ω(H32-O31-C23-C22-O33-H34), Def(R2)
ν2238244 (w)τ(H32-O31-C23-C18-H24, O33-C22-C21-H26)
ν3316343 (w)δ(H34-O33-C22-C23-O31-H32), ω(H36-O35-H37)
ν4354387 (w)ω(H16-N14-H15), ρ(H36-O35-H37), δ(C2-C3-C4)
ν5416433 (vw)ω(H34-O33-C22), ρ(H36-O35-H37, H15-
N14-H16)
ν6456483 (vw)Def(R2), ρ(H36-O35-H37)
ν7513516 (vw)ω(H17-N13-C12), ρ(H36-O35-H37, H15-N14-H16)
ν8538537 (w)τ(H38-O29-C27, H25- C19-C18-H24)
ν9558557 (vw)ω(H36-O35-H37, H34-O33-C22)
ν10595595 (w)Def(R2), τ(36-O35-H37)
ν11682663 (m)Def(R1)
ν12728700 (vw)δ(H30-O29-C27-O28), Def(R2)
ν13798775 (s)δ(O33-C22-C21-H26), Def(R2), ω(H36-O35)
ν14857859 (vw)ω(H7-C1-C2-H8, H9-C4-C5-H10)
ν15903895 (m)δ(O11-C12-N13), ω(H15-N14-H16), ρ(H36-O35-H37)
ν16950950 (m)δ(C19-C20-C21, H24-C18-C19-H25)
ν1710121023 (m)Def(R1), ω(H15-N14-H16)
ν1810921093 (w)ρ(H26-C21-C20), Def(R2)
ν1911161115 (w)δ(H25-C19-C18-H24, H32-O31-C23-C18-H24)
ν2011831184 (m)ρ(H30-O29-C27), δ(H25-C19-C18-H24, H26-C21-C20-C27)
ν2112191213 (m)ρ(H8-C2-C3), ω(H15-N14-H16)
ν2212441245 (w)δ(H7-C1-C2-H8, H10-C4-C5-H9), ω(H15-N14-H16)
ν2312701292 (w)ρ(H26-C21-C22, H25-C19-C18-H24)
ν2413071319 (m)ρ(H24-C18-C19-H25), Def(R2)
ν2513181336 (s)ρ(H17-N13-C12), ω(H15-N14-H16), θ(N14-N13)
ν2613661390 (vw)δ(H30-O29-C27), ρ(H25-C19-C18-H24)
ν2713861409 (vw)δ(H34-O33-C22, H32-O31-C23)
ν2814401440 (vw)δ(H7-C1-C2-H8, H10-C4-C5-H9)
ν2915121524 (vw)ρ(H17-N13-C12, H7-C1-C2-H8, H10-C4-C5-H9)
ν30
ν31
ν32
1530
1632
1655
1542 (vw)
1610 (s)
1656 (m)
δ(H17-N13-C12), ρ(H7-C1-C2-H8, H10-C4-C5-H9)
Def(R1), δ(H36-O35-H37)
Def(R2), δ(H34-O33-C22, H36-O35-H37)
Vw—very weak; w—weak; m—medium; s—strong, θ—stretching; ρ—in-plane rocking vibration; δ—scissor; ω—out of plane rocking vibration; τ—torsion; Def—deformation.
Table 3. Changes of typical chemical bond lengths between INH, PA and the INH-PA monohydrate cocrystal.
Table 3. Changes of typical chemical bond lengths between INH, PA and the INH-PA monohydrate cocrystal.
Chemical BondBond Length/Å
INHPACocrystal
O11-C12
C12-C3
1.478
1.588
-
1.227
1.506
C12-N131.470-1.359
N13-N141.477-1.406
N14-H151.000-1.016
N14-H161.000-1.022
H34-O33 0.960.967
O33-C22-1.431.369
C22-C21
C21-H26
-1.394
1.100
1.384
1.081
C22-C23-1.3951.402
C23-O31-1.4301.370
O31-H32-0.9600.962
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bo, Y.; Fang, J.; Zhang, Z.; Xue, J.; Liu, J.; Hong, Z.; Du, Y. Terahertz and Raman Spectroscopic Investigation of Monohydrate Cocrystal of Antitubercular Isoniazid with Protocatechuic Acid. Pharmaceutics 2021, 13, 1303. https://0-doi-org.brum.beds.ac.uk/10.3390/pharmaceutics13081303

AMA Style

Bo Y, Fang J, Zhang Z, Xue J, Liu J, Hong Z, Du Y. Terahertz and Raman Spectroscopic Investigation of Monohydrate Cocrystal of Antitubercular Isoniazid with Protocatechuic Acid. Pharmaceutics. 2021; 13(8):1303. https://0-doi-org.brum.beds.ac.uk/10.3390/pharmaceutics13081303

Chicago/Turabian Style

Bo, Yanhua, Jiyuan Fang, Ziming Zhang, Jiadan Xue, Jianjun Liu, Zhi Hong, and Yong Du. 2021. "Terahertz and Raman Spectroscopic Investigation of Monohydrate Cocrystal of Antitubercular Isoniazid with Protocatechuic Acid" Pharmaceutics 13, no. 8: 1303. https://0-doi-org.brum.beds.ac.uk/10.3390/pharmaceutics13081303

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop