Literature DB >> 25705449

Crystal structure of zwitterionic 3-(2-hy-droxy-2-phospho-nato-2-phosphono-eth-yl)imidazo[1,2-a]pyridin-1-ium monohydrate (minodronic acid monohydrate): a redetermination.

Annalisa Airoldi1, Piergiorgio Bettoni1, Monica Donnola1, Gianluca Calestani2, Corrado Rizzoli2.   

Abstract

In a previous study, the X-ray structure of the title compound, C9H12N2O7P2·H2O, was reported [Takeuchi et al., (1998 ▶). Chem. Pharm. Bull. 46, 1703-1709], but neither atomic coordinates nor details of the geometry were published. The structure has been redetermined with high precision as its detailed knowledge is essential to elucidate the presumed polymorphism of minodronic acid monohydrate at room temperature. The mol-ecule crystallizes in a zwitterionic form with cationic imidazolium[1,2a]pyridine and anionic phospho-nate groups. The dihedral angle formed by the planes of the pyridine and imidazole rings is 3.55 (9)°. A short intra-molecular C-H⋯O contact is present. In the crystal, mol-ecules are linked by O-H⋯O, N-H⋯O and C-H⋯O hydrogen bonds and π-π inter-actions [centroid-to-centroid distance = 3.5822 (11) Å], forming a three-dimensional structure.

Entities:  

Keywords:  bis­phospho­nate; crystal structure; hydrogen bonds; minodronic acid; polymorphism; redetermination

Year:  2015        PMID: 25705449      PMCID: PMC4331878          DOI: 10.1107/S2056989014026863

Source DB:  PubMed          Journal:  Acta Crystallogr E Crystallogr Commun


Chemical context

Minodronic acid, (1-hy­droxy-2-(imidazo[1,2-a]pyridin-3-yl)ethane-1,1-bis­(phospho­nic acid), has excellent bone resorption inhibitory activity, as well as anti-inflammatory, analgesic and anti­pyretic activities and is useful for the treatment of diseases in which an increased bone resorption participates (Tanishima & Morio, 2013 ▸; Yamane et al., 2003 ▸; Ito et al., 2010 ▸; Sato et al., 2006 ▸). For practical uses and solid pharmaceutical preparations, the monohydrate form is preferred, giving more stable crystals with respect to the anhydrous and polyhydrate forms. By developing a pharmaceutical preparation (Takeuchi et al., 1993 ▸), it was found that monohydrate crystals include two crystal forms, named D and E, having the same powder X-ray diffraction pattern but different dehydration temperatures. In fact, according to TG–DSC thermogravimetric analyses, crystal D has a dehydration peak temperature of 408 to 422 K, whereas for E it is 433 to 443 K. The monohydrate crystals of D and E are both produced by recrystallization of the free acid from aqueous hydro­chloric acid solutions by gradual cooling of a heat-dissolved solution under mild stirring conditions, followed by drying the crystals at 313–333 K under reduced pressure. The stirring mode and cooling conditions are the key factors in determining the final crystal form. The relation between the two crystal forms deviates from the conventional concept of polymorphism by the similarity of the powder X-ray diffraction patterns. In principle, different mechanisms could be invoked to justify the different thermal behavior, amongst which a small difference in the crystal packing or in the atomic inter­actions could lead to it. The crystal structure of a monohydrate form of minodronic acid has been published earlier (Takeuchi et al., 1998 ▸), but neither the atomic coordinates nor details of the mol­ecular and crystal geometry were reported. In the absence of detailed information, the present study of the crystal structures of both forms, D and E, was undertaken with the aim of finding a reasonable solution to the unusual ‘polymorphism’ problem. Batches of the D and E crystals were prepared as described in the Experimental section and characterized by thermogravimetric analysis. Powder X-ray diffraction was performed with a Thermo X’tra Diffractometer equipped with a Si(Li) solid state detector directly on the as-prepared samples without grinding (Fig. 1 ▸). This confirmed that the diffraction patterns of the two forms involve peaks occurring at the same 2θ angles, the unique difference consisting of a different intensity distribution probably originating from preferential orientation. It is noteworthy that both diffraction patterns are fully compatible with the lattice parameters reported by Takeuchi et al. (1998 ▸). Small single crystals suitable for XRD experiments were selected from both batches and their structures solved and refined. The analyses established that the crystal structures are perfectly comparable within experimental error, and compatible with that previously reported.
Figure 1

Comparison of the room-temperature powder X-ray diffraction patterns of forms D (red line) and E (blue line) of the title compound.

We report herein the redetermination of the crystal structure of form D of minodronic acid monohydrate, whereas that of form E has been deposited at the CCDC (Rizzoli & Calestani, 2014 ▸). On the basis of the present study, the peculiar thermal behaviour of the D and E forms of minodronic acid cannot be ascribed to structural differences, but it is probably due to morphological or microstructural features induced by the crystallization procedure and the influence of the dehydration process.

Structural commentary

The asymmetric unit of the title compound, Fig. 2 ▸, consists of minodronic acid and a water mol­ecule of crystallization. The acid mol­ecule crystallizes in a zwitterionic form with cationic imidazolium[1,2a]pyridine and anionic phospho­nate groups. The fused-ring system deviates from planarity, with the dihedral angle formed by the pyridine and imidazole rings being 3.55 (9)°. An analysis of the bond lengths within the fused-ring system, indicates that the C=N and C=C double-bond distribution shown in the Scheme is the most probable, but a resonant form involving the aromatic character of the pyridine ring and the localization of the positive change on the atom N2 also exists. The values of the P—O bond lengths indicate that the negative charge is delocalized on atoms O2 and O3, whose distances [P1—O2 = 1.5197 (10), P1—O3 = 1.5001 (10) Å] are inter­mediate between those observed for the protonated atoms O4, O5 and O7 [mean value 1.564 (7) Å] and for the localized P=O double bond [P2—O6 = 1.4817 (10) Å]. The observed distribution of the hydrogen atoms and of the charge on the phospho­nate groups differs from that deducible from the mol­ecular plot reported previously by Takeuchi et al. (1998 ▸). An intra­molecular C—H⋯O hydrogen bond is present (Table 1 ▸).
Figure 2

The mol­ecular structure of the title compound, showing the atom labelling. Displacement ellipsoids are drawn at the 50% probability level.

Table 1

Hydrogen-bond geometry (, )

DHA DHHA D A DHA
O1H1OO80.78(3)1.87(3)2.6468(18)169(2)
O4H4OO2i 0.90(3)1.73(3)2.6116(16)169(3)
O5H5OO3ii 0.88(3)1.62(3)2.4973(14)175(3)
O7H7OO2iii 0.83(2)1.80(2)2.6183(14)170(2)
O8H81O6iv 0.82(2)1.99(2)2.8011(15)169(3)
O8H82O6v 0.91(2)1.83(2)2.7334(16)175(2)
N1H1NO5vi 0.85(2)2.18(3)2.9930(15)158(2)
C2H2BO20.972.503.0224(17)114
C4H4O80.932.403.2079(19)145
C8H8O4vii 0.932.513.1704(17)128

Symmetry codes: (i) ; (ii) ; (iii) ; (iv) ; (v) ; (vi) ; (vii) .

Supra­molecular features

The crystal structure of the title compound is characterized by a very effective hydrogen-bonding network which is responsible for the unusually high value of the calculated density of the crystal (1.818 g cm−3). As shown in Fig. 3 ▸, the crystal packing may be described as an alternate stacking along the b axis of phospho­nate and organic layers forming a three-dimensional network through O—H⋯O, N—H⋯O and C—H⋯O hydrogen bonds (Table 1 ▸ and Fig. 3 ▸). Adjacent phospho­nate layers are bridged via hydrogen bonds involving the water mol­ecules, which are hosted inside channels parallel to the a axis and running between pairs of fused-ring systems connected by π–π inter­actions [Cg1⋯Cg1i = 3.5822 (11) Å; Cg1 is the centroid of the N2/C5–C9 ring; symmetry code: (i) −x, −y + 1, −z + 1]. The O—H⋯O hydrogen-bonding system within the phospho­nate layer (Fig. 4 ▸) generates rings arranged in (8), (12), (16) and (16) graph-set motifs (Bernstein et al., 1995 ▸).
Figure 3

The crystal packing of the title compound viewed along the a axis, showing the hydrogen-bonding network (violet dashed lines) and π–π inter­actions (blue dashed lines); see Table 1 ▸ for details.

Figure 4

A view along the b axis of the hydrogen-bonding network (dashed lines) in the phospho­nate layer (see Table 1 ▸ for detail). Nitro­gen, carbon (except C1) and C- and N-bound H atoms have been omitted for clarity.

Synthesis and crystallization

Crystal form D: crystals of minodronic acid (150 g) were added to 5.6 l of 1 N hydro­chloric acid (about 37.5 volumes) in a reactor equipped with mechanical stirring, and dissolved with heating. The dissolution temperature (368 K) was maintained for at least 30 min, followed by filtration of a small amount of undissolved particles. The clear filtrate was stirred at 110 r.p.m. and slowly cooled to 328 K in 4 h, and then to 295 K overnight. The crystalline precipitate was collected by filtration, washed with 300 ml of water and 300 ml of ethanol, and dried at 318 K to obtain 135 g of pure form D of minodronic acid monohydrate, as shown by the dehydration peak at 415 K (Takeuchi et al., 1993 ▸). The colourless crystal sample used for X-ray analysis was selected from the batch. Crystal form E: crystals of minodronic acid (15 g) were added to 0.6 l of 1 N hydro­chloric acid (about 40 volumes) in a reactor equipped with magnetic stirring and dissolved with heating. The dissolution temperature (368 K) was maintained for at least 30 min followed by filtration of a small amount of undissolved particles. The clear filtrate was stirred slowly so that the liquid surface could not create a vortex (<110 r.p.m.) and gradually cooled down to 295 K overnight. The crystal precipitate was collected by filtration, washed with 30 ml of water and 30 ml of ethanol, and dried at 318 K to obtain 14 g of pure E form of minodronic acid monohydrate, as shown by the dehydration peak at 438 K (Takeuchi et al., 1993 ▸). The colourless crystal sample used for X-ray analysis was selected from the batch.

Refinement

Crystal data, data collection and structure refinement details are summarized in Table 2 ▸. All the H atoms were located in a difference Fourier map. The N- and O-bound hydrogen atoms were freely refined. The C-bound H atoms were refined using a riding model approximation, with C—H = 0.93–0.97 Å and with U iso(H) = 1.2U eq(C).
Table 2

Experimental details

Crystal data
Chemical formulaC9H12N2O7P2H2O
M r 340.16
Crystal system, space groupTriclinic, P
Temperature (K)294
a, b, c ()7.3668(11), 8.9833(12), 9.9733(18)
, , ()75.0136(17), 77.2716(17), 88.5706(18)
V (3)621.54(17)
Z 2
Radiation typeMo K
(mm1)0.40
Crystal size (mm)0.18 0.15 0.10
 
Data collection
DiffractometerBruker SMART 1000 CCD
Absorption correctionMulti-scan (SADABS; Bruker, 2008)
T min, T max 0.640, 0.746
No. of measured, independent and observed [I > 2(I)] reflections10459, 4054, 3413
R int 0.020
(sin /)max (1)0.752
 
Refinement
R[F 2 > 2(F 2)], wR(F 2), S 0.035, 0.101, 1.07
No. of reflections4054
No. of parameters218
H-atom treatmentH atoms treated by a mixture of independent and constrained refinement
max, min (e 3)0.55, 0.28

Computer programs: APEX2 and SAINT (Bruker, 2008 ▸), SIR97 (Altomare et al., 1999 ▸), SHELXL2014 (Sheldrick, 2008 ▸), ORTEP-3 for Windows (Farrugia, 2012 ▸), SCHAKAL99 (Keller, 1999 ▸) and publCIF (Westrip, 2010 ▸).

Crystal structure: contains datablock(s) I. DOI: 10.1107/S2056989014026863/su5037sup1.cif Structure factors: contains datablock(s) I. DOI: 10.1107/S2056989014026863/su5037Isup2.hkl Click here for additional data file. Supporting information file. DOI: 10.1107/S2056989014026863/su5037Isup3.cml CCDC reference: 1038087 Additional supporting information: crystallographic information; 3D view; checkCIF report
C9H12N2O7P2·H2OZ = 2
Mr = 340.16F(000) = 352
Triclinic, P1Dx = 1.818 Mg m3
a = 7.3668 (11) ÅMo Kα radiation, λ = 0.71073 Å
b = 8.9833 (12) ÅCell parameters from 227 reflections
c = 9.9733 (18) Åθ = 8.2–21.3°
α = 75.0136 (17)°µ = 0.40 mm1
β = 77.2716 (17)°T = 294 K
γ = 88.5706 (18)°Block, colourless
V = 621.54 (17) Å30.18 × 0.15 × 0.10 mm
Bruker SMART 1000 CCD diffractometer3413 reflections with I > 2σ(I)
ω scanRint = 0.020
Absorption correction: multi-scan (SADABS; Bruker, 2008)θmax = 32.3°, θmin = 2.2°
Tmin = 0.640, Tmax = 0.746h = −11→10
10459 measured reflectionsk = −12→13
4054 independent reflectionsl = −14→14
Refinement on F20 restraints
Least-squares matrix: fullHydrogen site location: mixed
R[F2 > 2σ(F2)] = 0.035H atoms treated by a mixture of independent and constrained refinement
wR(F2) = 0.101w = 1/[σ2(Fo2) + (0.060P)2 + 0.1062P] where P = (Fo2 + 2Fc2)/3
S = 1.07(Δ/σ)max = 0.001
4054 reflectionsΔρmax = 0.55 e Å3
218 parametersΔρmin = −0.28 e Å3
Geometry. All e.s.d.'s (except the e.s.d. in the dihedral angle between two l.s. planes) are estimated using the full covariance matrix. The cell e.s.d.'s are taken into account individually in the estimation of e.s.d.'s in distances, angles and torsion angles; correlations between e.s.d.'s in cell parameters are only used when they are defined by crystal symmetry. An approximate (isotropic) treatment of cell e.s.d.'s is used for estimating e.s.d.'s involving l.s. planes.
xyzUiso*/Ueq
P10.51328 (4)−0.00236 (4)0.28423 (3)0.01430 (9)
P20.17313 (4)0.13064 (4)0.15681 (3)0.01514 (9)
O10.49943 (14)0.25308 (12)0.08484 (10)0.0206 (2)
H1O0.576 (3)0.303 (3)0.098 (2)0.039 (6)*
O20.67125 (13)0.05492 (12)0.33555 (10)0.0208 (2)
O30.56833 (15)−0.09174 (12)0.17446 (11)0.0232 (2)
O40.37247 (14)−0.10822 (11)0.41416 (11)0.0214 (2)
H4O0.349 (4)−0.080 (3)0.496 (3)0.070 (8)*
O50.22204 (14)0.01195 (11)0.06740 (11)0.0213 (2)
H5O0.298 (4)0.044 (3)−0.016 (3)0.072 (8)*
O60.10307 (14)0.27806 (12)0.08260 (12)0.0258 (2)
O70.03415 (14)0.04399 (13)0.29628 (11)0.0240 (2)
H7O−0.079 (3)0.054 (3)0.300 (3)0.050 (7)*
O80.77583 (16)0.42957 (14)0.09502 (14)0.0301 (3)
H810.800 (3)0.515 (3)0.038 (3)0.049 (7)*
H820.884 (3)0.379 (3)0.086 (2)0.049 (6)*
N10.25482 (18)0.67430 (14)0.19284 (13)0.0225 (2)
H1N0.277 (3)0.769 (3)0.148 (2)0.047 (6)*
N20.09082 (15)0.46468 (12)0.31543 (12)0.0158 (2)
C10.38347 (17)0.16943 (14)0.21539 (13)0.0138 (2)
C20.34325 (19)0.26352 (15)0.32898 (14)0.0189 (2)
H2A0.25070.20610.40980.023*
H2B0.45650.27050.36190.023*
C30.27645 (18)0.42267 (15)0.28236 (14)0.0176 (2)
C40.3748 (2)0.55521 (16)0.20766 (15)0.0219 (3)
H40.50300.56370.17230.026*
C50.08150 (19)0.62025 (15)0.25700 (14)0.0182 (2)
C6−0.0888 (2)0.69424 (16)0.27288 (16)0.0236 (3)
H6−0.09650.79840.22950.028*
C7−0.2424 (2)0.60818 (17)0.35401 (17)0.0260 (3)
H7−0.35720.65410.36700.031*
C8−0.2293 (2)0.44988 (17)0.41874 (16)0.0240 (3)
H8−0.33500.39370.47660.029*
C9−0.06486 (19)0.37863 (15)0.39786 (15)0.0200 (3)
H9−0.05760.27370.43850.024*
U11U22U33U12U13U23
P10.01367 (15)0.01350 (15)0.01500 (15)0.00313 (11)−0.00231 (11)−0.00340 (11)
P20.01371 (15)0.01476 (16)0.01648 (16)0.00114 (11)−0.00367 (11)−0.00301 (11)
O10.0207 (5)0.0207 (5)0.0165 (4)−0.0069 (4)−0.0002 (4)−0.0004 (4)
O20.0134 (4)0.0274 (5)0.0221 (5)0.0026 (4)−0.0045 (4)−0.0072 (4)
O30.0304 (5)0.0192 (5)0.0202 (5)0.0075 (4)−0.0029 (4)−0.0084 (4)
O40.0239 (5)0.0183 (5)0.0185 (4)−0.0037 (4)−0.0016 (4)−0.0008 (4)
O50.0248 (5)0.0190 (5)0.0213 (5)−0.0009 (4)−0.0038 (4)−0.0084 (4)
O60.0243 (5)0.0199 (5)0.0327 (6)0.0059 (4)−0.0123 (4)−0.0016 (4)
O70.0138 (5)0.0324 (6)0.0216 (5)−0.0010 (4)−0.0005 (4)−0.0024 (4)
O80.0204 (5)0.0258 (6)0.0404 (7)−0.0019 (4)−0.0044 (5)−0.0036 (5)
N10.0259 (6)0.0131 (5)0.0244 (6)−0.0003 (4)0.0008 (5)−0.0028 (4)
N20.0186 (5)0.0112 (5)0.0179 (5)0.0024 (4)−0.0037 (4)−0.0044 (4)
C10.0132 (5)0.0131 (5)0.0142 (5)0.0004 (4)−0.0017 (4)−0.0028 (4)
C20.0235 (6)0.0156 (6)0.0192 (6)0.0063 (5)−0.0070 (5)−0.0060 (5)
C30.0188 (6)0.0149 (6)0.0198 (6)0.0035 (4)−0.0034 (5)−0.0070 (5)
C40.0206 (6)0.0193 (6)0.0247 (7)0.0010 (5)−0.0007 (5)−0.0077 (5)
C50.0236 (6)0.0121 (5)0.0183 (6)0.0021 (4)−0.0033 (5)−0.0040 (4)
C60.0277 (7)0.0164 (6)0.0273 (7)0.0077 (5)−0.0080 (6)−0.0059 (5)
C70.0221 (7)0.0264 (7)0.0324 (7)0.0079 (5)−0.0074 (6)−0.0121 (6)
C80.0197 (6)0.0240 (7)0.0271 (7)−0.0017 (5)−0.0015 (5)−0.0075 (5)
C90.0222 (6)0.0153 (6)0.0211 (6)−0.0018 (5)−0.0029 (5)−0.0038 (5)
P1—O31.5001 (10)N2—C91.3743 (17)
P1—O21.5197 (10)N2—C51.3744 (16)
P1—O41.5745 (10)N2—C31.4004 (17)
P1—C11.8534 (13)C1—C21.5558 (17)
P2—O61.4817 (10)C2—C31.4914 (18)
P2—O51.5468 (10)C2—H2A0.9700
P2—O71.5625 (10)C2—H2B0.9700
P2—C11.8426 (13)C3—C41.3580 (19)
O1—C11.4233 (15)C4—H40.9300
O1—H1O0.78 (2)C5—C61.4022 (19)
O4—H4O0.90 (3)C6—C71.360 (2)
O5—H5O0.88 (3)C6—H60.9300
O7—H7O0.83 (2)C7—C81.412 (2)
O8—H810.83 (3)C7—H70.9300
O8—H820.91 (2)C8—C91.356 (2)
N1—C51.3363 (18)C8—H80.9300
N1—C41.3700 (18)C9—H90.9300
N1—H1N0.85 (2)
O3—P1—O2116.20 (6)P2—C1—P1115.68 (7)
O3—P1—O4107.20 (6)C3—C2—C1116.39 (11)
O2—P1—O4109.85 (6)C3—C2—H2A108.2
O3—P1—C1110.32 (6)C1—C2—H2A108.2
O2—P1—C1106.68 (6)C3—C2—H2B108.2
O4—P1—C1106.17 (6)C1—C2—H2B108.2
O6—P2—O5115.07 (6)H2A—C2—H2B107.3
O6—P2—O7113.46 (6)C4—C3—N2105.72 (11)
O5—P2—O7105.64 (6)C4—C3—C2129.60 (13)
O6—P2—C1109.47 (6)N2—C3—C2124.62 (12)
O5—P2—C1107.81 (6)C3—C4—N1109.01 (13)
O7—P2—C1104.75 (6)C3—C4—H4125.5
C1—O1—H1O111.6 (16)N1—C4—H4125.5
P1—O4—H4O116.9 (17)N1—C5—N2107.39 (11)
P2—O5—H5O116.7 (18)N1—C5—C6131.56 (13)
P2—O7—H7O118.0 (17)N2—C5—C6121.04 (12)
H81—O8—H82105 (2)C7—C6—C5117.60 (13)
C5—N1—C4109.26 (12)C7—C6—H6121.2
C5—N1—H1N121.0 (15)C5—C6—H6121.2
C4—N1—H1N129.6 (15)C6—C7—C8120.72 (13)
C9—N2—C5120.93 (11)C6—C7—H7119.6
C9—N2—C3130.37 (11)C8—C7—H7119.6
C5—N2—C3108.59 (11)C9—C8—C7120.97 (13)
O1—C1—C2112.48 (10)C9—C8—H8119.5
O1—C1—P2101.50 (8)C7—C8—H8119.5
C2—C1—P2112.69 (8)C8—C9—N2118.63 (12)
O1—C1—P1106.38 (8)C8—C9—H9120.7
C2—C1—P1107.91 (8)N2—C9—H9120.7
D—H···AD—HH···AD···AD—H···A
O1—H1O···O80.78 (3)1.87 (3)2.6468 (18)169 (2)
O4—H4O···O2i0.90 (3)1.73 (3)2.6116 (16)169 (3)
O5—H5O···O3ii0.88 (3)1.62 (3)2.4973 (14)175 (3)
O7—H7O···O2iii0.83 (2)1.80 (2)2.6183 (14)170 (2)
O8—H81···O6iv0.82 (2)1.99 (2)2.8011 (15)169 (3)
O8—H82···O6v0.91 (2)1.83 (2)2.7334 (16)175 (2)
N1—H1N···O5vi0.85 (2)2.18 (3)2.9930 (15)158 (2)
C2—H2B···O20.972.503.0224 (17)114
C4—H4···O80.932.403.2079 (19)145
C8—H8···O4vii0.932.513.1704 (17)128
  6 in total

1.  A short history of SHELX.

Authors:  George M Sheldrick
Journal:  Acta Crystallogr A       Date:  2007-12-21       Impact factor: 2.290

2.  Effect of minodronic acid (ONO-5920) on bone mineral density and arthritis in adult rats with collagen-induced arthritis.

Authors:  Itsuro Yamane; Hiroshi Hagino; Toru Okano; Makoto Enokida; Daisuke Yamasaki; Ryota Teshima
Journal:  Arthritis Rheum       Date:  2003-06

3.  Effect of minodronic acid hydrate on hip geometry in Japanese women with postmenopausal osteoporosis.

Authors:  Masako Ito; Teruki Sone; Masao Fukunaga
Journal:  J Bone Miner Metab       Date:  2009-11-25       Impact factor: 2.626

4.  Studies on novel bone resorption inhibitors. II. Synthesis and pharmacological activities of fused aza-heteroarylbisphosphonate derivatives.

Authors:  M Takeuchi; S Sakamoto; K Kawamuki; H Kurihara; H Nakahara; Y Isomura
Journal:  Chem Pharm Bull (Tokyo)       Date:  1998-11       Impact factor: 1.645

5.  A third-generation bisphosphonate, minodronic acid (YM529), successfully prevented the growth of bladder cancer in vitro and in vivo.

Authors:  K Sato; T Yuasa; M Nogawa; S Kimura; H Segawa; A Yokota; T Maekawa
Journal:  Br J Cancer       Date:  2006-10-17       Impact factor: 7.640

Review 6.  A review of minodronic acid hydrate for the treatment of osteoporosis.

Authors:  Shinji Tanishima; Yasuo Morio
Journal:  Clin Interv Aging       Date:  2013-02-15       Impact factor: 4.458

  6 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.