Literature DB >> 19223961

Reaction of chromium(III) with 3,4-dihydroxybenzoic acid: kinetics and mechanism in weak acidic aqueous solutions.

Kimon Zavitsanos1, Konstantinos Tampouris, Athinoula L Petrou.   

Abstract

The interactions between chromium(III) and 3,4-dihydroxybenzoic acid (3,4-DHBA) were studied resulting in the formation of oxygen-bonded complexes upon substitution of water molecules in the chromium(III) coordination sphere. The experimental results show that the reaction takes place in at least three stages, involving various intermediates. The first stage was found to be linearly dependent on ligand concentration k(1(obs))' = k(0) + k(1(obs))[3, 4-DHBA], and the corresponding activation parameters were calculated as follows: DeltaH(1(obs)) ( not equal) = 51.2 +/- 11.5 kJ mol(-1), DeltaS(1(obs)) ( not equal) = -97.3 +/- 28.9 J mol(-1) K(-1) (composite activation parameters) . The second and third stages, which are kinetically indistinguishable, do not depend on the concentrations of ligand and chromium(III), accounting for isomerization and chelation processes, respectively. The corresponding activation parameters are DeltaH(2(obs)) ( not equal) = 44.5 +/- 5.0 kJ mol(-1), DeltaS(2(obs)) ( not equal) = -175.8 +/- 70.3 J mol(-1) K(-1). The observed stages are proposed to proceed via interchange dissociative (I(d), first stage) and associative (second and third stages) mechanisms. The reactions are accompanied by proton release, as is shown by the pH decrease.

Entities:  

Year:  2009        PMID: 19223961      PMCID: PMC2635946          DOI: 10.1155/2008/212461

Source DB:  PubMed          Journal:  Bioinorg Chem Appl            Impact factor:   7.778


1. INTRODUCTION

The ligand 3,4-dihydroxybenzoic acid (3,4-DHBA) (Figure 1), 1, is known to be produced in the reaction of radicals, which are formed in pathophysiological cases (e.g., ischemic stroke [1, 2], traumatic brain injury [3], and Huntington's disease [4]), the 4-hydroxybenzoic acid or salicylic acid acting as radical-trapping agent [5]. It is used for assisting the Fenton reaction in effluent treatment [6] and in dechlorination of polychlorinated dioxins [7].
Figure 1

3,4-DHBA.

3,4-DHBA is employed in the preparation of resins having antioxidant properties, acting as radical scavenger [8] as well as in the preparation of composite polymer modified electrodes [9]. It was also found to interfere in the nucleation and crystal growth of radial alumina trihydrate particles [10]. Phenolic acids, in general, are present in fruits and plants and are participating in the chemical structure of humic substances, which can coordinate with nutrient ions, especially metal ions, increasing thus their bioavailability [11-14]. 3,4-dihydroxybenzoic acid, in particular, bearing both catecholic and carboxylic sites, shows special complexing properties. Cr(III), although present in traces in biological systems, is considered to play a role in the activation of some enzymes [15]. It has also been identified as a reactive component of an oligopeptide known as low-molecular-weight chromium binding substance or chromoduline [16, 17]. Yet, the biological role of Cr(III) remains mostly unclarified, one of the main reasons being the lack of intense characteristics like charge-transfer bands in the spectra with the only exception being the organochromium complexes [18, 19]. In the present study, the reaction of Cr(III) and 3,4-DHBA in weak acidic aqueous solutions is investigated. The mechanism of the complex formation and the stability of the complexes formed are studied and presented, in order to contribute in the clarification of the role of complexation in the uptake of metals in various biological systems (e.g., plants).

2. EXPERIMENTAL RESULTS

2.1. Reagents and materials

All reagents employed were of analytical grade and were used as received. Aqueous solutions containing 3,4-DHBA (Alfa Aesar) in concentrations ranging from 7.45 × 10−3 to 1.62 × 10−2 M were prepared using dilute (0.1 M) KOH solutions for pH adjustment in order for the ligand to be dissolved. Stock solutions (0.2 M) of Cr(III) were prepared from Cr(NO3)2·9H2O. The ionic strength was adjusted using KNO3. All solutions used in the present study were freshly prepared in order to avoid side reactions (transformation and decomposition). The addition of the Cr(III) solution kept the pH below 4 due to its acidic hydrolysis.

2.2. Kinetic experiments

Electronic spectra were recorded on a Varian Cary 3E, UV-vis spectrophotometer. The kinetic experiments were also followed at the above instrument. All kinetic experiments were performed at pH values below 4 in the presence of air. Pseudo-first-order conditions were employed for most of the kinetic experiments. For the first stage of the reactions which was studied at 279–295 K, the duration of the step was ~16 000–5000 seconds, the concentrations were of the order of 10−2 M for the ligand and 10−3 M for Cr(III), and plots of ln(A − A ∞) against time (the absorbance is decreasing) were set up, where A and A ∞ are absorbances at time t and at infinite time (after the completion of this step). The plots were found to be linear (Figure 2) for at least three half lives. The rate constants were calculated from the constant slope of the line. At the higher temperatures where the second and third stages were studied, the first stage being faster is not observed. For the second and third (consecutive) stages which were studied at 303–323K, the duration of the steps was ~20 000–8000 seconds, the concentrations were of the order of 10−2 M for the ligand and 10−3 M for Cr(III), and plots of ln(A ∞ − A ) against time (the absorbance is increasing) were also found to be linear. This is the case of a polyfunctional compound acting sequentially; a linear plot of ln(A ∞ − A ) = f (t) could be obtained when a biphasic reaction takes place resembling to a single-stage reaction [20]. This could also happen when k 3 ≫ k 2 or k 2 ≫ k 3. In the case of our system, k 2 could be a lot smaller than k 3. The k 1 and k 2 (k 2 ≪ k 3) values at various temperatures are given in Table 1.
Figure 2

Typical linear plot of ln(A − A ∞) versus time at 279 K and [3,4-DHBA] = 1.50 × 10−2 M, [Cr(III)] = 2.50 × 10−3 M, I = 0.02 M.

Table 1

Values of k 1(obs) and k 2 (k 2 ≪ k 3) at various temperatures.

Temperature (K)Rate constants
k 1(obs) × 102 (M−1 s−1) k 2 × 104 (s−1)
2791.10
2862.52
2903.15
2953.74
3030.91
3091.35
3131.54
3182.11
3233.11
The kinetics was followed at various wavelengths yielding identical results though changes in absorbance were in all cases small. Uncomplexed Cr(III) species does not interfere in the (absorbance) measurements since it is included in both A and A ∞ and is thus eliminated. The activation parameters (ΔH and ΔS ) were calculated from the linear Eyring plots according to activated complex theory. The activation parameters ΔH 1(obs) and ΔS 1(obs) corresponding to k 1(obs) and ΔH 2(obs) , ΔS 2(obs) corresponding to k 2 (k 2 ≪ k 3) were thus estimated and presented in Table 2.
Table 2

Activation parameters for steps 1 (k 1(obs)) and 2 (k 2 ≪ k 3).

ΔH 1(obs) (kJ mol−1)ΔS 1(obs) (J mol−1 K−1)ΔH 2(obs) (kJ mol−1)ΔS 2(obs) (J mol−1 K−1)
51.2 ± 11.5 −97.3 ± 28.944.5 ± 5.0 −175.8 ± 70.3
The A ∞ values were obtained from the kinetic measurements and from the A = f (t) plots (at the certain wavelengths) making it possible to check if the reaction was run to completion.

3. DISCUSSION

3.1. Kinetics and mechanisms

The UV-vis spectra of Cr(III) and 3,4-DHBA solutions of concentrations 5.5 × 10−3 M are presented in Figure 3. At pH < 4, where all the kinetic experiments were conducted, Cr(III) exists mainly in the hexaaqua monomeric form. However, reaction with Cr(H2O)5(OH)2+ should be considered at the pH range 3 to 4, since a small amount of Cr(H2O)5(OH)2+ is present due to the equilibrium which is characterized by a pKa value (Cr3+/Cr(OH)2+) of about 4 [21]. The visible spectrum of Cr(H2O)6 3+ exhibits two maxima in the region of 410 and 575 nm (Figure 3) accounting for the 4A2g → 4T1g and the 4A2u → 4T2g Cr(III) transitions, respectively.
Figure 3

UV-vis spectra of Cr(III) (spectrum 1) and 3,4-DHBA (spectrum 2) solutions of 5.5 × 10−3 M concentrations and spectrophotometric cell path d = 1 cm.

Under the applied experimental conditions, the ligand exists as neutral molecule and monoanion [22, 23] abbreviated as DHBA and DHBA−, respectively. In the ligand molecules, intramolecular hydrogen bonding between the two hydroxyl groups occurs, favored by the formation of a five-membered ring. The acid ligand and the Cr(H2O)6 3+ complex dissociation equilibria established are as follows: The corresponding pK values for the above acid dissociation constants, at 25°C, are given pK1 = 4.5, pK2 = 8.7, pK3 = 12.8, and pKa = 4.0. Upon mixing of the reactants, violet Cr(III) and light-brown ligand solutions, formation of a light green complex, 2, assigned to be oxygen-bound Cr(III) compound takes place. The substitution of water molecules in the Cr(III) coordination sphere by 3,4-DHBA molecule results in a change of the ligand field. Decrease in absorbance at 575 nm is at first observed. The formation of 2 is consistent with the formation and subsequent transformation (substitution) kinetics. Experiments conducted at various temperatures (279, 286, 290, and 295 K) yielded identical kinetic behavior. However, the attacks by Cr(III) can take place only by releasing protons because the hydroxyls and the carboxylic groups are efficiently blocked (are protonated). This results in the pH decrease of the solution (Figure 4). Thus, the k 1 pathway occurs by attack of Cr(H2O)5OH2+ at the carboxylic ligand group leading to a complex which, in acidic solutions, is protonated, 2.
Figure 4

Typical pH versus time plot of 3,4-DHBA/Cr(III) mixtures at 313 K. Conditions: [3,4-DHBA]0 = 1.62 × 10−2 M, [Cr(III)]0 = 2.50 × 10−3 M, [KOH] = 4.00 × 10−3 M, I = 0.02 M.

The kinetics of the first stage of the reaction (k 1(obs)′) followed a first-order rate law, k 1(obs)′ = k 0 + k 1(obs)[3, 4-DHBA]. The rate constant k 0 corresponds to a reaction of Cr(III) when [3,4-DHBA] = 0.0, that is, to a reaction which does not involve the ligand. This may be any reaction of Cr(III) with KOH or with other Cr(III) ions, that is, an oligomerization reaction which always takes place in Cr(III) aqueous solutions. This does not affect the reactions and the mechanism of the reaction with 3,4-DHBA. That is, in order to obtain information on how many ligands react in the first step, the dependence of k 1(obs)′ on the ligand concentration was studied (Figure 5). In previous studies [24, 25], where the reactions of Cr(III) with 3,4-dihydroxyphenylpropionic acid (dihydrocaffeic acid) and 3,4-dihydroxyphenylpropenoic acid (caffeic acid) were investigated kinetically, the experimental data led to the conclusion that the reactive form of the metal ion is Cr(OH)2+ and the reactive form of the ligand is the neutral species.
Figure 5

Dependence of k 1(obs)′ on ligand concentration at 279 K and [Cr(III)] = 2.50 × 10−3 M, I = 0.02 M.

A dissociative mechanism, I is expected for step 1 since, in the case of the conjugate base Cr(H2O)5OH2+, the dissociative mechanism I is favored [26]. This is because of the strong labilizing effect, which is induced by the coordinated OH−, presumably, on the trans H2O molecule. This leads to a 102-103 fold enhanced reaction rate for the hydroxy-aqua over the hexaaqua ion. In conclusion, for the first attack, a reaction of Cr(OH)2+ can be proposed (Scheme 1). This would suggest an inverse dependence on [H+] since [Cr(OH)2+] = Ka[Cr(aq) 3 +]/[H+], which is not observed, implying that a protonation of the hydrogen-bonded ligand molecule exists; a breakout of the hydrogen bonding must take place before the attack by Cr(OH)2+. We suggest attack on the –OH of the protonated –COOH group (–HOC=OH+) since there the proton is loosely attached to oxygen due to the fact that withdrawing of negative charge by the H+ of the C=OH+ is taking place. The protonation of the –COOH group though in small extent makes the H of the OH more labile and so the attack by Cr(H2O)5(OH)2+ is easier. The small extent of the protonation of the –COOH group is the reason for the small yield of the reaction which is expressed by the small absorbance change ΔA.
Scheme 1
Scheme 1 presents the mechanism proposed, according to the above experimental results. Therefore, it holds that Rate = k 1[Cr(OH)2+][3, 4-DHBAH+] = k 1K0[3, 4-DHBA][H+](Ka[Cr(aq) 3 +]/[H+]) = k 1K0Ka[3, 4-DHBA][Cr(aq) 3 +] = k 1(obs)[3, 4-DHBA][Cr(aq) 3 +], where k 1(obs) = k 1K0Ka. In Figure 6, spectra of the reaction mixture at 288 K are recorded at various times after mixing. The decrease in absorbance at short reaction times and the increase in longer reaction times suggest that the reaction is a multistep process. Spectra are shown only for the decrease in absorbance since at longer reaction time the increase in absorbance that follows would cause overlapping of the spectra.
Figure 6

Typical UV-vis spectra of 3,4-DHBA/Cr(III) mixture at various times after mixing at 288 K. Conditions: [3,4-DHBA] = 6.80 × 10−3 M, [Cr(III)] = 3.75 × 10−2 M.

The calculated negative value of ΔS 1(obs) would suggest an associative mechanism for the first stage of the reaction. Cr(H2O)5OH2+, however, as stated above, bears a dissociative, (I ), mechanism. This suggests, as presented above, that a complex reaction is taking place and composite activation parameters exist, that is, ΔH 1(obs) = ΔH 0 + ΔH + ΔH 1 and ΔS 1(obs) = ΔS 0 + ΔS + ΔS 1 , where ΔH 0 and ΔS 0 correspond to an equilibrium established prior to step 1. ΔH 1 and ΔS 1 correspond to the first step (k 1). Therefore, the resulting negative value of ΔS as well as the resulting value of ΔH do not correspond only to step 1 reaction which is actually taking place by an I mechanism due to the reactive species Cr(H2O)5OH2+. Dependence on ligand and Cr(III) concentrations was studied in order to find if a second or third ligand molecule or Cr(III) species enters the coordination sphere of the already formed complex 2. Therefore, the possibility of two consecutive steps involving transformations of 2 to 3 and of 3 to 4 is investigated. The above scheme should give a biphasic reaction resembling a single-stage reaction when a polyfunctional compound reacts sequentially [20]. Thus, in our case, the polyfunctional ligand 3,4-DHBA reacting with Cr(III) gives the polyfunctional compound 2 which reacts sequentially and the biphasic reaction resembles a single-stage reaction giving a linear plot of ln(A ∞ − A ) = f (t). This could also happen when k 3 ≫ k 2. The experimental results (Figure 7) clearly demonstrate that the concentrations of 3,4-DHBA and Cr(III) have no effect on the observed rate constants k 2(obs) (≪k 3(obs)) in the concentration range studied. However, study at higher concentrations was restricted due to the low solubility of 3,4-DHBA in weak acidic aqueous solutions. Values of k 2(obs) (≪k 3(obs)) are presented in Table 1 and the corresponding ΔH 2(obs) , ΔS 2(obs) values are listed in Table 2.
Figure 7

Dependence of k 2(obs) upon (a) ligand concentration at [Cr(III)] = 2.50 × 10−3 M, I = 0.02 M, 313 K and (b) Cr(III) concentration at [3,4-DHBA] = 7.45 × 10−3 M, I = 0.5 M, 318 K.

The fact that the second and third steps (k 2, k 3) were found to be independent both on ligand and Cr(III) concentrations, that is, the rate exhibits a first-order dependence on the product of step 1, suggests that transformations are taking place within the already formed complex, 2. The activation parameters deduced from the temperature-dependence experiments may be used for proposing structures of the activated complexes (Scheme 2), and the taking place mechanism.
Scheme 2

Activated complexes 2 and 3 .

The negative value of ΔS 2 , the independence of k 2(obs) (k 2 ≪ k 3) on both ligand and Cr(III) concentrations, and the increase in absorbance (i.e., of the extinction coefficients) led to the assignment of the observed transformations as associatively activated substitution of water molecules from the Cr(III) coordination sphere by 3,4-DHBA through isomerization and chelation in two consecutive steps (). The pH decrease suggests a concomitant proton release as proposed in the mechanism (Scheme 1). Scheme 1 shows that 2, the first formed species, undergoes isomerization to give another oxygen-bonded species, 3. Complex 3 chelates in a k 3 step, (k 2 ≪ k 3) chelation, to produce 4, the final chelated complex [Cr(3,4-DHBA−3H)(H2O)4]. Isolation in the solid form of the final product gives a compound the elemental analyses of which (15.4% C and 3.7% H) correspond to the formula [Cr(3,4-DHBA−3H)(H2O)4]·2KNO3·4H2O, for which the calculated percentages for C and H are 15.3% and 3.5%, respectively. Alternatively, the first attack could occur at the hydroxyl group (step k 1) and a following chelation could account for a second step (k 2) only. We prefer though the already suggested mechanism which is in accord with the mechanisms of the reactions of Cr(III) with 3,4-dihydroxyphenylpropionic acid and 3,4-dihydroxyphenylpropenoic acid where the three stages (i.e., complexation, isomerization, and chelation) were kinetically distinguishable [24, 25]. In our present case, three steps are also suggested because of the necessity for isomerization prior to chelation. Chelation between the groups carboxylic and phenolic is unfavorable due to the formation of seven-membered ring. We exclude the possibility of a reversible reaction, where the final absorbance could correspond to Cr(III) and complex (in equilibrium) resulting in k obs = k + k based on the inertness of the various species (small values of reaction rates) and the existence of further reactions which, then, would shift an equilibrium well to the right. Another possibility is attack at the carboxylic group and subsequent chelation at the same group through the oxygens. We are discarding this alternative because a four-membered ring would result if chelation at the carboxylic group took place whereas a more preferable situation, that is, a five-membered ring is formed according to the suggested mechanism. Since Cr, Mo, and W are in the same group of the periodic table, the following comparison can be made: [Cr(H2O)5(OH)]2+ in our system reacts with a rate constant two orders of magnitude slower than that of [W3O4(H2O)9]4+ and three orders of magnitude slower than that of [Mo3O4(H2O)9]4+ indicating an extremely less labile species [27]. The comparison is made for substitution reactions of the conjugate base forms of the above aqua ions.

3.2. Structure of the complexes—mode of binding

In the UV-vis spectra, the changes are caused only by the change of the ligand field. Since oxidation does not take place under the experimental conditions, only a shift of the maximum and disappearance of the shoulder due to the high absorbtivities in the UV region of the produced complexes take place. A 1 : 1 stoichiometry for the reaction of Cr(III) with 3,4-DHBA is proposed because of the observed k 1 dependence on ligand concentration. The elemental analyses of the isolated final product also suggest a 1 : 1 stoichiometry. The consecutive isomerization and chelation reactions taking place in the Cr(III) center do not cause any change in stoichiometry. The proposed structures of the activated complexes 2 and 3 are given in Scheme 2. Associative mechanism has been supported to operate in reactions of Cr(III) [18, 28, 29]. In the above suggested mechanism, the phenolic groups act as internal attacking groups to the H2O molecules of the Cr(III) coordination sphere supplying thus a proton, which is then released as H3O+ (Scheme 1). The driving force for the production of the final chelated product 4 must be the most stabilized chelated form compared to forms 2 or 3. The suggested catecholic mode of binding was also found to operate in the coordination complexes of 3,4-dihydroxyphenylpropionic acid (dihydrocaffeic acid) and 3,4-dihydroxyphenylpropenoic acid (caffeic acid) with Cr(III) [24, 25] as already mentioned. This type of binding was also reported for complexes of dihydrocaffeic, caffeic, and ferulic acids with Co(II), Ni(II), Cu(II), Fe(III), Mn(II), Mn(III), V(V), V(IV,V), and Zn(II) [22, 23, 30–32]. Catecholic type of coordination was also suggested for the Fe(III)-2,3-DHBA complex [33, 34].

4. CONCLUSIONS

In the present study, the reaction between Cr(III) and 3,4-DHBA in weak acidic aqueous solutions was investigated. The experimental results are consistent with a three-step mechanism in which an initial attack (step 1) between the acid molecule (ligand) and the Cr(H2O)5OH2+ complex giving a carboxylate bound Cr(III), 2, is followed by two consecutive kinetically indistinguishable nonligand and non-Cr(III) dependent steps. The two consecutive steps are assigned as isomerization and chelation step (steps 2 and 3). The reactions are followed by a pH decrease because proton release is taking place, according to an associative mode of activation (steps 2 and 3). The negative value of the entropy of activation of step 2, the independence on ligand and Cr(III) concentrations, the increase of the extinction coefficients, and the pH decrease due to release of protons upon complexation led to the proposed mechanism (Scheme 1). The observed transformations were assigned as substitution of water molecules from the coordination sphere of Cr(III) by the ligand through complexation (step 1) following an I mechanism, isomerization, and chelation in two consecutive kinetically indistinguishable steps supported to follow associative mechanisms (step 2 and step 3).
  10 in total

1.  Isolation and characterization of a biologically active chromium oligopeptide from bovine liver.

Authors:  C M Davis; J B Vincent
Journal:  Arch Biochem Biophys       Date:  1997-03-15       Impact factor: 4.013

2.  Monitoring of reactive oxygen species production after traumatic brain injury in rats with microdialysis and the 4-hydroxybenzoic acid trapping method.

Authors:  N Marklund; F Clausen; T Lewander; L Hillered
Journal:  J Neurotrauma       Date:  2001-11       Impact factor: 5.269

3.  Hydroxyl radical formation is greater in striatal core than in penumbra in a rat model of ischemic stroke.

Authors:  Shimin Liu; Miao Liu; Steven Peterson; Minoru Miyake; Val Vallyathan; Ke Jian Liu
Journal:  J Neurosci Res       Date:  2003-03-15       Impact factor: 4.164

4.  On the application of 4-hydroxybenzoic acid as a trapping agent to study hydroxyl radical generation during cerebral ischemia and reperfusion.

Authors:  Miao Liu; Shimin Liu; Steven L Peterson; Minoru Miyake; Ke Jian Liu
Journal:  Mol Cell Biochem       Date:  2002 May-Jun       Impact factor: 3.396

5.  Kinetics and mechanism of the reaction between chromium(III) and 3,4-dihydroxyphenylpropionic (dihydrocaffeic) acid in weak acidic aqueous solutions.

Authors:  Athinoula L Petrou; Patrina Paraskevopoulou; Maria Chrysikopoulou
Journal:  J Inorg Biochem       Date:  2004-01       Impact factor: 4.155

6.  Isolation of a biologically active low-molecular-mass chromium compound from rabbit liver.

Authors:  A Yamamoto; O Wada; T Ono
Journal:  Eur J Biochem       Date:  1987-06-15

7.  Increased vulnerability to 3-nitropropionic acid in an animal model of Huntington's disease.

Authors:  M B Bogdanov; R J Ferrante; S Kuemmerle; P Klivenyi; M F Beal
Journal:  J Neurochem       Date:  1998-12       Impact factor: 5.372

8.  The electrochemical copolymerization of 3,4-dihydroxybenzoic acid and aniline at microdisk gold electrode and its amperometric determination for ascorbic acid.

Authors:  J J Sun; D M Zhou; H Q Fang; H Y Chen
Journal:  Talanta       Date:  1998-03       Impact factor: 6.057

9.  An in vitro hydroxyl radical generation assay for microdialysis sampling calibration.

Authors:  Rui Chen; Julie A Stenken
Journal:  Anal Biochem       Date:  2002-07-01       Impact factor: 3.365

10.  Kinetics and mechanism of the reaction between chromium(III) and 3,4-dihydroxy-phenyl-propenoic acid (caffeic acid) in weak acidic aqueous solutions.

Authors:  Vladimiros Thoma; Konstantinos Tampouris; Athinoula L Petrou
Journal:  Bioinorg Chem Appl       Date:  2008       Impact factor: 7.778

  10 in total
  2 in total

1.  Kinetics and Mechanisms of the Chromium(III) Reactions with 2,4- and 2,5-Dihydroxybenzoic Acids in Weak Acidic Aqueous Solutions.

Authors:  Kimon Zavitsanos; Athinoula L Petrou
Journal:  Bioinorg Chem Appl       Date:  2010-07-11       Impact factor: 7.778

2.  Kinetics and Mechanism of the Reaction between Chromium(III) and 2,3-Dihydroxybenzoic Acid in Weak Acidic Aqueous Solutions.

Authors:  Athinoula L Petrou; Vladimiros Thoma; Konstantinos Tampouris
Journal:  Bioinorg Chem Appl       Date:  2010-05-31       Impact factor: 7.778

  2 in total

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