Terahertz identification and quantification of penicillamine enantiomers

Terahertz identification and quantification of penicillamine enantiomers 1 first-author 1 2 1 corresp Corresponding author. E-mail address: chenmin@sinap.ac.cn Corresponding author. E-mail address: chenmin@sinap.ac.cn chenmin@sinap.ac.cn 1 Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Jiading Campus, Shanghai 201800, China Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Jiading Campus, Shanghai 201800, China Center for Terahertz Research, Rensselaer Polytechnic Institute, 110 8th Street, Troy, NY 12180, USA Center for Terahertz Research, Rensselaer Polytechnic Institute, 110 8th Street, Troy, NY 12180, USA Identification and characterization of L -, D - and DL - penicillamine were demonstrated by Terahertz time-domain spectroscopy (THz-TDS). To understand the physical origins of the low frequency resonant modes, the density functional theory (DFT) was adopted for theoretical calculation. It was found that the collective THz frequency motions were decided by the intramolecular and intermolecular hydrogen bond interactions. Moreover, the quantification of penicillamine enantiomers mixture was demonstrated by a THz spectra fitting method with a relative error of less than 3.5%. This technique can be a valuable tool for the discrimination and quantification of chiral drugs in pharmaceutical industry. Terahertz spectroscopy Vibration modes Quantitative analysis 1 Introduction 1 Terahertz (THz) time-domain spectroscopy (TDS) is an extremely promising technique for chemical and biomedical applications since rotational and vibrational transitions of molecules and the low frequency vibration of crystal lattice lie within the far- and mid-infrared spectral range [ 1 1 ] . At present, a wide range of applications of THz-TDS technology have been reported [ 2 2 , 3 3 , 4 4 , 5 5 , 6 6 , 7 7 ] . Numerous structural biomolecules and chemicals have been investigated and the results indicate that the absorption characteristics in THz range are directly related to the composition and structure of molecules. The highly-sensitive THz spectroscopy provides us with a particular fingerprint method to discriminate molecules. At present, there is an increasing interest in studying low frequency vibration modes in chiral drugs because these modes can be used to identify isomers. The recent trend in pharmaceutical industry is to market the drugs in a pure enantiomeric form because different chiral forms of a drug may have different biological properties and clinical effect [ 8 8 ] . For some drugs, one kind of isomer is biologically-active and effective for therapy, while another has no effect or even toxicity. Hydrogen bonding and van der waals force play important roles in the interaction between drugs and organism systems. The corresponding vibrational modes of these weak interaction forces usually lie in THz range, which provides the theoretical basis for using THz technology to achieve molecular discrimination and study their interactions. Penicillamine is a thiol drug used in the treatment of rheumatoid arthritis. Only pure D -penicillamine is used in clinic since the L - form and DL -racemate are much more toxic, as shown by severe adverse reactions such as neuritis in patients treated with the DL -penicillamine [ 9 9 ] . A high degree of purification of D -penicillamine is essential for drug products. The X-ray crystal structure and IR spectra of a racemic mixture of D - and L - penicillamine have been studied. However, there is a lack of study on their vibrational modes [ 10 10 ] . In this paper, the low frequency absorption spectra of L -, D - and DL - Penicillamine are measured and characterized by THz-TDS. DFT theoretical calculations are carried out at the B3LYP/6-311G** level using Gaussian 03 packages and Gauss-View visualization program [ 11 11 ] . Excellent agreement with experiments is achieved in the values of calculated vibrational mode frequencies. In addition, quantification of penicillamine enantiomers in the mixtures is determined by fitting the obtained absorption THz spectra. 2 Experimental 1 2.1 Setup 2 The setup for THz-TDS system is the same as the one described in detail in Ref.[ 12 12 ]. The whole system is placed in a closed box purged with dry nitrogen gas in order to minimize the absorption of water vapor. The setup allows spectroscopic data to be recorded in the frequency range from 200 GHz to 2.2 THz with a dynamic range of about 1000 and a spectral resolution of better than 40 GHz. The measurements are considered accurate between 0.2 and 1.9 THz because of the excessive attenuation caused by the samples at higher frequency. 2.2 Data analysis 2 The method detailed in the paper by Duvillaret et al. [ 13 13 ] is used to extract the complex refractive index from the THz-TDS data. In short, the extraction of parameters is performed in two steps. Firstly, a reference spectrum is obtained in the absence of the sample, and then a spectrum in the presence of the sample is measured. The ratio of the sample data to the reference data is the complex transmission coefficient, which is a function of frequency. As noted in Ref.13, the error function is gradually minimized using a standard algorithm, which yields the desired values for complex refractive index. 2.3 Sample preparation 2 All samples were crystalline powder with purity higher than 98%, purchased from Sigma-Aldrich. Fig.1 Fig.1 shows the molecular structures of L -, D - and DL -penicillamine. All samples were carefully mixed with polyethylene (PE) powder with a weight ratio of 1:1 and a total weight of 200 mg, and then pressed into 1.554−1.642 mm thick pellets. PE was chosen because of low absorption coefficients ( < 5 cm -1 below 4 THz). 2.4 Computational methods 2 Theoretical calculation was performed by density functional theory (DFT). DFT method is widely used in molecular spectra calculations. The amount of its computation is considerably smaller than HF (Hartree-Fock) and MP2 (Moller-Plesset) methods under the same conditions, and its requirements are relatively low for computer system performance. Moreover, the DFT method can obtain more accurate calculation results relative to the HF method because of considering the electron correlation energy [ 14 14 ] . Geometry optimization and frequency analysis of the title compounds were performed using DFT B3LYP with 6-311++G** basis set [ 15 15 , 16 16 ] . The chiral penicillamine and their racemate exist with a zwitterion structure in the solid state. The two C-O bond lengths in carboxylate group are equal and three hydrogen atoms are attached to a nitrogen atom [ 10 10 ] . The X-ray structure of racemate DL -penicillamine provides useful information of the interaction between the enantiomers. The isolated and dimer models of L - and D -penicillamine were also investigated. 3 Results and discussion 1 The absorption spectra of L -, D - and DL -penicillamine measured at room temperature are shown in Fig.2 Fig.2 . The experimental results show significant difference in the absorption spectra of enantiomers ( L - and D -penicillamine) and their racemate ( DL -penicillamine). The absorption spectra of L - and D -penicillamine also show some differences in THz range. L - penicillamine has two absorption peaks located at 1.52 and 1.88 THz, while the peaks lie at 1.58 and 1.87 THz for D -penicillamine. This is in good agreement with the Raman spectrum of D -penicillamine measured by Howard-lock H E, et al [ 10 10 ] . In addition, two weak absorption peaks appear at 1.00 and 1.32 THz for D -penicillamine. The crystal structures of enantiomers ( L - and D -penicillamine) are orthorhombic with the space group P222 1 , and DL -penicillamine is monoclinic with the space group P2 1 /c [ 10 10 ] . It indicates that THz-TDS is sensitive to the change of the crystal structure and can be used for identification application. The vibrational frequencies were calculated by quantum chemical theory using the packages of Gaussion 03. The molecular structures of L -, D - and DL -penicillamine were made through geometrical optimization. No negative frequency mode was found in the results of calculation. Table 1 Table 1 shows the experimental modes and the calculated modes values at B3LYP/6-311++G** level for L -, D -, DL -penicillamine in 0.2−2.0 THz range. Fewer frequency modes were obtained from single molecule model compared with the dimer and racemic molecules models. It suggests that the molecular interactions are necessary to be considered in calculating low frequency vibration. 1 2 Sample 1 2 Exp./THz 1 2 Ref.[ 10 10 ] 2 1 B3LYP/6-311++G** 1 2 Approximate Mode Assignment 1 2 Vibrational 1 1 Single 1 1 Dimer (IR intensity) 1 1 1 1 1 1 1 1 1 1 0.50 (2.6857) 2 1 Butterfly (face to face) 1 1 1 1 1 1 1 1 1 1 1.05 (6.57) 1 1 Torsion 1 1 1 1 1 1 1 1 1 1 1 1 1.20 (3.2327) 1 1 Wiggle 1 1 1 1 L - 1 1 1 1 1 1 1 1 1.39 (0.1632) 1 1 Torsion 1 1 1 1 1 1 1.52 1 1 1 1 1.61 1 1 − 1 1 1 1 1 1 1 1 1.88 1 1 1 1 1 1 1.90 (1.3298) 1 1 Torsion 1 1 1 1 1 1 1 1 1 1 1 1 0.49 (2.7132) 2 1 Butterfly (face to face) 1 1 1 1 1.00 1 1 1 1 1 1 1.04 (6.4357) 1 1 Torsion 1 1 1 1 D - 1 1 1 1 1 1 1 1 1.19 (3.3303) 1 1 Wiggle 1 1 1 1 1.32 1 1 1 1 1 1 1.38 (0.1759) 1 1 Torsion 1 1 1 1 1 1 1.58 1 1 1.56 1 1 1.61 1 1 − 1 1 1 1 1 1 1 1 1.87 1 1 1.80 1 1 1 1 1.89 (1.2677) 1 1 Torsion 1 1 1 1 1 1 1 1 1 1 1 1 0.63 (1.42) 2 1 Butterfly (face to face) 1 1 1 1 1.36 1 1 1 1 1 1 1.29 (0.00) 1 1 Wiggle 1 1 1 1 DL - 1 1 1.53 1 1 1 1 1 1 1.41 (3.65) 1 1 Torsion 1 1 1 1 1 1 1.65 1 1 1 1 1.61 1 1 1.67 (2.20) 1 1 Torsion 1 1 According to the calculation by Gauss-View visualization program, approximate structural and low frequency vibrational modes were suggested. Two hydrogen atoms of the protonated amine take part in the formation of intra- and intermolecular hydrogen bonds. Two strong intermolecular hydrogen bonds N-H···O=C (0.152 nm) with N···O distance equal to 0.259 nm exist in a pair of racemate DL-penicillamine molecules. These interactions play an important role in maintaining the molecular structure in the solid state. The corresponding modes of racemate DL - penicillamine are shown in Fig.3 Fig.3 . The DFT calculation predicts four bands at 0.63, 1.29, 1.41 and 1.67 THz for DL -penicillamine. The vibrational mode at 0.63 THz is assigned as butterfly mode with face to face of D - and L -molecules of penicillamine. The mode at 1.29 THz comes from wiggle and modes at 1.41 and 1.67 THz are regarded as torsions which cause weak deformation vibrations of whole molecules. Therefore, the intra- and intermolecular hydrogen bonds are involved in the molecular motion. The low frequency modes obtained from the theoretical calculation indicate that the vibration modes are not localized the atomic motions but the whole molecular motions. When the dimer model is adopted to simulate the low frequency motion of L -, D - and DL - molecules, there are some differences of vibrational frequency between racemate and its enantiomers, which is caused by the different crystal structures of racemate (monoclinic P2 1 /c) and the enantiomer (orthorhombic P222 1 ). In addition, the vibrational frequencies of the enantiomers in L - and D - molecules are somewhat different. This can be understood from the following discussion. On one hand, the enantiomers have the same crystal structure parameters except the configuration, which would form different hydrogen bonds or relative motions with their respective adjacent molecules. It suggests that the frequency difference of penicillamine enantiomers obtained by experiment originates from their different configurations. On the other hand, theoretical model does not simulate the actual experimental condition perfectly and gives almost the same vibrational frequencies for L - and D -penicillamine. For example, multiple molecular interactions are not considered. This causes some discrepancy between the theoretical prediction and experimental results, such as in the value of peak frequency. To demonstrate the capability of quantification using THz-TDS, we measured the absorption spectra of solid powder mixtures of L - and D -penicillamine with different weight ratio. The weight content of L -penicillamine in the mixtures was 40%, 60% and 80%, respectively. Fig.4 Fig.4 shows the absorption spectra of L -, D -penicillamine and their mixture (60% L - and 40% D -). According to the Lambert-Beer law, if the absorption spectra of pure sample and their mixture are measured, the relative percentage of pure sample in the mixture can be obtained using the least square fit method. In this work, the absorption spectra of the penicillamine enantiomers and their mixture were fitted in the range from 1.2 to 1.7 THz, and the relative content of penicillamine enantiomers was obtained. A series of mixtures of penicillamine enantiomers with different weight ratio were measured and analyzed. The calculated weight percentage concentrations of L -penicillamine in the mixture are 40.9%, 62.1% and 81.5%, which are quite close to the actual concentrations of 40%, 60% and 80%. The analytic errors are 2.5%, 3.5% and 1.9%, respectively. The concentration of penicillamine enantiomers in their mixture can be quantitatively determined using the THz spectra fitting method. This provides a valuable tool for quantitative analysis of medicine and potential applications in pharmaceutical industry. 4 Conclusion 1 We have characterized the far infrared absorption properties of L -, D - and DL -penicillamine by THz-TDS in the region of 0.2-1.9 THz. The chiral molecules show distinct THz fingerprints due to their different crystal structure and molecular conformation. The DFT theory with the 6-311++G** basis set has been employed to calculate the low frequency resonance. More accurate information has been obtained from the dimer model than from the isolated molecule model, which suggests that weak molecular interaction needs to be taken into account in calculating low frequency vibration. The calculation results indicate that the characteristic vibrational frequencies come from various collective motions and both intra- and inter-molecular hydrogen bond interactions are involved. The quantification of penicillamine enantiomers mixture has been achieved by the THz spectra fitting method with a relative error less than 3.5%. The study shows that THz-TDS is a capable tool for qualitative and quantitative detection of chiral drugs. References Ferguson B , Zhang X C . Nat Mater , 2002 , 1 : 26 - 33 . Nat Mater Melinger J S , Laman N , Grischkowsky D . Appl Phys Lett , 2008 , 93 : 011102 . Appl Phys Lett Jameson A D , Kevek J W , Tomaino J L , et al . Appl Phys Lett , 2011 , 98 : 221111 . Appl Phys Lett Walker G C , Bowen J W , Naftaly M , et al . Meas Sci Technol , 2011 , 22 : 037003 . Meas Sci Technol Reid C B , Pickwell-Macpherson E , Laufer J G , et al . Phys Med Biol , 2010 , 55 : 4825 - 4838 . Phys Med Biol Kawase K , Shibuya T , Hayashi S , et al . C R Physique , 2010 , 11 : 510 - 518 . C R Physique Wang Weining , Wang Guo , Zhang Yan . Chin Phys B , 2011 , 20 : 123301 . Chin Phys B Kean W F , Lock C J L , Howard-Lock H E . Lancet , 1991 , 338 : 1565 - 1568 . Lancet Gotti R , Pomponio R , Andrisano V , et al . J Chromatogr A , 1999 , 844 : 361 - 369 . J Chromatogr A Howard-Lock H E , Lock C J L , Smalley P S . J Cryst Spectr Res , 1983 , 13 : 333 - 353 . J Cryst Spectr Res Frisch M J , Trucks G W , Schlegel H B , et al . GAUSSIAN03, Revision C. 02 , Gaussian Inc ., Pittsburgh, Pa. , 2003 . GAUSSIAN03, Revision C. 02 Han J G , Xu H , Zhu Z Y , et al . Chem Phys Lett , 2004 , 392 : 348 - 351 . Chem Phys Lett Duvillaret L , Garet F , Coutaz J L . IEEE J Sel Top. Quantum Electron , 1996 , 2 : 739 - 746 . IEEE J Sel Top. Quantum Electron Becke A D . J Chem Phys , 1993 , 98 : 5648 - 5652 . J Chem Phys Lee C , Yang W , Parr R G . Phys Rev B , 1988 , 37 : 785 - 789 . Phys Rev B Dewar M J S , Zoebisch E G , Healy E F , et al . J Am Chem Soc , 1985 , 107 : 3902 - 3909 . J Am Chem Soc 10.13538/j.1001-8042/nst.2013.01.005 Supported by National Natural Science Foundation of China (No.60907044) , the Major State Basic Research Development Program in China (No.2010CB834301) , Chinese Academy of Sciences Key Project of International Co-operation (No.GJHZ09058) and Science Foundation of the Chinese Academy of Sciences 2012-12-05 2013-02-01 2021-04-11 2013 Nuclear Science and Techniques 1 24