Next Article in Journal
Angiogenesis in Regenerative Dentistry: Are We Far Enough for Therapy?
Next Article in Special Issue
Supramolecular Zn(II)-Dipicolylamine-Azobenzene-Aminocyclodextrin-ATP Complex: Design and ATP Recognition in Water
Previous Article in Journal
Assessment of the Efficacy and Clinical Utility of Different Circulating Tumor Cell (CTC) Detection Assays in Patients with Chemotherapy-Naïve Advanced or Metastatic Non-Small Cell Lung Cancer (NSCLC)
Previous Article in Special Issue
On the Single-Crystal Structure of Tenofovir Alafenamide Mono-Fumarate: A Metastable Phase Featuring a Mixture of Co-Crystal and Salt
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Enhancing the Physiochemical Properties of Puerarin via L-Proline Co-Crystallization: Synthesis, Characterization, and Dissolution Studies of Two Phases of Pharmaceutical Co-Crystals

1
Department of Chemistry, Zhejiang University, Hangzhou 310028, China
2
College of Chemistry, Chemical Engineering and Materials Science, Soochow University, Suzhou 215123, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2021, 22(2), 928; https://doi.org/10.3390/ijms22020928
Submission received: 2 December 2020 / Revised: 8 January 2021 / Accepted: 14 January 2021 / Published: 18 January 2021
(This article belongs to the Special Issue Macrocycles as Catalysts and Drug Carriers)

Abstract

:
Puerarin (PUE) is a Chinese traditional medicine known to enhance glucose uptake into the insulin cells to downregulate the blood glucose levels in the treatment of type II diabetes. Nevertheless, the bioavailability of pristine PUE is limited due to its poor solubility and low intestinal permeability. In this work, we demonstrate that the solubility of PUE can be significantly enhanced via its co-crystallization with L-Proline (PRO). Two crystalline phases, namely, the solvate-free form [PUE][PRO] (I) and the solvated form [PUE]2[PRO]∙EtOH∙(H2O)2 (II) are isolated. These two phases are characterized by single-crystal X-ray diffraction (SCXRD), powder X-ray diffraction (PXRD), Fourier-transformed infrared (FT-IR) spectra, nuclear magnetic resonance (NMR), and thermogravimetric analysis in association with differential scanning calorimetry (TGA-DSC). The solubility and dissolution rate of both I and II in water, gastrointestinal tract at pH 1.2, and phosphate buffer at pH 6.8 indicates a nearly doubled increase as compared to the pristine PUE. 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay of pristine PUE, I and II against murine colon cancer cell lines CT-26 and human kidney cell lines HEK-293 indicated that neither compound exhibits obvious cytotoxicity after 24 h. This work showcases that the readily available and biocompatible PRO can be a promising adjuvant to enhance the physicochemical properties of PUE toward orally administered drug formulation with improved pharmacokinetics.

Graphical Abstract

1. Introduction

One of the primary goals in pharmaceutical crystallization is to increase the solubility of the otherwise less soluble drugs to improve the bioavailability of the active pharmaceutical ingredients (APIs) toward the better pharmacokinetic outcome [1,2,3]. The core of such endeavor underlies the exploitation of different forms of APIs, such as co-crystals, salts, solvates (including hydrates), and polymorphs, followed by their physicochemical properties optimization [4]. Among these forms, co-crystallization is widely used to optimize the physicochemical properties of APIs due to the retention of its intrinsic pharmaceutical properties [5,6,7,8,9,10,11,12]. The pharmacokinetic properties can be regulated through shaping the co-crystals to improve its solubility and dissolution rate [13,14,15,16], mechanical stability [17], hygroscopicity [18], and ultimately bioavailability [10].
Puerarin (PUE, Scheme 1) is a major active ingredient in Chinese traditional medicine (CTM) that originates from the kudzu root known as Gegen [19]. PUE has been utilized for the treatment of diabetes mellitus in China since the 1990s to enhance the glucose uptake into the insulin cell, and subsequently downregulate the blood glucose levels in the type II diabetes treatment [20]. Notably, PUE has also been reported to inhibit the proliferation of a variety of cancer cells, such as colorectal cancer [21,22,23], colon cancer [23,24,25], and breast cancer [26]. However, because of its sizable conjugated skeleton, PUE exhibits poor solubility, leading to low oral bioavailability [27]. According to the biopharmaceutics classification system (BCS), PUE can be categorized as IV class due to its low solubility and limited intestinal permeability [28]. Besides, PUE also exhibits a short elimination half-life which further hindered its clinical adaption [29]. To improve the oral bioavailability and circulation time of PUE, a paucity of methodologies have been developed, such as the PUE encapsulation by HP-β-CD-PLGA nanoparticles [30], lipid nanoparticles [31], as well as RGD-appended PEGylated lipid nanoparticles [29].
In this work, we demonstrate that the solubility of PUE can be readily improved via its co-crystallization with L-proline (PRO, Scheme 1). PRO is well-known to exhibit good potential for co-crystal formation due to the presence of both hydrogen-bonding donor and acceptor sites in the molecule [32,33,34]. Two different phases of PUE-PRO co-crystals, namely the solvate-free form [PUE][PRO] (I) and the solvated form [PUE]2[PRO]∙EtOH∙(H2O)2 (II), have been prepared under different crystallization conditions and characterized spectroscopically. Compounds I and II represent the first cases of PUE-based co-crystals. As compared to those of the pristine PUE, the solubility and dissolution rates of I and II in water (pH 7.0), gastrointestinal tract at pH 1.2, and phosphate buffer at pH 6.8 indicated a nearly doubled increase. 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay of pristine PUE, I and II against murine colon cancer cell lines CT-26 and human kidney cell lines HEK-293 indicated that neither compound exhibits obvious cytotoxicity upon culturing the cells for an extended time of 24 h and 48 h, respectively. This work demonstrates that PUE co-crystal forms of I and II with the presence of readily available and biocompatible PRO can be a promising alternative to pristine PUE with better pharmacokinetic outcomes.

2. Results and Discussion

2.1. Synthesis and Co-Crystal Characterizations

The preferred route for the uptake of drugs, including Chinese traditional medicine (CTM), is the oral administration because of its convenience, low cost, and high persevering acquiescence compared to other routes. However, many of the bioactive ingredients of CTM are poorly soluble, leading to low oral bioavailability and delivery problems. Specifically for PUE, the oral bioavailability is only about 7% [32,33,34,35], and the dominant drug uptake at present is through intravenous injection. Thus, in the current study, the co-crystallization method using L-proline (PRO) as the partner excipient is employed to improve the water solubility of PUE. PRO was employed as a model excipient due to its good biocompatibility, its presence of both donor and acceptor sites for hydrogen-bonding, and its good solubility in water to guarantee the overall solubility of the product [36,37].
Both I and II were prepared as single crystals from equivalent molar of PUE and PRO in EtOH/H2O (9:1, v/v) under different precipitation temperatures. In the 1H nuclear magnetic resonance (NMR) spectra of I and II, peaks in the down-field area (δ 6.79–8.35 ppm) and the up-field area (δ 1.65–2.07 ppm) are diagnostic of the presence of PUE and PRO (Supplementary Figure S1). The triplet at δ 1.05 ppm in II further indicates the inclusion of EtOH solvate within the crystalline sample. The experimental powder X-ray diffraction (PXRD) patterns of I and II agree well with those simulated from single-crystal X-ray diffraction data (Figure 1), indicating highly phase purity of both complexes. In the Fourier-transform infrared (FT-IR) spectra (Supplementary Figure S2), peaks at 1605/1419 cm−1 (I) and 1590/1398 cm−1 (II) were assignable as the asymmetric/symmetric vibrations of COO in PRO [38,39,40]. Notably, the FT-IR spectra for both I and II feature no peak at around 1700 cm−1, indicating that the carboxylate remains deprotonated and the proton migrates to the –NH– to give a quaternary amine–carboxylate zwitterionic pair [41], as typically found in proline-based compounds [42,43].
The thermogravimetric analysis coupled with differential scanning calorimetry (TGA-DSC) reveals that I exhibit no obvious solvate loss until 224 °C (endothermic), followed by dramatic decomposition (Supplementary Figure S3a). This strongly suggests that I crystallizes as a solvate-free form. Meanwhile, the TGA-DSC curves of II show smooth and continuous solvate loss before 150 °C (Supplementary Figure S3b). This is followed by a plateau until 211 °C, at which point the PUE and PRO pair start to decompose. The total weight loss of 5.8% before 211 °C corresponds to the loss of two aqua and one EtOH molecules (calculated 7.8%), and the slight deviation might due to part of the solvate loss (mainly EtOH) during sample preparation.

2.2. Crystal Structure Analysis

Phase I crystallizes in the monoclinic space group P21 (Table 1) with its asymmetric unit featuring one PUE and one PRO molecules (Figure 2a). In the PRO moiety of I, the two C–O distances are nearly identical (1.249(3) Å and 1.255(3) Å), indicating that the carboxylate is in a conjugated and uncoordinated fashion, and the proton on the –COOH migrates to the –NH to give a zwitterion, corroborating that observed in the FT-IR spectrum. Such an averaged C–O distances are similar to those zwitterions based on metal-carboxylate pairs, such as trans-[Zn(iso-hmnH)2(H2O)2] and cis-[Zn(hmnH)2(H2O)2]·2H2O [44,45], and hemi-zwitterions such as [CuCl(iso-hmnH2)(iso-hmnH)]·H2O and [CuCl(hmnH2)(hmnH)] (iso-hmnH = 2-(hydroxymethyl) isonicotinate; hmnH = 2-(hydroxymethyl)nicotinate) [46].
Compound I is rich hydrogen-bonding interactions (Supplementary Table S1), and as a result, the PUE molecules are arranged to give a layered structure within the crystallographic ab plane (Figure 2b), and these layers are further stacked along the c direction in such a way that one-dimensional channels along the b direction is generated to accommodate PRO molecules.
Crystal structure analysis revealed that the solvated form of II belongs to the triclinic P1 space group (Table 1), and its asymmetric unit consists of two PUE, one PRO, one EtOH, and two aqua molecules (Figure 3a). In II, the PRO moiety also remains in its zwitterionic form similar to that identified in I. The two PUE molecules are in a head-to-tail arrangement as supported by hydrogen-bonding. With the additional presence of one PUE, and crystalline solvates EtOH and H2O, the hydrogen-bonding interactions in II are more complex than that of I (Supplementary Table S1). Similar to that found in I, PUE is also found to stack within the crystallographic ab plane (Figure 3b). Nevertheless, the presence of EtOH and aqua solvates marginally increased the separations between these layers to finally yield a sandwich-type structure in which the PRO and its hydrogen-bonded solvates reside in-between these layers. It is notable that EtOH is associated with PRO via a single hydrogen bond and can be readily extruded from the network without disturbing the overall structural integrity of II. This may serve to explain the slight deviation of experimental weight loss (5.8%) from that of the calculated (7.8%) observed in the TGA-DSC diagram (Supplementary Figure S3b).

2.3. Hirshfeld Surface Analysis of TAF, I and II

Hirshfeld surface (HS) analysis serves as a powerful tool to figure out the packing modes and intermolecular interactions in molecular crystals [47]. To evaluate and compare the nature of close interactions between the different lattice components that contributed to the supramolecular assembly of I and II, the three-dimensional (3D) HS and the two-dimensional (2D) fingerprint maps were generated using the Crystal Explorer program [48]. As illustrated in Figure 4, the deep-red spots on the surface portrait the shortest N–H∙∙∙N, N–H∙∙∙O, O–H∙∙∙N, and O–H∙∙∙O interactions, while other visible spots on the surfaces represent subtle H∙∙∙H contacts. Overall, the HS of I and II differ from each other in shape, reflecting the different weights of intermolecular contacts.
A summary of intermolecular contacts for PUE and PRO are illustrated with decomposed fingerprint plots (Figure 5) to showcase detailed atom-pair close contacts, including H∙∙∙H, O∙∙∙H/H∙∙∙O, and C∙∙∙H/H∙∙∙C contacts. Intermolecular interactions appear as discrete spikes in the 2D fingerprint plots. The quantitative analysis (Figure 6) shows that the H∙∙∙H and O∙∙∙H/H∙∙∙O interactions are the largest contributors to the total HS for I and II as compared to C∙∙∙H/H∙∙∙C and other contacts, represented by the large surfaces (blue) in the fingerprint plots.

2.4. Solubility and Dissolution Rates of PUE, I and II

It is well-established that the solid form of APIs has extensive influence on the solubility and dissolution rate of the drug. Consequently, it is vital to choose a suitable API solid form for successful drug expansion. According to the biopharmaceutics classification system (BCS) system, PUE is categorized as a class IV drug, meaning low solubility and low intestinal permeability. The supersaturation solid form of API may increase the bioavailability and permeability.
The comprehensive solubility of PUE, I, and II in the gastrointestinal tract at pH 1.2, phosphate buffer at pH 6.8, and H2O at pH 7.0 are depicted in Figure 7 as well as in Tables S2 and S3 (Supplementary Material). Both I and II exhibit nearly doubled solubility as compared to PUE with the solvated-free form I been comparable to the solvated form of II. In contrast, II also exhibits higher dissolution rate than both I and PUE. The exceptional solubility of I and II, and dissolution rate of II as compared to I and PUE might due to the presence of rich hydrogen-bonding interactions that enables simultaneously more competitive interactions with the solution media to facilitate the dissociation of the subcomponents of the co-crystal. Besides, the sandwich-type structure of II with solvate-accessible layer separations (Figure 3b) as compared to the channel structure of I may also result in the quick entrance of solution media into the crystal lattice to accelerate the dissolution of II. Thus, the increase of dissolution can be observed in the order II > I > PUE, and the solubility and dissolution of PUE were successfully improved via co-crystallization with PRO.

2.5. MTT Assays

PUE is widely reported to inhibit the proliferation of various cancer cells such as colon cancer [23,24,25]. It is reported that PUE is capable of inducing the apoptosis of colon cancer HT-29 cells [49]. To further evaluate the potential of PUE, I, and II as anti-cancer agents, we evaluated the biological activities of the pristine PUE, I, and II against the proliferation of CT-26 cancer cell lines (murine origin) and HEK-293 (human origin) cell lines using the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. As shown in Figure 8, PUE, I, and II exhibit negligible cytotoxicity toward CT-26 and HEK-293 cell lines upon culturing the cells with the drugs for 24 h. These complexes are also comparatively more toxic toward HEK-293 than CT-26 cell lines, and the CT-26 cells were also well survived upon incubating for additional 24 h (Supplementary Figure S4).

3. Materials and Methods

3.1. General

Puerarin (PUE) and L-Proline (PRO) were respectively purchased from J&K Scientific Ltd (Beijing, China) and Aladdin (Shanghai, China). Roswell Park Memorial Institute (RPMI) 1640 cell culture medium (for CT-26) and Dulbecco’s modified Eagle medium (DMEM) for HEK-293 were purchased from Bristol-Myers Squibb Trading Co. Ltd (New York, NY, USA). 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide (MTT) was obtained from Sigma (St. Louis, MO, USA). All the other chemicals were obtained directly from commercial sources and used as received. The 1H nuclear magnetic resonance (NMR) spectra were recorded on a Bruker DRX-400 spectrometer (Bruker, Ettlingen, Germany) and chemical shifts (δ) are referenced to the residual solvent peaks or internal standard TMS. FT-IR spectra were measured on a Thermo Scientific Nicolet iS50 FT-IR spectrometer (Thermo Fisher Scientific Co., Waltham, MA, USA) as KBr disks (400–4000 cm−1). Powder X-ray diffraction (PXRD) patterns were recorded on a Rigaku D/Max-2550PC micro-diffractometer (Rigaku Corporation, Tokyo, Japan). A rotating-anode Cu-target X-ray (λ= 1.5406 Å) was used, which was worked at 40 kV, 250 mA with the scanning ranges of 3.0 to 40.0° and scanning speed of 5° min−1 with an increasing step size of 0.02° and count time of 0.5–2 s. The TGA-DSC investigation was performed on a TA DSC Q100 differential filtering calorimeter (TA Instruments, New Castle, DE, USA) at the heating rate of 10 °C min−1 under a nitrogen stream of 50 cm3 min−1. Ultraviolet-visible (UV–vis) spectra were collected on a Thermo Scientific Evolution 300 UV spectrometer (Thermo Scientific, Waltham, MA, USA). The cytotoxicity data were evaluated on a Bio-Rad Elisa Plate Reader 680 (Bio-Rad, Hercules, CA, USA).

3.2. Synthesis of Co-Crystals I and II

The co-crystal was synthesized via the recrystallization method [50]. PUE (0.832 g, 2.0 mmol) and PRO (0.230 g, 2.0 mmol) were dissolved in 0.6 mL of EtOH (90%) and the solution was stirred at 90 °C for 30 min and then cooled to 50°C. The plate-shaped colorless crystals of I were formed after 24 h. The same amount of PUE and PRO were stirred at 50°C followed by cooling to room temperature yielded colorless fine needle crystals of II after two days.
Characterization Data for I.1H NMR (400 MHz, DMSO-d6) δ 9.60 (br, 1H), 8.35 (s, 1H), 7.93 (d, J = 8.8 Hz, 1H), 7.39 (d, J = 8.7 Hz, 2H), 7.01 (d, J = 8.8 Hz, 1H), 6.80 (d, J = 8.6 Hz, 2H), 5.01 (br, 2H), 4.80 (d, J = 9.8 Hz, 2H), 4.54 (br, 1H), 4.03 (t, J = 9.2 Hz, 1H), 3.63–3.73 (m, 2H), 3.18–3.43 (m, 6H), 2.97–3.04 (m, 1H), 1.98–2.05 (m, 1H), 1.85–1.98 (m, 1H), 1.75–1.82 (m, 1H), 1.65–1.74 (m, 1H). IR (KBr pellet): 3398 (br), 3248 (sh), 2923 (m), 1632 (s), 1592 (s), 1515 (s), 1446 (s), 1398 (s), 1268 (s), 1214 (s), 1176 (m), 1085 (s), 1061 (s), 1014 (s), 890 (m), 837 (m), 797 (m), 632 (m), 602 (m), 546 (m) cm−1.
Characterization Data for II.1H NMR (400 MHz, DMSO-d6) δ 9.56 (br, 3H), 8.35 (s, 2H), 7.93 (d, J = 8.8 Hz, 1H), 7.39 (d, J = 8.5 Hz, 2H), 7.06–6.91 (m, 1H), 6.80 (d, J = 8.6 Hz, 2H), 5.00 (br, 4H), 4.81 (br, 4H), 4.53 (br, 2H), 4.37 (br, 2H), 4.02 (br, 2H), 3.61–3.73 (m, 4H), 3.41–3.45 (m, 5H), 3.18–3.28 (m, 8H), 3.00–3.03 (m, 2H), 1.89–2.04 (m, 3H), 1.66–1.79 (m, 3H), 1.05 (t, J = 7.0 Hz, 3H). IR (KBr pellet): 3361 (br), 3266 (sh), 3084 (m), 3018 (m), 2924 (w), 2855 (w), 2736 (w), 2673 (w), 2607 (w), 1646 (s), 1605 (s), 1515 (s), 1445 (s), 1387 (s), 1358 (s), 1339 (s), 1258 (s), 1280 (s), 1180 (s), 1085 (s), 1043 (s), 1015 (m), 989 (s), 943 (s), 888 (w), 834 (m), 795 (w), 773 (w), 737 (w), 657 (w), 630 (m), 545 (m), 514 (m), 490 (w), 440 (w) cm−1.

3.3. Single-Crystal X-ray Crystallography

Single-crystal X-ray diffraction data were collected on a Bruker SMART APEX II diffractometer (Bruker, Ettlingen, Germany) using graphite-monochromated Cu Kα radiation (λ = 1.54178 Å). The collected data were corrected for absorption with SADABS [51]. The structures were solved by direct methods and refined by full-matrix least-squares on F2 using SHELXTL-2016 [52]. Non-hydrogen atoms were refined anisotropically. The H atoms on the phenolic and alcoholic –OH, water, and cationic NH2 were placed in calculated positions either by HFIX instructions or by Calc-OH program in the WinGX suite [53]. Crystallographic data for the crystal structures have been deposited in the Cambridge Crystallographic Data Centre (CCDC) with supplementary numbers of 2041676 (I) and 2041677 (II). These data can be obtained free of charge from the CCDC via www.ccdc.cam.ac.uk/data_request/cif. A summary of the key crystallographic data are listed in Table 1.

3.4. Solubility and Dissolution Measurement

The solubility and dissolution studies of pure PUE, I, and II were measured using a Thermo Scientific Evolution 300 UV spectrometer in deionized water at pH 7.0, gastrointestinal tract at pH 1.2, and phosphate buffer at pH 6.8 solution at 37 °C. The concentration was calculated utilizing a standard curve, which was made by measuring the absorbance of different concentrations at their λmax. A superfluity sum of the sample was added to 10 mL of H2O medium, pH 1.2 and pH 6.8 phosphate, respectively.
The oversaturated solution was stirred at 500 RPM by magnetic stirring at 37 °C for 24 h and the suspension was filtered through Whatman’s 0.45 mm syringe channel, diluted sufficiently and the absorbance was measured at their λmax. For the dissolution rate, the sample was specifically poured into 250 mL pH 7.0 water, pH 1.2 gastrointestinal tract, and pH 6.8 phosphates medium and then stirred at 150 RPM for 2h at 37 °C. Each solution was measured by the absorbance at λmax. At customary interims, the solution was replenished with 3 mL of the respective blank solution to maintain an overall steady volume.

3.5. Cytotoxicity Evaluation by MTT Assay

CT-26 and HEK-293 cell line was purchased from American Type Culture Collection (Manassas, VA, USA). The cell lines were cultured in RPMI 1640 (for CT-26) or DMEM (for HEK-293) containing 10% fetal bovine serum (FBS, Zhejiang Tianhang Biotechnology Co. Ltd., Hangzhou, China) and 1% penicillin/streptomycin (Zhejiang Tianhang Biotechnology Co. Ltd., Hangzhou, China). Cells grew as a monolayer and were detached upon confluence using trypsin (0.5% w/v in PBS, Zhejiang Tianhang Biotechnology Co. Ltd., Hangzhou, China). The cells were harvested from the cell culture medium by incubating in trypsin solution for 3 min. The cells were centrifuged, and the supernatant discarded. A 3 mL portion of serum-supplemented cell culture medium was added to neutralize any residual trypsin. The cells were re-suspended in serum-supplemented RPMI 1640 (for CT-26) or DMEM (for HEK-293) at a concentration of 5 × 104 cells per 1 mL. Cells were cultured at 37 °C and 5% CO2 for MTT studies.
The CT-26 or HEK-293 cells were seeded at a density of 1 × 104 cells per well in 100 µL of complete RPMI 1640 (CT-26) or DMEM (for HEK-293), and cultured for 16 h for attachment. The culture medium was then replaced by a serum-free medium containing various concentrations of pristine PUE, I, and II. After incubation with a period of 24 h and 48 h, the MTT solution (100 μL, 0.5 mg/mL in serum-free RPMI 1640 medium for CT-26 and DMEM for HEK-293) was added to replace the cell culture medium. After incubating the cells at 37 °C for 4 h, the MTT solution was removed and DMSO (100 μL) added to dissolve the formazan crystals formed, and the microplates were agitated for 5 min at a medium rate before the spectrophotometric measurement at a wavelength of 570 nm on a microplate reader. The untreated cells served as the 100% cell viability control, while the completely dead cells served as the blank. All experiments were carried out with four replicates. The relative cell viability (%) related to control cells was calculated by the formula below:
V % = [ A ] e x p e r i m e n t a l [ A ] b l a n k [ A ] c o n t r o l [ A ] b l a n k × 100 %
where V% is the percentage of cell viability, [A]experimental is the absorbance of the wells culturing the treated cells, [A]blank is the absorbance of the blank, and [A]control is the absorbance of the wells culturing untreated cells.

4. Conclusions

Two new forms of PUE have been successfully obtained and structurally authenticated via its co-crystallization with abundant and naturally-occurring PRO. The solubility and dissolution data suggested that the obtained co-crystal forms I and II exhibit marginal enhancement as compared to the pristine PUE. The solvated form of II further outperforms that of I probably due to its rich hydrogen-bonding interactions and sandwich-type structure that facilitate the quick diffusion of solution media. The use of the bio-compatible and non-toxic natural product as a co-delivery reagent may point to new directions in pharmaceutical co-crystallization to produce drug formulations with synergistic functions.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/22/2/928/s1.

Author Contributions

Conceptualization, X.H.; Methodology, M.I. and L.L.; Investigation, J.-W.W., L.L., K.-X.Y., C.-U.P. and J.S.; Data curation, X.H., K.-X.Y. and W.-H.Z.; Writing—original draft preparation, M.I. and L.L.; Writing—review and editing, G.T., X.H. and W.-H.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant No. 51873185.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are contained within this article and supplementary materials.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Suzuki, H.; Yakushiji, K.; Matsunaga, S.; Yamauchi, Y.; Seto, Y.; Sato, H.; Onoue, S. Amorphous solid dispersion of meloxicam enhanced oral absorption in rats with impaired gastric motility. J. Pharm. Sci. 2018, 107, 446–452. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Chen, Y.; Li, L.; Yao, J.; Ma, Y.-Y.; Chen, J.-M.; Lu, T.-B. Improving the solubility and bioavailability of apixaban via apixaban–oxalic acid cocrystal. Cryst. Growth Des. 2016, 16, 2923–2930. [Google Scholar] [CrossRef]
  3. Ulrich, J.; Frohberg, P. Problems, potentials and future of industrial crystallization. Front. Chem. Sci. Eng. 2013, 7, 1–8. [Google Scholar] [CrossRef] [Green Version]
  4. Duggirala, N.K.; Perry, M.L.; Almarsson, Ö.; Zaworotko, M.J. Pharmaceutical cocrystals: Along the path to improved medicines. Chem. Commun. 2016, 52, 640–655. [Google Scholar] [CrossRef]
  5. Surov, A.O.; Voronin, A.P.; Manin, A.N.; Manin, N.G.; Kuzmina, L.G.; Churakov, A.V.; Perlovich, G.L. Pharmaceutical cocrystals of diflunisal and diclofenac with theophylline. Mol. Pharm. 2014, 11, 3707–3715. [Google Scholar] [CrossRef] [PubMed]
  6. Steed, J.W. The role of co-crystals in pharmaceutical design. Trends Pharmacol. Sci. 2013, 34, 185–193. [Google Scholar] [CrossRef] [Green Version]
  7. Karimi-Jafari, M.; Padrela, L.; Walker, G.M.; Croker, D. Creating cocrystals: A Review of pharmaceutical cocrystal preparation routes and applications. Cryst. Growth Des. 2018, 18, 6370–6387. [Google Scholar] [CrossRef]
  8. Qiao, N.; Li, M.; Schlindwein, W.; Malek, N.; Davies, A.; Trappitt, G. Pharmaceutical cocrystals: An overview. Int. J. Pharm. 2011, 419, 1–11. [Google Scholar] [CrossRef]
  9. Shan, N.; Zaworotko, M.J. The role of cocrystals in pharmaceutical science. Drug Discov. Today 2008, 13, 440–446. [Google Scholar] [CrossRef]
  10. Korotkova, E.I.; Kratochvil, B. Pharmaceutical cocrystals. Procedia Chem. 2014, 10, 473–476. [Google Scholar] [CrossRef] [Green Version]
  11. Sun, C.C. Cocrystallization for successful drug delivery. Expert Opin. Drug Deliv. 2012, 10, 201–213. [Google Scholar] [CrossRef] [PubMed]
  12. Shan, N.; Perry, M.L.; Weyna, D.R.; Zaworotko, M.J. Impact of pharmaceutical cocrystals: The effects on drug pharma-cokinetics. Expert Opin. Drug Metab. Toxicol. 2014, 10, 1255–1271. [Google Scholar] [CrossRef] [PubMed]
  13. Childs, S.L.; Chyall, L.J.; Dunlap, J.T.; Smolenskaya, V.N.; Stahly, B.C.; Stahly, G.P. Crystal engineering approach to forming cocrystals of amine hydrochlorides with organic acids. Molecular complexes of fluoxetine hydrochloride with benzoic, succinic, and fumaric acids. J. Am. Chem. Soc. 2004, 126, 13335–13342. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Blagden, N.; De Matas, M.; Gavan, P.T.; York, P. Crystal engineering of active pharmaceutical ingredients to improve solubility and dissolution rates. Adv. Drug Deliv. Rev. 2007, 59, 617–630. [Google Scholar] [CrossRef] [PubMed]
  15. Thakuria, R.; Delori, A.; Jones, W.; Lipert, M.P.; Roy, L.; Rodríguez-Hornedo, N. Pharmaceutical cocrystals and poorly soluble drugs. Int. J. Pharm. 2013, 453, 101–125. [Google Scholar] [CrossRef]
  16. Elder, D.P.; Holm, R.; De Diego, H.L. Use of pharmaceutical salts and cocrystals to address the issue of poor solubility. Int. J. Pharm. 2013, 453, 88–100. [Google Scholar] [CrossRef]
  17. Aitipamula, S.; Wong, A.B.H.; Chow, P.S.; Tan, R.B. Novel solid forms of the anti-tuberculosis drug, isoniazid: Ternary and polymorphic cocrystals. CrystEngComm 2013, 15, 5877. [Google Scholar] [CrossRef]
  18. Tao, Q.; Chen, J.-M.; Ma, L.; Lu, T.-B. Phenazopyridine cocrystal and salts that exhibit enhanced solubility and stability. Cryst. Growth Des. 2012, 12, 3144–3152. [Google Scholar] [CrossRef]
  19. Zhou, Y.-X.; Zhang, H.; Peng, C. Puerarin: A review of pharmacological effects. Phytother. Res. 2014, 28, 961–975. [Google Scholar] [CrossRef]
  20. Kato, E.; Kawabata, J. Glucose uptake enhancing activity of puerarin and the role of C-glucoside suggested from activity of related compounds. Bioorg. Med. Chem. Lett. 2010, 20, 4333–4336. [Google Scholar] [CrossRef] [Green Version]
  21. Li, J.; Guo, C.; Lu, X.; Tan, W. Anti-colorectal cancer biotargets and biological mechanisms of puerarin: Study of molecular networks. Eur. J. Pharmacol. 2019, 858, 172483. [Google Scholar] [CrossRef] [PubMed]
  22. Li, R.; Ma, X.; Song, Y.; Zhang, Y.; Xiong, W.; Li, L.; Zhou, L. Anti-colorectal cancer targets of resveratrol and biological molecular mechanism: Analyses of network pharmacology, human and experimental data. J. Cell. Biochem. 2019, 120, 11265–11273. [Google Scholar] [CrossRef] [PubMed]
  23. Deng, X.-Q.; Zhang, H.-B.; Wang, G.-F.; Xu, D.; Zhang, W.-Y.; Wang, Q.-S.; Cui, Y.-L. Colon-specific microspheres loaded with puerarin reduce tumorigenesis and metastasis in colitis-associated colorectal cancer. Int. J. Pharm. 2019, 570, 118644. [Google Scholar] [CrossRef] [PubMed]
  24. Zhang, X.; Wang, B.; Mo, J. Puerarin 6″-O-xyloside possesses significant antitumor activities on colon cancer through inducing apoptosis. Oncol. Lett. 2018, 16, 5557–5564. [Google Scholar] [CrossRef] [Green Version]
  25. Zhang, W.-G.; Yin, X.-C.; Liu, X.-F.; Meng, K.-W.; Tang, K.U.N.; Huang, F.-L.; Xu, G.; Gao, J.I.E. Puerarin induces hepa-tocellular carcinoma cell apoptosis modulated by mapk signaling pathways in a dose-dependent manner. Anticancer Res. 2017, 37, 4425–4431. [Google Scholar]
  26. Liu, X.; Zhao, W.; Wang, W.; Lin, S.; Yang, L. Puerarin suppresses LPS-induced breast cancer cell migration, invasion and adhesion by blockage NF-κB and Erk pathway. Biomed. Pharmacother. 2017, 92, 429–436. [Google Scholar] [CrossRef]
  27. Yan, J.; Guan, Z.; Zhu, W.; Zhong, L.; Qiu, Z.-Q.; Yue, P.; Wu, W.-T.; Liu, J.; Huang, X. Preparation of puerarin chitosan oral nanoparticles by ionic gelation method and its related kinetics. Pharmaceutics 2020, 12, 216. [Google Scholar] [CrossRef] [Green Version]
  28. Li, H.; Dong, L.; Liu, Y.; Wang, G.; Wang, G.; Qiao, Y. Biopharmaceutics classification of puerarin and comparison of perfusion approaches in rats. Int. J. Pharm. 2014, 466, 133–138. [Google Scholar] [CrossRef]
  29. Dong, Z.; Guo, J.; Xing, X.; Zhang, X.; Du, Y.; Lu, Q. RGD modified and PEGylated lipid nanoparticles loaded with puerarin: Formulation, characterization and protective effects on acute myocardial ischemia model. Biomed. Pharmacother. 2017, 89, 297–304. [Google Scholar] [CrossRef]
  30. Tao, H.-Q.; Meng, Q.; Li, M.-H.; Yu, H.; Liu, M.-F.; Du, D.; Sun, S.-L.; Yang, H.-C.; Wang, Y.-M.; Ye, W.; et al. HP-β-CD-PLGA nanoparticles improve the penetration and bioavailability of puerarin and enhance the therapeutic effects on brain ischemia–reperfusion injury in rats. Naunyn-Schmiedeberg’s Arch. Pharmacol. 2012, 386, 61–70. [Google Scholar] [CrossRef]
  31. Luo, C.-F.; Yuan, M.; Chen, M.-S.; Liu, S.-M.; Zhu, L.; Huang, B.-Y.; Liu, X.-W.; Xiong, W. Pharmacokinetics, tissue dis-tribution and relative bioavailability of puerarin solid lipid nanoparticles following oral administration. Int. J. Pharm. 2011, 410, 138–144. [Google Scholar] [CrossRef] [PubMed]
  32. Song, Y.; Wang, L.-Y.; Liu, F.; Li, Y.-T.; Wu, Z.-Y.; Yan, C.-W. Simultaneously enhancing the In Vitro/In Vivo performances of acetazolamide using proline as a zwitterionic coformer for cocrystallization. CrystEngComm 2019, 21, 3064–3073. [Google Scholar] [CrossRef]
  33. Liu, M.; Hong, C.; Yao, Y.; Shen, H.; Ji, G.; Li, G.; Xie, Y. Development of a pharmaceutical cocrystal with solution crys-tallization technology: Preparation, characterization, and evaluation of myricetin-proline cocrystals. Eur. J. Pharm. Bio-pharm. 2016, 107, 151–159. [Google Scholar] [CrossRef] [PubMed]
  34. Jensen, K.T.; Löbmann, K.; Rades, T.; Grohganz, H. Improving Co-Amorphous Drug Formulations by the Addition of the Highly Water Soluble Amino Acid, Proline. Pharmaceutics 2014, 6, 416–435. [Google Scholar] [CrossRef] [PubMed]
  35. Yi, T.; Tang, D.; Wang, F.; Zhang, J.; Zhang, J.; Wang, J.; Xu, X.; Zhang, J. Enhancing both oral bioavailability and brain penetration of puerarin using borneol in combination with preparation technologies. Drug Deliv. 2017, 24, 422–429. [Google Scholar] [CrossRef] [Green Version]
  36. Dou, M.; Lu, C.; Sun, Z.; Rao, W. Natural cryoprotectants combinations of l-proline and trehalose for red blood cells cryopreservation. Cryobiology 2019, 91, 23–29. [Google Scholar] [CrossRef]
  37. Thangavel, P.; Ramachandran, B.; Kannan, R.; Muthuvijayan, V. Biomimetic hydrogel loaded with silk and l -proline for tissue engineering and wound healing applications. J. Biomed. Mater. Res. Part B Appl. Biomater. 2017, 105, 1401–1408. [Google Scholar] [CrossRef]
  38. Papageorgiou, S.K.; Kouvelos, E.P.; Favvas, E.P.; Sapalidis, A.A.; Romanos, G.E.; Katsaros, F.K. Metal–carboxylate inter-actions in metal–alginate complexes studied with FTIR spectroscopy. Carbohydr. Res. 2010, 345, 469–473. [Google Scholar] [CrossRef]
  39. Mary, Y.S.; Ushakumari, L.; Harikumar, B.; Varghese, H.T.; Panicker, C.Y. FT-IR, FT-raman and SERS spectra of L-proline. J. Iran. Chem. Soc. 2009, 6, 138–144. [Google Scholar] [CrossRef]
  40. Yuan, F.-L.; Yuan, Y.-Q.; Chao, M.-Y.; Young, D.J.; Zhang, W.-H.; Lang, J.-P. Deciphering the structural relationships of five cd-based metal–organic frameworks. Inorg. Chem. 2017, 56, 6522–6531. [Google Scholar] [CrossRef]
  41. Ratajczak, H.; Barycki, J.; Pietraszko, A.; Baran, J.; Debrus, S.; May, M.; Venturini, J. Preparation and structural study of a novel nonlinear molecular material: The l-histidinum dihydrogenarsenate orthoarsenic acid crystal. J. Mol. Struct. 2000, 526, 269–278. [Google Scholar] [CrossRef]
  42. Abu-Nawwas, A.-A.H.; Cano, J.; Christian, P.; Mallah, T.; Rajaraman, G.; Teat, S.J.; Winpenny, R.E.P.; Yukawa, Y. An Fe (III) wheel with a zwitterionic ligand: The structure and magnetic properties of [Fe(OMe)2(proline)]12[ClO4]12. Chem. Commun. 2004, 314–315. [Google Scholar] [CrossRef] [PubMed]
  43. An, H.; Lan, Y.; Li, Y.; Wang, E.; Hao, N.; Xiao, D.; Duan, L.; Xu, L. A novel chain-like polymer constructed from heter-opolyanions covalently linked by lanthanide cations: (C5H9NO2)2[La(H2O)7CrMo6H6O24]·11H2O (Proline=C5H9NO2). Inorg. Chem. Commun. 2004, 7, 356–358. [Google Scholar] [CrossRef]
  44. Armaghan, M.; Shang, X.J.; Yuan, Y.Q.; Young, D.J.; Zhang, W.H.; Hor, T.S.A.; Lang, J.P. Metal−organic frameworks via emissive metal-carboxylate zwitterion intermediates. ChemPlusChem 2015, 80, 1231–1234. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Liu, Y.; Lin, S.-X.; Niu, R.-J.; Liu, Q.; Zhang, W.-H.; Young, D.J. Zinc and cadmium complexes of pyridinemethanol car-boxylates: Metal carboxylate zwitterions and metal–organic frameworks. ChemPlusChem 2020, 85, 832–837. [Google Scholar] [CrossRef] [PubMed]
  46. Armaghan, M.; Young, D.J.; Wu, D.; Wei, Y.; Yuan, F.-L.; Ng, S.W.; Amini, M.M.; Zhang, W.-H.; Young, D.J.; Chen, H.W.; et al. Isolation of first row transition metal-carboxylate zwitterions. RSC Adv. 2015, 5, 42978–42989. [Google Scholar] [CrossRef] [Green Version]
  47. Spackman, M.A.; Jayatilaka, D. Hirshfeld surface analysis. CrystEngComm 2009, 11, 19–32. [Google Scholar] [CrossRef]
  48. Wolff, S.; Grimwood, D.; McKinnon, J.; Turner, M.; Jayatilaka, D.; Spackman, M. Crystal Explorer; The University of Western Australia: Crawley, Australia, 2012. [Google Scholar]
  49. Yu, Z.; Li, W. Induction of apoptosis by puerarin in colon cancer HT-29 cells. Cancer Lett. 2006, 238, 53–60. [Google Scholar] [CrossRef]
  50. Sugandha, K.; Kaity, S.; Mukherjee, S.; Isaac, J.; Ghosh, A. Solubility enhancement of ezetimibe by a cocrystal engineering technique. Cryst. Growth Des. 2014, 14, 4475–4486. [Google Scholar] [CrossRef]
  51. Bruker. APEX2, SAINT and SADABS; Bruker ACS Inc.: Madison, WI, USA, 2014. [Google Scholar]
  52. Sheldrick, G.M. Crystal structure refinement with SHELXL. Acta Cryst. 2015, 71, 3–8. [Google Scholar] [CrossRef]
  53. Farrugia, L.J. WinGX suite for smallmolecule single-crystal crystallography. J. Appl. Cryst. 1999, 32, 837–838. [Google Scholar] [CrossRef]
Scheme 1. The structure of Puerarin and l-Proline.
Scheme 1. The structure of Puerarin and l-Proline.
Ijms 22 00928 sch001
Figure 1. PXRD patterns of I and II showing a good agreement between the experimental and those simulated from single-crystal diffraction data and thus high phase purity of I and II.
Figure 1. PXRD patterns of I and II showing a good agreement between the experimental and those simulated from single-crystal diffraction data and thus high phase purity of I and II.
Ijms 22 00928 g001
Figure 2. The asymmetric unit (a) and hydrogen-bonding network (b) of I. In b, the structure is propagated along the crystallographic a and c directions showing that PRO (distinguished by dark magenta for clarity purpose) resides inside the 1D channels generated by stacking of the PUE molecules. Color codes: O (red), N (blue), C (black), H (light yellow).
Figure 2. The asymmetric unit (a) and hydrogen-bonding network (b) of I. In b, the structure is propagated along the crystallographic a and c directions showing that PRO (distinguished by dark magenta for clarity purpose) resides inside the 1D channels generated by stacking of the PUE molecules. Color codes: O (red), N (blue), C (black), H (light yellow).
Ijms 22 00928 g002
Figure 3. The asymmetric unit (a) and hydrogen-bonding network (b) of II. In b, the structure is propagated along the crystallographic b and c directions showing that PRO, one EtOH, and a pair of H2O molecules are sandwiched in-between the layers generated by PUE. The PRO and crystalline solvates were distinguished by dark magenta for clarity purposes. Color codes: O (red), N (blue), C (black), H (light yellow).
Figure 3. The asymmetric unit (a) and hydrogen-bonding network (b) of II. In b, the structure is propagated along the crystallographic b and c directions showing that PRO, one EtOH, and a pair of H2O molecules are sandwiched in-between the layers generated by PUE. The PRO and crystalline solvates were distinguished by dark magenta for clarity purposes. Color codes: O (red), N (blue), C (black), H (light yellow).
Ijms 22 00928 g003
Figure 4. Hirshfeld surfaces for the PUE with its neighboring molecules in the crystals of I (a), II-A (b), and II-B (c).
Figure 4. Hirshfeld surfaces for the PUE with its neighboring molecules in the crystals of I (a), II-A (b), and II-B (c).
Ijms 22 00928 g004
Figure 5. Two-dimensional fingerprint plots of I (a), II-A (b), and II-B (c), showing the contribution of each type of contact within the molecule.
Figure 5. Two-dimensional fingerprint plots of I (a), II-A (b), and II-B (c), showing the contribution of each type of contact within the molecule.
Ijms 22 00928 g005
Figure 6. Relative contribution ratios of each type of contacts within the structure of I, and two parts of II (denoted as II-A and II-B), showing that the presence of water molecules in II raises the collective contribution of O∙∙∙H/H∙∙∙O amongst these interactions.
Figure 6. Relative contribution ratios of each type of contacts within the structure of I, and two parts of II (denoted as II-A and II-B), showing that the presence of water molecules in II raises the collective contribution of O∙∙∙H/H∙∙∙O amongst these interactions.
Ijms 22 00928 g006
Figure 7. A summary of the solubility profile of PUE, I and II in the gastrointestinal tract at pH 1.2, phosphate buffer at pH 6.8, and H2O at pH 7.0 (a), as well as the time-dependent dissolution profile of these compounds under respective conditions: gastrointestinal tract at pH 1.2 (b), phosphate buffer at pH 6.8 (c), and H2O at pH 7.0 (d), reflecting the dissolution rate of each compound. Data are expressed as average values from three replicates.
Figure 7. A summary of the solubility profile of PUE, I and II in the gastrointestinal tract at pH 1.2, phosphate buffer at pH 6.8, and H2O at pH 7.0 (a), as well as the time-dependent dissolution profile of these compounds under respective conditions: gastrointestinal tract at pH 1.2 (b), phosphate buffer at pH 6.8 (c), and H2O at pH 7.0 (d), reflecting the dissolution rate of each compound. Data are expressed as average values from three replicates.
Ijms 22 00928 g007
Figure 8. The cell viability data of PUE, I, and II against murine colon cancer CT-26 (a) and human kidney cell lines HEK-293 (b) upon culturing 24 h.
Figure 8. The cell viability data of PUE, I, and II against murine colon cancer CT-26 (a) and human kidney cell lines HEK-293 (b) upon culturing 24 h.
Ijms 22 00928 g008
Table 1. Summary of crystallographic data for I and II.
Table 1. Summary of crystallographic data for I and II.
CompoundIII
Empirical formulaC21H20O9∙C5H9NO22(C21H20O9)·C5H9NO2·2H2O·C2H5OH
Formula Weight531.501029.97
Crystal systemmonoclinicTriclinic
Space groupP21P1
a10.2463 (4)6.3123 (2)
b7.8656 (3)12.1389 (4)
c15.4520 (6)15.7274 (6)
α9081.7330 (10)
β97.7830 (10)89.7330 (10)
γ9087.4870 (10)
V31233.86 (8)1191.43 (7)
Z21
Dc/(g cm−3)1.4311.436
μ (Mo-Kα)/mm−10.9520.974
F(000)560544
Total reflections13,35740,998
Unique reflections44149385
No observations44119322
No parameters349673
Flack parameter0.03 (3)−0.01 (3)
Rint0.02540.0293
R a0.02930.0332
wR b0.07740.0957
GOF c1.0571.080
a R1 = Σ||Fo| − |Fc||/Σ|Fo|, b wR2 = {Σ[ω(Fo2 − Fc2)2]/Σ[ω(Fo2)2]}1/2, and c GOF = {Σ[ω(Fo2 − Fc2)2]/(n − p)}1/2, where n is the number of reflections, and p the total number of parameters refined.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Inam, M.; Liu, L.; Wang, J.-W.; Yu, K.-X.; Phan, C.-U.; Shen, J.; Zhang, W.-H.; Tang, G.; Hu, X. Enhancing the Physiochemical Properties of Puerarin via L-Proline Co-Crystallization: Synthesis, Characterization, and Dissolution Studies of Two Phases of Pharmaceutical Co-Crystals. Int. J. Mol. Sci. 2021, 22, 928. https://doi.org/10.3390/ijms22020928

AMA Style

Inam M, Liu L, Wang J-W, Yu K-X, Phan C-U, Shen J, Zhang W-H, Tang G, Hu X. Enhancing the Physiochemical Properties of Puerarin via L-Proline Co-Crystallization: Synthesis, Characterization, and Dissolution Studies of Two Phases of Pharmaceutical Co-Crystals. International Journal of Molecular Sciences. 2021; 22(2):928. https://doi.org/10.3390/ijms22020928

Chicago/Turabian Style

Inam, Muhammad, Lu Liu, Jian-Wei Wang, Ka-Xi Yu, Chi-Uyen Phan, Jie Shen, Wen-Hua Zhang, Guping Tang, and Xiurong Hu. 2021. "Enhancing the Physiochemical Properties of Puerarin via L-Proline Co-Crystallization: Synthesis, Characterization, and Dissolution Studies of Two Phases of Pharmaceutical Co-Crystals" International Journal of Molecular Sciences 22, no. 2: 928. https://doi.org/10.3390/ijms22020928

APA Style

Inam, M., Liu, L., Wang, J. -W., Yu, K. -X., Phan, C. -U., Shen, J., Zhang, W. -H., Tang, G., & Hu, X. (2021). Enhancing the Physiochemical Properties of Puerarin via L-Proline Co-Crystallization: Synthesis, Characterization, and Dissolution Studies of Two Phases of Pharmaceutical Co-Crystals. International Journal of Molecular Sciences, 22(2), 928. https://doi.org/10.3390/ijms22020928

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