Original Research ARTICLE
Application of Solid-State NMR to Reveal Structural Differences in Cefazolin Sodium Pentahydrate From Different Manufacturing Processes
- 1National Institutes for Food and Drug Control, Beijing, China
- 2State Key Laboratory of Applied Organic Chemistry, College of Chemistry and Chemical Engineering, Lanzhou University, Lanzhou, China
- 3Zhejiang Institute for Food and Drug Control, Hangzhou, China
Solid-state Nuclear magnetic resonance, thermogravimetric analysis, X-ray diffraction, and Fourier-transform infrared spectroscopy were combined with theoretical calculation to investigate different crystal packings of α-cefazolin sodium obtained from three different vendors and conformational polymorphism was identified to exist in α-cefazolin sodium. Marginal differences observed among cefazolin sodium pentahydrate 1, 2, and 3 were speculated as being caused by the proportion of conformation 2.
Understanding the solid form of an active pharmaceutical ingredient (API) is a key part of its development. The differential and specific arrangements of molecules, such as co-crystals, solvates, polymorphs, and hydrates can determine a new chemical entity's crystal form. The mechanical and physical characteristics of pharmaceutical powders and the manufacturing process depend on their crystal morphology and structure, which are basically determined by the directionality and strength of interactions among molecules within the crystal lattice (Hancock et al., 2000). For organic crystals, hydrogen-bonding frameworks are often formed among molecules. Hydrogen bonding has received special focus because its directional nature causes anisotropy in molecular crystals (Payne et al., 1996). Therefore, the types and numbers of donors and acceptors of hydrogen bonds and the directionality of molecular interactions should have a large effect on crystal structures.
Cefazolin sodium (hereafter referred to as CEZ-Na) is a first-generation cephalosporin antibiotic (Figure 1). It has been widely used for treating various bacterial infections of the respiratory tract, skin, skin structure, bone and joints, and genital, urinary, and biliary tracts since the early 1970s (Nishida et al., 1969; Kariyone et al., 1970; Kusaba, 2009). It exists in four crystalline forms: α-crystal (pentahydrate, C14H13N8NaO4S3·5H2O), β-crystal (sesquihydrate, C14H13N8NaO4S3·1.5H2O), γ-crystal [C14H13N8NaO4S3·(HOCH2)2], and amorphous (Pikal and Delleman, 1989; Hu et al., 2002; Jacob et al., 2011). Interconversion occurs between these forms under certain conditions (Kamat et al., 1988; Osawa et al., 1988; Yang and Hu, 2005). The crystal structure of α-CEZ-Na was identified in 1983 (Stephenson and Diseroad, 2000) and manufactured by Fujisawa Co. Ltd., Japan; it was later transferred to Kyongbo Pharm. Co., Ltd., South Korea. Hu et al. identified another new form of α-CEZ-Na, which they obtained by using an isopropyl alcohol–water system (Hu et al., 2008); this was introduced to the Chinese market by Shenzhen Gosun Pharmaceutical Co., Ltd. These two α-CEZ-Na forms are similar: two symmetry-independent molecules are present in an asymmetric unit, where both molecules are well-ordered in the lattice, with large solvent tunnels. The thermogravimetric analysis (TGA) profiles are similar; however, their powder X-ray diffraction (PXRD) patterns (Figures S1, S2 in the electronic supplementary material) have certain differences. This indicates that the different manufacture processes affect the crystallization behavior of α-CEZ-Na, which results in different polymorphic forms.
The quality target product profile (QTPP) incorporates the intended quality attributes of the drug product to ensure efficacy and safety. The QTPP can be used as an effective tool to define the critical characteristics of the drug product and manufacturing process and thus ensure the intended performance of the product (Yu, 2008). Nuclear magnetic resonance (NMR) spectroscopy is a technique that exploits the magnetic attributes of nuclei and has been used to investigate the atomic environments of molecules. Solid-state NMR (SS-NMR) can provide detailed information on not just molecular structures but also intra- and intermolecular interactions, particularly molecular dynamics (Brown and Spiess, 2001). Thus, SS-NMR measurements are extremely useful for characterizing the crystal forms of pharmaceutical solids (Harris, 2007; Bai et al., 2010). It has been employed to characterize two amorphous forms and provide insight into the structures and molecular mobility of their amorphous forms (Frederick et al., 2013), detect interactions between the drug and substrate excipients (Remenar et al., 2004; Skotnicki et al., 2016), and gain information about the conformations of APIs and the interactions with water molecules (Chen et al., 2000). In this study, SS-NMR was combined with theoretical calculation to reveal the differences in the polymorphic structures of α-cefazolin sodium from different vendors.
Reagents and Materials
CEZ-Na pentahydrate 1 (C14H13N8NaO4S3·5H2O, batch no: L100100, L100200, L100300) was obtained from Fujisawa Co., Ltd., Kirihara-cho, Fujisawa-shi, Kanagawa, Japan. CEZ-Na pentahydrate 2 (C14H13N8NaO4S3·5H2O, batch no: 1368, 1391, 1401) was obtained from Kyongbo Pharm. Co., Ltd., Silok-ro, Ahsan-si, Choongchungnam-do, South Korea. CEZ-Na pentahydrate 3 (C14H13N8NaO4S3·5H2O, batch no: 1203283, 1203423, 1203403) was obtained from Shenzhen Gosun Pharmaceutical Co., Ltd., Guangdong, China. CEZ-Na sesquihydrate was prepared by crystal transformation of pentahydrate as described previously (Yang and Hu, 2005). Amorphous cefazolin sodium was prepared from the powder for injection produced by the freeze-drying technique. The purity of all samples was >99.1%. All samples were verified by solution 1H and 13C NMR and used in all measurements without further purification.
Powder X-Ray Diffraction
The PXRD patterns of the samples were recorded at room temperature (298 K) on a SmartLab diffractometer (Rigaku Corporation, Tokyo, Japan) with Cu Kα radiation (1.54 Å) at 45 kV and 200 mA. PXRD scans were recorded in continuous mode with a step size of 0.02° and scan rate of 1° per minute over a 2θ range of 2–60°.
TGA measurements were performed on a TAQ 500 (TA Instruments, Inc., New Castle, DE, USA). The α-CEZ-Na samples (3.40 mg) was weighed and analyzed at a heating rate of 20°C/min from 20 to 160°C and purged with nitrogen gas (25 mL/min).
Solid-State Nuclear Magnetic Resonance Spectroscopy
All SS-NMR spectra were obtained on a Bruker AVANCE II WB400 NMR spectrometer at 298 ± 2 K by using a 4 mm double-resonance (HX) probe at a 13C frequency of 100.61 MHz. The 13C chemical shifts were referenced externally to tetramethylsilane (TMS) of 0 ppm. The 13C CP/MAS experiments were carried out at a sample spinning rate of 10,000 ± 2 Hz with 4 mm BRUKER rotors and a contact time of 2 ms. The TPPM decoupling during data acquisition was provided by a proton decoupling field of 104.2 kHz. The recycle delays and number of scans of 13C CP/MAS experiments were 29 s and 320, respectively.
Two-dimensional 1H-13C heteronuclear correlation (HETCOR) experiments were performed with a Lee–Goldburg 1H-1H decoupling during the evolution period, which consisted of a 2.3 μs π/2 1H pulse and four Lee–Goldburg cycles per evolution increment. Different contact times were optimized and used to assign all chemical shifts of the complex 13C spectrum. The short contact-time HETCOR experiment (i.e., 50 μs) provided correlations between directly bonded nuclei, and the long contact-time experiment (i.e., 300 μs) provided possible long-range correlations. Eighty evolution points were collected for the indirect proton dimension, each with 34 scans. The 1H dimension was referenced internally to the methyl group of CEZ-Na at 2.73 ppm. All reported proton shifts were scaled by 0.578, which is the expected scaling factor for Lee–Goldburg homonuclear decoupling. The data were processed by multiplication by a cosine-square function in the t1 and t2 dimensions prior to the Fourier transformation.
Fourier Transform Infrared Spectroscopy
Fourier transform infrared (FTIR) spectra were recorded on an EQUINOX 55 spectrophotometer (Bruker, Zurich, Switzerland). Spectra were recorded at 26°C from 4,000 to 400 cm−1 with a resolution of 1 cm−1. Measurement samples were prepared by grinding α-CEZ-Na powder with KBr and compressing the mixture into thin pellets. Temperature-programmed IR was conducted as follows. (a) Method 1: IR spectra were recorded at 26, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, and 105°C at a heating rate of 5°C/min from 26 to 105°C. (b) Method 2: IR spectra were recorded at 50, 55, 56, 58, 60, 62, 64, 65, 66, 68, 70, 71, 72, 73, 74, 75, 76, 77, and 80°C at a heating rate of 0.5°C/min from 26 to 80°C.
Molecular Conformational Analysis
Conformational analysis was carried out by using the CONFORMER module of BIOVIA Material Studio with the COMPASS II force field. A torsion energy contour map was built up for the selected compounds, where the selected torsion angles were varied from 0° to 360° with a step size of 60° and constrained at each grid point whilst all other degrees of freedom were optimized for the molecule.
The cell structures were optimized by using the Forcite module of BIOVIA Materials Studio with the COMPASS II force field.
Results and Discussion
NMR Data for Cefazolin
The pentahydrate, sesquihydrate, and amorphous forms of CEZ-Na were distinctly distinguished by 13C SS-NMR spectral patterns with better resolved signals than those reported previously (Byrn et al., 1985). The carbon nuclei of the pentahydrate and sesquihydrate exhibited single peaks with no splitting (Figure 2) as compared with that of the amorphous form, which has broad resonances due to the distribution of conformers. Furthermore, comparison of the NMR data of pentahydrate and sesquihydrate showed that the chemical shifts and number of resonances for β-lactam (55–58 ppm) and carbonyl group (158–170 ppm) of the pentahydrate were distinctly different from those of the sesquihydrate.
Figure 2. 13C NMR spectra of CEZ-Na: (1) 13C NMR data in DMSO-d6, (2) 13C SS NMR spectra of α-CEZ-Na 1, (3) 13C SSNMR spectra of β-CEZ-Na, and (4) 13C SSNMR spectra of amorphous CEZ-Na.
The 13C NMR data obtained from CEZ-Na pentahydrate in inert DMSO-d6 was typically similar to the data obtained from SSNMR. Hence, the data assigned from 13C NMR served as guidance for the assignment of the 13C SS-NMR data. Distortionless enhancement of polarization transfer (DEPT), 1H–1H homonuclear correlation spectroscopy (COSY), heteronuclear single quantum coherence (HSQC), and heteronuclear multiple bond correlation (HMBC) were explored for accurate assignment of all the signals, especially the ambiguous signals at 163.5, 165.0, 166.4, 166.6, and 167.0 ppm. In particular, by the correlations between H6 and C8; H7 and C8, and C19; H10 and C12; H317 and C14; H20 and C19 and C25. Thus, all protons and carbons were assigned. The assignments of C12, 14, and C19 were corrected (Tori et al., 1981; Figures S3-1–S3-4 and Table 1).
The batch-to-batch consistency of the 13C SS-NMR data for the three forms of CEZ-Na pentahydrate 1, 2, and 3, respectively, indicated the stability of the manufacturing process and the reproducibility of the method (Figures S4–S6 in the electronic Supplementary Material). Distinct differences among CEZ-Na pentahydrate 1, 2, and 3 were observed with respect to the shape of peak C19 (i.e., intensity of the shoulder peak of C19) and the relative intensity of peak C14 (Figure 3). The 1H–13C HETCOR spectra were recorded at two contact times (50 and 300 μs) so that all chemical shifts observed in the complex 13C NMR spectrum of CEZ-Na could be assigned. The short contact-time sequence provided correlations between directly bonded nuclei, while the long-contact-time sequence provided possible long-range correlations, such as the dipolar coupling of the hydroxyl proton to the carbon to which the hydroxyl group is bonded (Figures S7–S8 and Figure 4). All carbon chemical shifts of CEZ-Na were assigned based on the significant correlations observed in the 1H–13C HETCOR spectra (300 μs) between H-4/C-6 and C-10; H-7/C-19; H-10/C-12; H3-17/C-14; H-20/C-19, and H-25/C-20 as well as the data obtained from solution NMR. Thus, all carbon chemical shifts of CEZ-Na were assigned (Table 2). Differences in signals among CEZ-Na pentahydrate 1, 2, and 3 were further ascertained at C14 and C19; this may be attributed to the minor configurational differences in the packing of molecules in a unit cell (Bai et al., 2010; Wang et al., 2012).
FTIR Data for Cefazolin
Different solid-state packing interactions provide small wavenumber shifts, which can be assessed if high-resolution vibrational measurements are applied (Arrais et al., 2003; Arrais and Savarino, 2009). The IR spectra of CEZ-Na pentahydrate 1, 2, and 3 were examined at 26°C with a high-resolution of 1 cm−1, and the results were undistinguished. Compared to the spectra reported in the literature (Pedroso and Salgado, 2013), no significant differences were found (Table 2 and Figure S9 in the electronic Supplementary Material). The IR spectral trends of pentahydrate 1, 2, and 3 were examined and analyzed in the range of 26–105°C (method 1). The general trends of the IR spectrum significantly changed at 75°C: [νOH(crystal water) and νNH(CONH)] were shifted by Δν of +17 and +11 cm−1, respectively, whereas νCO(β−lactam C = O) was shifted by Δν of −9 cm−1. In addition, the decreasing intensity of the νNH(CONH I band) was replaced by the increasing intensity of the absorption at 1,689 cm−1. This indicated changes in the intermolecular interactions of water molecules at approximately 75°C (Figure S10 in the electronic Supplementary Material). For further investigation of the differences between pentahydrate 1, 2, and 3 and the effect of temperature on the crystal form of CEZ-Na, a more detailed comparison was made (method 2). During the characteristic stage of the decreasing intensity of the νNH(CONHI band), pentahydrate 1 clearly changed at 76°C, and pentahydrate 2 and 3 accordingly changed at 71°C.
The different physical and chemical properties of a drug substance with polymorphic forms result in different values for the chemical reactivity, dissolution rate, melting point apparent solubility, vapor pressure, optical and mechanical properties, and density. These characteristics not only decide how or if a drug substance can be manufactured, but also affect the bioavailability, stability and dissolution of drug products. Therefore, the safety, efficacy, and quality of a drug product are largely depend on the polymorphism of its drug substance.
In conformational polymorphism, a molecule in an asymmetric unit can adopt different conformations through a controlled crystallization process. This is a critical part of polymorphism. Because SS-NMR is sensitive to variations in internuclear distances and local electronic structures resulting from various conformations, it is a powerful probe for characterization. In this study, the similarity of the TGA, FTIR, and SS 13C NMR spectra indicated a similar water content and crystal structure for CEZ-Na pentahydrate 1, 2, and 3. Two symmetry-independent molecules are present in an asymmetric unit, and both molecules are well-ordered in the lattice with large solvent tunnels that were identified by single-crystal XRD. XRD is a common experimental method for obtaining a detailed crystal structure in the ideal state. Thus, it failed to provide the integral solid state of the API and falsely indicated the presence of only a single structure for the polymorph. Although two independent molecules are present in an asymmetric unit of CEZ-Na pentahydrate, the symmetry of the two molecules should lead to single peaks without splitting in the 13C SS NMR spectra. All signals were displayed as single peaks, except the C-19 resonances; this indicates that molecules present in an asymmetric unit may adopt different conformations (Mimura et al., 2002; Liu et al., 2005).
Computational chemistry methods were utilized to study the conformation of CEZ-Na for exploring the differences with conformational polymorphism (Figure 5). The initial conformation 1 was reported by us earlier (Hu et al., 2008). The geometry of the conformation was energy-minimized by the CONFORMER module of the BIOVIA Material Studio with the COMPASS II force field to obtain conformation 2. The cell structures were optimized by using the Forcite module of BIOVIA Materials Studio with the COMPASS II force field. Similar crystal-lattice energies were generated (1: −817.8 kcal/mol; 2: −830.1 kcal/mol), which indicates that conformational polymorphism exists in CEZ-Na. The shoulder peak of C19 seen in SS-NMR data may be related to intramolecular hydrogen bond formation as seen in conformation 2. Although the intermolecular arrangements are quite different and should also affect C20 and C25, no significant change was found for the resolution of SSNMR.
Figure 5. Comparison of the molecular conformations of CEZ-Na (A) the conformations of α-CEZ-Na: 1. from crystal structure of α-CEZ-Na; 2. from calculation of molecular conformational analysis; (B) Global molecular arrangements: 1. from crystal structure of α-CEZ-Na; 2. from calculation of molecular conformational analysis).
Marginal differences were observed among CEZ-Na pentahydrate 1, 2, and 3 produced by different pharmaceutical companies, with PXRD, 13C SS-NMR spectra, temperature-programmed IR, and high-resolution TGA observations (Figure S11 in the electronic Supplementary Material) displaying subtle differences in the crystal structure of the pentahydrates. In particular, it was speculated that the shape of peak C19 (intensity of shoulder peak of C19) and the relative intensity of peak C14 were influenced by the proportion of conformation 2 differing among pharmaceutical companies (the proportion was greatest for CEZ-Na 2).
Understanding the solid-state structure of API at the local and long-range levels is crucial to developing strategies for the proper administration of these drugs. SS-NMR spectra and X-ray crystallography are the two techniques of choice to investigate the solid-state structure of synthetic pharmaceuticals and natural products. SS-NMR supplies structural information on powder amorphous solids for which single-crystal diffraction structures cannot be obtained. And SSNMR is uniquely sensitive to the local structure in polymorphic systems. Integrating all of these methods, such as TG, XRPD, and IR can provide insight into the local and long-range structure of polymorphic materials.
By comparing and analyzing of the SS-NMR data, it was speculated that conformation 2 was exhibited in three forms of CEZ-Na pentahydrate, and CEZ-Na pentahydrate 2 exhibited the highest percentage of conformation 2. These differences with respect to the percentage of conformation 2 in the molecular packing in the unit cell were attributed to a range of manufacturing processes, such as drying, wet granulation, spray-drying, and compaction. Future study will involve investigating the consistency in quality and the relationship between minor differences of polymorphs and the stability and efficacy.
The Food and Drug Agency (FDA) defines QTPP as the quality attributes related to the safety and efficacy of a product. After QTPP has been ascertained, the identification of the relevant critical quality attributes (CQAs) to ensure the intented product quality is the next important work. Understanding the solid form of an active pharmaceutical ingredient (API) is a key part of its development. So understanding the differences in the polymorphic structure of α-cefazolin sodium may help with making modifications to incorporate new knowledge with the product's development.
YT: design of experiments, analysis of experimental data; WW: SS-NMR data acquisition; W-BZ: PXRD and TGA data acquisition; J-QQ: theoretical study on molecular conformational analysis; C-QH: design of experiments, analysis of experimental data.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
WW thanks the Fundamental Research Funds for the Central Universities (lzujbky-2015-268) for the financial support.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2018.00113/full#supplementary-material
Arrais, A., Boccaleri, E., Croce, G., Milanesio, M., Orlando, R., and Diana, E. (2003). Synthesis, structural and spectroscopic study of the donor–acceptor complexes between fluorene and D2h cyano molecular building blocks. Cryst. Eng. Comm. 5, 388–394. doi: 10.1039/B308380B
Arrais, A., and Savarino, P. (2009). Raman spectroscopy is a convenient technique for the efficient evaluation of cyclodextrin inclusion molecular complexes of azo-dye colorants and largely polarisable guest molecules. J. Incl. Phenom. Macrocycl. Chem. 64, 73–81. doi: 10.1007/s10847-009-9538-5
Brown, S. P., and Spiess, H. W. (2001). Advanced solid-state NMR methods for the elucidation of structure and dynamics of molecular, macromolecular, and supramolecular systems. Chem. Rev. 101, 4125–4156. doi: 10.1021/cr990132e
Byrn, S. R., Gray, G., Pfeiffer, R. R., and Frye, J. (1985). Analysis of solid-state carbon-13 NMR spectra of polymorphs (benoxaprofen and nabilone) and pseudopolymorphls (cefazolin). J. Pharm. Sci. 74, 565–568. doi: 10.1002/jps.2600740516
Chen, L. R., Padden, B. E., Vippagunta, S. R., Munson, E. J., and Grant, D. J. (2000). Nuclear magnetic resonance and infrared spectroscopic analysis of nedocromil hydrates. Pharm. Res. 17, 619–624. doi: 10.1023/A:1007533419711
Hancock, B. C., Clas, S. D., and Christensen, K. (2000). Micro-scale measurement of the mechanical properties of compressed pharmaceutical powders. 1: the elasticity and fracture behaviour of microcrystalline cellulose. Int. J. Pharm. 209, 27–35. doi: 10.1016/S0378-5173(00)00541-X
Jacob, S., Nair, A. B., Patil, P. N., and Panda, B. P. (2011). Solid state crystallinity, amorphous state, and its implication in the pharmaceutical process. Int. J. Pharm. Sci. Res. 2, 472–482. doi: 10.13040/IJPSR.0975-8232
Kamat, M. S., Osawa, T., DeAngelis, R. J., Koyama, Y., and DeLuca, P. P. (1988). Estimation of the degree of crystallinity of cefazolin sodium by X-ray and infrared methods. Pharm. Res. 5, 426–429. doi: 10.1023/A:1015984300692
Mimura, H., Gato, K., Kitamura, S., Kitagawa, T., and Kohda, S. (2002). Effect of water content on the solid-state stability in two isomorphic clathrates of cephalosporin: cefazolin sodium pentahydrate (α Form) and FK041 hydrate. Chem. Pharm. Bull. 50, 766–770. doi: 10.1248/cpb.50.766
Nishida, M., Matsubara, T., Murakawa, T., Mine, Y., Yokota, Y., and Kuwahara, S. (1969). In vitro and in vivo evaluation of cefazolin, a new cephalosporin C derivative. Antimicrob. Agents Chemother. 9, 236–243.
Osawa, T., Kamat, M. S., and DeLuca, P. P. (1988). Hygroscopicity of cefazolin sodium: application to evaluate the crystallinity of freeze-dried products. Pharm. Res. 5, 421–425. doi: 10.1023/A:1015932316622
Payne, R. S., Roberts, R. J., Rowe, R. C., McPartlin, M., and Bashal, A. (1996). The mechanical properties of two forms of primidone predicted from their crystal structures. Int. J. Pharm. 145, 165–173. doi: 10.1016/S0378-5173(96)04760-6
Pikal, M. J., and Delleman, K. M. (1989). Stability testing of pharmaceuticals by high-sensitivity isothermal calorimetry at 25 deg C: cephalosporins in the solid and aqueous solution states. Int. J. Pharm. 50, 233–252. doi: 10.1016/0378-5173(89)90127-0
Remenar, J. F., Wenslow, R., Ostovic, D., and Peresypkin, A. (2004). Solid-state nuclear magnetic resonance determination of the physical form of BHA on common pharmaceutical excipients. Pharm. Res. 21, 185–188. doi: 10.1023/B:PHAM.0000012168.44628.ed
Skotnicki, M., Apperley, D. C., Aguilar, J. A., Milanowski, B., Pyda, M., and Hodgkinson, P. (2016). Characterization of two distinct amorphous forms of valsartan by solid-state NMR. Mol. Pharm. 13, 211–222. doi: 10.1021/acs.molpharmaceut.5b00646
Stephenson, G. A., and Diseroad, B. A. (2000). Structural relationship and desolvation behavior of cromolyn, cefazolin and fenoprofen sodium hydrates. Int. J. Pharm. 198, 167–177. doi: 10.1016/S0378-5173(00)00331-8
Keywords: cefazolin sodium, conformation, polymorphism, solid-state NMR, solvate
Citation: Tian Y, Wang WD, Zou W-B, Qian J-Q and Hu C-Q (2018) Application of Solid-State NMR to Reveal Structural Differences in Cefazolin Sodium Pentahydrate From Different Manufacturing Processes. Front. Chem. 6:113. doi: 10.3389/fchem.2018.00113
Received: 09 October 2017; Accepted: 26 March 2018;
Published: 10 April 2018.
Edited by:Simone Brogi, University of Siena, Italy
Reviewed by:Manuel Etzkorn, Heinrich Heine Universität Düsseldorf, Germany
Aldo Arrais, Università degli Studi del Piemonte Orientale, Italy
Copyright © 2018 Tian, Wang, Zou, Qian and Hu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chang-Qin Hu, firstname.lastname@example.org