Literature DB >> 25141257

Pharmacokinetics-pharmacodynamics analysis of bicyclic 4-nitroimidazole analogs in a murine model of tuberculosis.

Suresh B Lakshminarayana1, Helena I M Boshoff2, Joseph Cherian1, Sindhu Ravindran1, Anne Goh1, Jan Jiricek1, Mahesh Nanjundappa1, Amit Nayyar2, Meera Gurumurthy1, Ramandeep Singh2, Thomas Dick1, Francesca Blasco1, Clifton E Barry2, Paul C Ho3, Ujjini H Manjunatha1.   

Abstract

PA-824 is a bicyclic 4-nitroimidazole, currently in phase II clinical trials for the treatment of tuberculosis. Dose fractionation pharmacokinetic-pharmacodynamic studies in mice indicated that the driver of PA-824 in vivo efficacy is the time during which the free drug concentrations in plasma are above the MIC (fT>MIC). In this study, a panel of closely related potent bicyclic 4-nitroimidazoles was profiled in both in vivo PK and efficacy studies. In an established murine TB model, the efficacy of diverse nitroimidazole analogs ranged between 0.5 and 2.3 log CFU reduction compared to untreated controls. Further, a retrospective analysis was performed for a set of seven nitroimidazole analogs to identify the PK parameters that correlate with in vivo efficacy. Our findings show that the in vivo efficacy of bicyclic 4-nitroimidazoles correlated better with lung PK than with plasma PK. Further, nitroimidazole analogs with moderate-to-high volume of distribution and Lung to plasma ratios of >2 showed good efficacy. Among all the PK-PD indices, total lung T>MIC correlated the best with in vivo efficacy (rs = 0.88) followed by lung Cmax/MIC and AUC/MIC. Thus, lung drug distribution studies could potentially be exploited to guide the selection of compounds for efficacy studies, thereby accelerating the drug discovery efforts in finding new nitroimidazole analogs.

Entities:  

Mesh:

Substances:

Year:  2014        PMID: 25141257      PMCID: PMC4139342          DOI: 10.1371/journal.pone.0105222

Source DB:  PubMed          Journal:  PLoS One        ISSN: 1932-6203            Impact factor:   3.240


Introduction

Every year nearly 8 million new cases of Tuberculosis (TB) are reported globally resulting in 1.4 million deaths [1]. Poor treatment compliance – due to the requirement for prolonged multidrug therapy – as well as the use of inadequate regimens has fueled the emergence of multi-drug-resistant and extensively-drug-resistant (MTD-TB and XDR-TB) TB strains. MDR-TB is resistant to at least isoniazid and rifampicin and XDR-TB is resistant to isoniazid, rifampicin, fluoroquinolones and at least one of the injectables [2]. TB control programs are further complicated in settings where the incidence of co-infection with HIV is high because drug-drug interactions with anti-retroviral therapy are difficult to avoid [2], [3]. Hence there is an urgent need to discover new TB drugs active against drug-resistant forms of TB and compatible with treatment against HIV. PA-824 [4] and OPC-67683 [5] are two bicyclic 4-nitroimidazoles currently in phase II clinical trials, representing a promising new class of therapeutics for TB [6]. Preclinical testing of PA-824 showed bactericidal activity in various in vitro and in vivo models [7], . PA-824 was shown to be well tolerated in healthy subjects, following oral daily doses for 7 days [9]. These results, combined with the demonstrated activity of PA-824 against drug-sensitive and multidrug-resistant Mtb, supported the progression of this compound and its evaluation as a novel treatment against TB. An early bactericidal activity (EBA) study of PA-824 revealed a lack of dose response between 200 and 1200 mg administered daily for 14 days [10]. Dose-fractionation PK-PD studies in mice showed the PK-PD driver of PA-824 to be the time during which the free drug concentrations in plasma were above the MIC (fT) [11]. In retrospect, clinical investigators established that fT was 100% at all doses between 200 and 1200 mg daily. An additional phase II trial between 50 and 200 mg was undertaken to establish the lowest efficacious dose [12]. 200 mg of PA-824 was found to be efficacious and used in combination with other anti-TB drugs [13]. Physico-chemical properties, in vitro potency, in vitro and in vivo pharmacokinetics (PK) are critical determinants for in vivo efficacy. PA-824 is highly lipophilic and exhibits poor aqueous solubility. To overcome the limitation of its low solubility and improve its oral bioavailability, a cyclodextrin formulation was developed and used for in vivo animal efficacy studies [4], [7]. Extensive lead optimization efforts were undertaken to improve aqueous solubility, metabolic stability, in vitro potency and in vivo efficacy of anti-tubercular nitroimidazoles and various analogs were synthesized [5], [14]–[23]. Comprehensive in vivo pharmacology studies are generally resource and time intensive. This is particularly true for TB because of the slow rate of M. tuberculosis growth, lengthy treatment duration and requirement of high containment facility. In this study, a panel of closely related potent bicyclic 4-nitroimidazoles (NI) was profiled both in vitro and in vivo. The data is retrospectively analyzed to identify the PK parameters that correlate with in vivo efficacy of a series of bicyclic 4-nitroimidazoles. The results of this investigation showed that PK properties such as volume of distribution and lung exposure predicts in vivo efficacy of bicyclic 4-nitroimidazoes better than other PK parameters. Thus, in vitro potency and lung PK could be used as surrogate to guide the prioritization of new pre-clinical candidates for lengthy efficacy studies, there by expediting drug discovery.

Materials and Methods

Bacteria, culture conditions and chemicals

M. tuberculosis (Mtb) (H37Rv, ATCC 27294) culture conditions have been described previously [16]. Synthesis of PA-824, Amino-824, Aminoethyl-824, NI-135, NI-147 and NI-136 have been previously reported [16], [20]. Other NI analogs NI-622, NI-644, NI-176, NI-269, NI-182, NI-145, NI-297 and NI-302 have been described in two patents [14], [18] and synthesis of these compounds to be described elsewhere. All solutions were made as 20 mM stocks in DMSO. Hydroxypropyl-β-cyclodextrin, Lecithin, granular was purchased from Acros/Organics, New Jersey, and USA. Minimum Inhibitory concentration (MIC99) against wild type Mtb H37Rv and cofactor F420 deficient (FbiC mutant) [24] was determined by the broth dilution method as described earlier [16].

In vitro physico-chemical properties

In vitro physicochemical and PK parameters like solubility, log P, PAMPA, Caco-2 permeability and mouse plasma protein binding were determined in-house in medium to high throughput format assays. Briefly, solubility was measured using a high throughput equilibrium-solubility (HT-Eq sol) assay using a novel miniaturized shake-flask approach and streamlined HPLC analysis [25]; lipophilicity determination was carried out in 96-well micro titer plates and the diffusion of compounds between two aqueous compartments separated by a thin octanol liquid layer was measured [26]; PAMPA permeability experiments were carried out in 96-well micro titer filter plates at absorption wavelengths between 260 and 290 nm [27]; Caco-2 assay was carried out in a 96-well format, and compound concentrations in each chamber were measured by LC/MS as described previously [28] and plasma protein binding was determined in mouse plasma using an ultra-filtration method [29].

Ethics Statement

All animal experimental protocols (protocol #023/2009 and #025/2009 for PK; protocol #004/2010 for efficacy) involving mice were reviewed and approved by the Institutional Animal Care and Use Committee (IACUC) of Novartis Institute for Tropical Diseases (NITD). The animal research complied with Singapore Animal Veterinary Authority and global Novartis policy on the care and use of animals. Experimental and control animals infected with Mtb were euthanized at the end of the experiment. All procedures during pharmacokinetic experiments were performed under isoflurane anesthesia and all efforts were made to minimize suffering.

In vivo pharmacokinetic (PK) studies

Female CD-1 mice obtained from Biological Resource Center in Singapore were used for in vivo PK studies. Mice were acclimatized before initiation of pharmacokinetic (PK) experiments. Feed and water were given ad libitum. The compounds were formulated at a concentration of 1, 2.5 or 5 mg/mL for a dose of 10 mg/kg (Amino-824, AminoEthyl-824) or 25 mg/kg (PA-824, NI-135, NI-147, NI-136, NI-176, NI-269, NI-182, NI-145, NI-297 and NI-302) or 50 mg/kg (NI-622 and NI-644) given orally and at 1 or 2 mg/mL concentration for a dose of 5 or 10 mg/kg given intravenously. The CM-2 formulations were prepared in 10% w/v hydroxypropyl-β-cyclodextrin and 10% v/v lecithin in water as described earlier [7], [8]. The formulation was centrifuged and the supernatant was collected for intravenous administration. After oral dosing, blood and lung samples from mice were collected at various time points ranging from 0.08 hrs (but 0.02 hrs for i.v dosing) to 48 hrs. Groups of three mice were used for each time point. Blood was centrifuged at 13,000 rpm for 7 min at 4°C, plasma was harvested and stored at –20°C until analysis. Lung tissue samples were excised, dipped in PBS, gently blotted with absorbent paper, dried, weighed and stored at –20°C until further analysis. For LC/MS/MS analysis, 50 µL of plasma samples were precipitated using 400 µL of acetonitrile:methanol:acetic acid (90∶9.8∶0.2) containing 500 ng/mL of either related compound or warfarin as internal standard. After vortexing and centrifuging the mixture, the supernatant was removed and 5 µL of sample analyzed. Whole lung tissue was homogenized in 2 mL of PBS. 50 µL of lung homogenate was taken and processed as described above for plasma samples. The standard calibration curve was prepared by spiking blank plasma and lung tissue with different concentrations of the compound. In addition, quality control samples with three different concentrations were prepared in respective blank matrix and analyzed together with the unknown samples for validation purposes. Analyte quantitation was performed by LC/MS/MS using optimized conditions for each compound. Liquid chromatography was performed using an Agilent 1100 HPLC system (Santa Clara, CA), with the Agilent Zorbax XDB Phenyl (3.5 µ, 4.6×75 mm) column at an oven temperature of 35°C and 45°C, coupled with a triple quadruple mass spectrometer (Applied Biosystems, Foster City, CA). Instrument control and data acquisition were performed using Applied Biosystems software Analyst 1.4.2. The mobile phases used were A: water-acetic acid (99.8∶0.2, v/v) and B either as: acetonitrile-acetic acid (99.8∶0.2, v/v) or methanol-acetic acid (99.8∶0.2, v/v), using a gradient, with a flow rate of 1.0 mL/min, and a run time of 6 to 8 min. Under these conditions the retention times of various compounds ranged between 3.2 and 6.5 minutes. Multiple reaction monitoring (MRM) was combined with optimized mass spectrometry parameters to maximize detection specificity and sensitivity. Most of the compounds were analyzed using electrospray ionization in the positive mode. The recovery of the compound from both plasma and lung tissue were good and consistent across the concentration range studied. The lower limit of quantification for different compounds ranged between 1 and 49 ng/mL in plasma and 1 and 132 ng/g in lungs. Calibration curve was prepared freshly and analyzed with every set of study samples. Intraday variability was established with triplicate quality control samples at three concentration levels. The results were accepted if relative standard deviation was less than 15%. Mean values of compound concentrations in plasma and lungs were obtained from three animals at each time point and plotted against time to generate concentration-time profiles. PK parameters were determined using Phoenix WinNonlin, version 6.3 (Pharsight, A Certara company, USA, www.pharsight.com), by non-compartmental modeling using built-in model (200–202) for extravascular and intravenous bolus dosing. The oral bioavailability (F) was calculated as the ratio between the area under the curve (AUCinf) following oral administration and the AUCinf following intravenous administration corrected for dose (F = AUCp.o*dosei.v./AUCi.v*dosep.o).

In vivo mouse efficacy studies

In vivo mouse efficacy studies were determined after intranasal infection of Balb/c mice with 103 cfu Mtb H37Rv. Treatment was initiated 4 weeks after infection. Compounds were orally administered in CM-2 formulation for 4 weeks daily. Bacterial loads were determined at 2 and 4 weeks post treatment [30]. Statistical analysis was done by a one-way analysis of variance, followed by a multiple comparison analysis of variance by a one-way Tukey test (GraphPad Prism, version 5.02, San Diago, California USA, www.graphpad.com). Differences were considered statistically significant at the 95% level of confidence [7].

Calculation of PK-PD parameters

The MIC against Mtb was used to calculate PK-PD indices (Cmax/MIC, AUC/MIC and T). The Cmax/MIC was defined as the ratio of peak plasma concentration (Cmax) to the MIC, the AUC/MIC was defined as the ratio of the AUC0–24 to the MIC, and the time above MIC (T) was defined as 24 h period during which the total compound concentration exceeded the MIC. Cmax/MIC and AUC/MIC were calculated as ratios from PK parameter obtained from non-compartmental analysis and MIC. T were derived from Phoenix WinNonlin software by specifying MIC as therapeutic response and time above therapeutic response was obtained. Using plasma protein binding, unbound concentrations in plasma were calculated, PK parameters were derived from Phoenix WinNonlin and PK-PD indices were defined as fCmax/MIC, fAUC/MIC and % fT where ‘f’ refers to free concentration. For calculation and plotting of mean concentration-time curve, concentrations indicated as below the lower limit of quantification (LLOQ) were replaced by 0.5*LLOQ. Ignoring the values here would impact some of the PK-PD parameters. This approach has no impact on pharmacokinetic parameter calculations [31].

PK-PD analysis

PK-PD indices were estimated from the plasma and lung drug concentrations, in vitro potency and plasma protein binding. A Spearman’s rank correlation [32], [33] was run to determine the relationship between various PK parameters and mean log CFU reduction using Prism software (GraphPad Prism, version 5.02, San Diago, California USA, www.graphpad.com).

Results

In vitro potency and physico-chemical properties

In an effort to improve the solubility and potency of PA-824, diverse structural analogs of PA-824 were synthesized and their in vitro activities reported [14], [16], [18]–[20]. A few potent bicyclic 4-nitroimidazole analogs were selected and characterized in detail (Figure 1). In vitro Mtb potency and physico-chemical properties of these nitroimidazole analogs are summarized in Table 1. All the NI analogs studied showed Mtb specific growth inhibitory activity and no cytotoxicity was observed in THP1 macrophage cell lines (Table 1). F420 deficient (FbiC) mutants were resistant to all these bicyclic 4-nitroimidazoles analogs (Table 1), suggesting a mechanism of action similar to PA-824, involving F420-dependent bio-activation [24]. Modifications on the benzyl ring (NI-135, NI-147, NI-136 and NI-176), and oxazine ring (NI-269, NI-182, NI-145, NI-297, NI-302 and NI-176) showed significant improvement of in vitro potency compared to PA-824. The nitroimidazole (NI) analogs tested in this study displayed a wide range of solubility (<2 to 127 µg/mL). Amino-nitroimidazoles showed improved solubility compared to their respective benzyl ether analogs (Amino-824 vs. PA-824 and NI-269 vs. NI-145). NI-297, a biphenyl derivative of NI-182, had very poor aqueous solubility (<2 µg/mL) due to its high lipophilicity (logP 6). In general, the logP of all the other NI derivatives ranged between 2.4 and 3.8 and all showed high permeability except NI-644, which had moderate permeability in the Caco-2 assay. Overall the compounds exhibited moderate-to-high mouse plasma protein binding (80 to 98%), except for Aminoethyl-824, which showed the lowest binding (45%).
Figure 1

Chemical structures of bicyclic 4-nitroimidazole analogs used in this study.

Table 1

In vitro potency and physicochemical properties for bicyclic 4-nitroimidazole analogs.

Compound IDH37Rv MIC99(mg/L)H37Rv::F420-mutantsMIC99 (mg/L)CytotoxicityCC50 (mg/L)SolubilitypH 6.8(µg/mL)Log PPAMPA(Log Pecm/sec, %FA)CaCo-2(Papp, 10−6 cm/sec)CaCo-2(Efflux ratio)Mice PPB (%)
A-BB-AB-A/A-B
PA-8240.30>10>1013.02.75–4.2, 9927.620.40.890
NI-6220.18>10>1028.22.9–5.3, 319.512.41.397
NI-6440.09>10>1014.02.6–4.9, 832.813.74.890
Amino-8240.14>10>101272.5–4.8, 9135.421.00.680
AminoEthyl-8240.06>10>101043.6–4.6, 8518.620.71.145
NI-1350.03>10>1077.03–4.5, 9023.415.30.691.5
NI-1470.03>10nd56.02.4–5.2, 43ndnd
NI-1360.03>10>1043.02.6–4.7, 8315.412.60.882.7
NI-1760.03>10>108.03.5–3.7, 99ndnd
NI-2690.03>10>1012.03.15–3.7, 99ndnd
NI-1820.02>10>1032.03.05–3.8, 9918.313.90.891.3
NI-1450.06>10>10<2.03.85–3.85, 99ndnd
NI-2970.02>10>10<2.06–4.1, 97nd98.2
NI-3020.03>10>1026.03.2–3.6, 9945.118.70.495.7

MIC99 = Minimum Inhibitory Concentration required to reduce the bacterial growth by 99%, MIC against both H37Rv (wild type) and F420 deficient (FbiC) mutants were tested. PAMPA = Parallel Artificial Membrane Permeability Assay, CaCo-2 = Permeability using colon carcinoma cell lines, PPB = Plasma Protein Binding, Pe = effective permeability, FA = fraction absorbed, Papp = apparent permeability, A-B = Apical to Basolateral, B-A = Basolateral to Apical, nd = not determined, CC50 = Cytotoxicity against THP1 macrophage cell lines was determined as described previously [76], with Puromycin as positive control (CC50 = 1.4 mg/L).

MIC99 = Minimum Inhibitory Concentration required to reduce the bacterial growth by 99%, MIC against both H37Rv (wild type) and F420 deficient (FbiC) mutants were tested. PAMPA = Parallel Artificial Membrane Permeability Assay, CaCo-2 = Permeability using colon carcinoma cell lines, PPB = Plasma Protein Binding, Pe = effective permeability, FA = fraction absorbed, Papp = apparent permeability, A-B = Apical to Basolateral, B-A = Basolateral to Apical, nd = not determined, CC50 = Cytotoxicity against THP1 macrophage cell lines was determined as described previously [76], with Puromycin as positive control (CC50 = 1.4 mg/L).

In vivo plasma PK properties

Each compound was subjected to intravenous and oral mouse PK in CM-2 formulation. The total compound concentration in mouse plasma was measured and free plasma concentrations were calculated using in vitro plasma protein binding (Table 1 and 2). NI-147 with a hydroxyl functional group on the benzyl ring, displayed markedly inferior PK properties (Table 2, Figure S1 in File S1). The poor PK is likely due to glucuronidation of the hydroxyl group as suggested by the presence of an extra peak corresponding to +176 Da in the mass spectrometry analysis. Hence NI-147 was not included in in vivo mouse efficacy studies. The NI analogs displayed a wide range of volume of distribution (Vss = 0.7 to 4.2 L/kg) corresponding to 1.1 to 7 times total body water. NI-135 showed higher Vss (4.2 L/kg), followed by NI-182, NI-297 and NI-269 (2.6 to 3 L/kg). All other analogs showed moderate Vss similar to PA-824, except for NI-644, which showed a low volume of distribution (Vss = 0.7 L/kg). The total systemic clearance was low to moderate (4 to 44 mL/min/kg) corresponding to 5 to 49% of hepatic blood flow. NI-135 and AminoEhtyl-824 showed moderate clearance (41 and 44 mL/min/kg respectively). All other analogs showed clearance similar to PA-824 (10 to 25 mL/min/kg), except for NI-297 and NI-302, which showed very low clearance (5 mL/min/kg). The elimination half-life ranged between 0.7 h and 6.7 h for the NI analogs studied. NI-297, NI-135 and NI-302 showed long half-life (3.7 to 6.7 h). All other analogs showed half-life similar to PA-824 (1.3 to 2.8 h), except for NI-644 and AminoEthyl-824, which showed short half-life (<1 h). Generally, the NI analogs at comparable doses displayed rapid absorption (Tmax, 0.3 to 4 h), except for NI-297, which showed delayed absorption with a Tmax of 8 h, possibly due to its higher lipophilicity and lower solubility. The peak plasma concentration (Cmax) ranged between 1.2 µg/mL and 12.9 µg/mL and exposure (AUC) ranged between 4.8 µg.h/mL and 144 µg.h/mL (Table 2). NI-135 had significantly lower plasma Cmax, exposure and oral bioavailability compared to PA-824, likely due to its three times higher in vivo clearance combined with higher volume of distribution. On the contrary, NI-302 and NI-297 had higher systemic plasma exposure mostly due to decreased in vivo clearance (Table 2, Figure S1 in File S1). At comparable dose, NI-622 and NI-644 showed similar plasma Cmax, exposure and oral bioavailability compared to PA-824. All other NI analogs showed moderate plasma exposure and oral bioavailability (64 to 88%) except for NI-135, NI-297 and NI-302. Despite, NI analogs displayed varied aqueous solubility (<2 to 127 ug/mL), in CM2 formulation all the analogs showed moderate to high oral bioavailability (Table 2). Interestingly, NI-145 and NI-297 had least solubility (<2 µg/mL), in CM-2 formulation both compounds showed good oral bioavailability. The use of CM-2 (lipid coated cyclodextrin complexation) formulation is known to improve solubility and bioavailability [34]. At 25 mg/kg, the free plasma Cmax and AUC parameters ranged from 0.1–0.8 µg/mL and 0.4–5.1 µg.h/mL respectively. NI-644 showed similar free plasma concentration as PA-824, whereas all other NI analogs showed relatively lower free plasma Cmax and exposure (Table S1 in File S1).
Table 2

In vivo pharmacokinetic parameters in plasma for bicyclic 4-nitroimidazole analogs.

CompoundIDDose(mg/kg)p.o. PK parametersi.v. PK parametersa
Cmax(µg/mL)AUC24(µg.h/mL)Tmax(hr)T1/2(p.o.)F (%)Vss(L/kg)CL(mL/min/kg)T1/2(i.v)
PA-824*256.050.922.71001.612.11.6
NI-622*5014.8108.442.11001.814.71.8
NI-644*5016.289.513.61000.79.50.9
Amino-824*101.76.00.32.0742.320.52.0
AminoEthyl-824*101.02.911.7762.044.00.7
NI-135∧251.24.80.52.9514.241.04.0
NI-147∧250.040.020.10.20.30.470.80.3
NI-136∧252.010.70.52.0641.825.11.3
NI-176∧252.213.714.3861.822.21.0
NI-269∧252.816.00.33.8882.619.72.1
NI-182∧253.522.50.52.0823.015.22.8
NI-145∧251.816.224.2681.717.02.0
NI-297∧256.099.184.91002.65.06.7
NI-302∧2512.9144.144.11001.24.33.7

i.v dosing at either 10 mg/kg* or 5 mg/kg∧.

Cmax = maximum concentration reached in plasma, AUC24 = exposure between 0 to 24 h, Tmax = time to reach maximum concentration, T1/2 = half-life, F = oral bioavailability, Vss = volume of distribution at steady state, CL = total systemic clearance.

i.v dosing at either 10 mg/kg* or 5 mg/kg∧. Cmax = maximum concentration reached in plasma, AUC24 = exposure between 0 to 24 h, Tmax = time to reach maximum concentration, T1/2 = half-life, F = oral bioavailability, Vss = volume of distribution at steady state, CL = total systemic clearance.

In vivo lung PK properties

The primary and the most important site of Mtb infection in patients is lung tissue. To understand the effect of structural changes in the NI molecules on lung PK parameters, we measured total compound concentration in mouse lungs (Table 3). The NI analogs at 25 mg/kg dose showed a wide range of values for lung Cmax (4.2–17.8 µg/g), Tmax (0.08–2 h) and exposure (18.6–233 µg*h/g). All NI analogs displayed near parallel concentration-time profile in plasma and lung tissue, suggesting a rapid equilibrium between these two tissues. Interestingly, the lung-to-plasma partitioning varied from 0.5 to 4.6 for Cmax and 0.4 to 3.9 for AUCs across the series in correlation with the observed volume of distribution (Table 3). NI-135, NI-136, NI-182 and NI-297 showed lung partitioning of 2.5 to 4.6 fold, and are comparable to PA-824. NI-622, NI-644, Amino-824 and NI-302 showed lower lung partitioning (<1) compared to PA-824. Although, NI-135 and NI-136 showed higher lung to plasma ratio (3.6 to 4.6), their absolute lung concentrations were 2.5 to 7.5 fold lower than PA-824. On the contrary, NI-302 showed lower lung partitioning, but its absolute concentrations were comparable to PA-824 (Table 3).
Table 3

In vivo pharmacokinetic parameters in lungs for bicyclic 4-nitroimidazole analogs.

CompoundDose(mg/kg)Lung PK parametersLung to Plasma ratio
Cmax (µg/g)AUC24 (µg.h/g)Tmax (hr)T1/2 (p.o.)Cmax AUC24
PA-8242517.8139.90.34.83.02.7
NI-6225010.271.10.51.80.70.7
NI-644507.538.112.90.50.4
Amino-824101.44.50.51.30.80.8
AminoEthyl-824101.65.20.51.21.61.8
NI-135255.518.60.53.34.63.9
NI-147250.52.60.1---
NI-136257.239.80.52.53.63.7
NI-176254.830.30.54.32.22.2
NI-269255.427.50.32.71.91.7
NI-1822511.473.20.51.93.33.3
NI-145254.234.624.12.32.1
NI-2972516.3233.484.62.72.4
NI-302259.7100.324.10.80.7

Established mouse efficacy

Based on the in vitro potency and the in vitro and in vivo PK results, ten compounds were selected for in vivo mouse efficacy studies with 4 weeks of daily oral treatment. The mean lung CFU reductions compared to untreated controls are summarized in Table 4. The efficacy ranged from 0.5 to 1.56 log at 25 mg/kg and 0.6 to 2.3 log at 100 mg/kg compared to vehicle treated animals. At 25 mg/kg, NI-622 and NI-644 were significantly (P<0.05) less efficacious than PA-824, however other NI analogs (NI-135, NI-136, NI-182 and NI-297) showed comparable efficacy to PA-824. At 100 mg/kg, AminoEthyl-824, NI-135, NI-136 and NI-302 showed comparable efficacy to PA-824, however NI-622, NI-644, Amino-824 and NI-182 were significantly (P<0.05) less efficacious than PA-824. For PA-824, NI-135 and NI-136 a dose dependent increase in efficacy was observed, on the contrary, no dose-dependent increase in bactericidal activity was observed for NI-622, NI-644 and NI-182. However, none of the selected bicyclic 4-nitroimidazole analogs showed significantly better efficacy than PA-824 at respective 25 and 100 mg/kg doses.
Table 4

In vivo pharmacodynamics of bicyclic 4-nitroimidazole analogs studied in mice.

Mean log lung CFU ± SEM ΔMean log lung CFU reduction ± SEM
Dose (mg/kg)Vehicle control2510025100
PA-8246.07±0.124.67±0.374.12±0.131.40±0.371.95±0.13
6.24±0.024.67±0.103.93±0.091.57±0.102.31±0.09
6.66±0.145.20±0.054.19±0.111.46±0.052.47±0.11
1.48±0.09$ 2.30±0.08$
NI-6226.66±0.145.77±0.085.85±0.110.89±0.08* 0.81±0.11*
NI-6446.66±0.146.18±0.036.09±0.070.48±0.03* 0.57±0.07*
Amino-8246.22±0.08nd5.14±0.08nd1.08±0.08*
AminoEthyl-8246.22±0.08nd4.53±0.07nd1.69±0.07ns
NI-1356.24±0.024.76±0.074.36±0.071.48±0.07ns 1.88±0.07ns
NI-1366.24±0.024.93±0.064.18±0.171.31±0.06ns 2.06±0.17ns
NI-1826.07±0.124.84±0.144.81±0.181.23±0.14ns 1.26±0.18*
NI-2976.07±0.124.51±0.11nd1.56±0.11ns nd
NI-3026.62±0.11nd4.45±0.22nd2.17±0.22ns

The mean log lung CFU’s in five independent in vivo efficacy studies ranged between 6.07 and 6.66 in untreated infected mice.

nd = not determined.

Mean log lung CFU reduction compared to untreated controls. Each data represents mean value ± SEM from 5 animals.

Mean value from three independent experiments (n = 15).

*Significant difference at P<0.05 compared to respective PA-824 doses.

No significant difference at P<0.05 compared to respective PA-824 doses.

The mean log lung CFU’s in five independent in vivo efficacy studies ranged between 6.07 and 6.66 in untreated infected mice. nd = not determined. Mean log lung CFU reduction compared to untreated controls. Each data represents mean value ± SEM from 5 animals. Mean value from three independent experiments (n = 15). *Significant difference at P<0.05 compared to respective PA-824 doses. No significant difference at P<0.05 compared to respective PA-824 doses.

Correlation of PK parameters with efficacy

Further, mouse PK and efficacy data was used to understand the relationship of PK parameters with in vivo efficacy for a series of NI analogs. Both PK and efficacy data at 25 mg/kg were available for only 7 compounds. PK parameters were obtained after a single oral dose (Table 2), while efficacy studies were performed at oral daily doses of 25 and 100 mg/kg for 4 weeks (Table 4). The relationship between mean log CFU reduction with PK parameters (Cmax and AUC) was analyzed in both plasma and lungs using the Spearman’s rank correlation (Figure 2, Table S1 in File S1). The free plasma concentrations were obtained by correcting for in vitro mouse plasma protein binding. As shown in Figure 2, with the limited set of compounds, the in vivo efficacy correlated well with lung PK parameters than plasma PK parameters. The Spearman’s rank correlation coefficient for lung Cmax and AUC were 0.76 and 0.52 respectively. Across the NI analogs studied, compounds with higher lung concentration (PA-824, NI-297 and NI-182) tended to achieve higher efficacy (Δlog CFUs ranging from 1.23 to 1.56), likewise compounds with lower lung concentration (NI-644 and NI-622) displayed only marginal efficacy (Δlog CFUs ranging from 0.48 to 0.89) (Table 3 and 4). In general, lung Cmax and exposure showed positive correlation with in vivo efficacy for bicyclic 4-nitroimidazoles.
Figure 2

Correlation of PK parameters (Cmax, AUC) with in vivo efficacy in mice for bicyclic 4-nitroimidazole analogs in total plasma concentration (A), free plasma concentration (B) and total lung concentration (C).

Each data point represents Δ mean log lung CFU reduction compared to untreated controls (mean value ± SEM from 5 animals). r is the Spearman’s rank correlations coefficient.

Correlation of PK parameters (Cmax, AUC) with in vivo efficacy in mice for bicyclic 4-nitroimidazole analogs in total plasma concentration (A), free plasma concentration (B) and total lung concentration (C).

Each data point represents Δ mean log lung CFU reduction compared to untreated controls (mean value ± SEM from 5 animals). r is the Spearman’s rank correlations coefficient.

Correlation of PK-PD indices with efficacy

In vitro activity against Mtb is one of the key determinants of in vivo efficacy, hence the relationship between mean log lung CFU reduction with three primary descriptive PK-PD indices (Cmax/MIC, AUC/MIC and T) was analyzed in both plasma and lungs (Figure 3, Table S2 in File S1). As observed above, over all, the in vivo efficacy seems to have strong positive correlation with lung PK parameters than plasma. Among all the PK-PD indices, total lung T correlated the best with in vivo efficacy (r = 0.88) than lung Cmax/MIC (r = 0.63) and AUC/MIC (r = 0.63) (Figure 3C, Table S2 in File S1). For all the compounds analyzed, the total lung T ranged between 64 and 100% resulting in 0.9–1.56 log lung CFU reduction. In this analysis, NI-644 was found to be an outlier, with T of 84% resulted in only 0.48 log CFU reduction. Overall, these results suggest that in vivo efficacy of bicyclic 4-nitroimidazole analogs correlates better with the time during which the total lung concentrations are above in vitro potency.
Figure 3

Correlation of PK-PD indices (Cmax/MIC, AUC/MIC and T) with in vivo efficacy in mice for bicyclic 4-nitroimidazole analogs in total plasma concentration (A), free plasma concentration (B) and total lung concentration (C).

Each data point represents Δ Mean log lung CFU reduction compared to untreated controls (mean value ± SEM from 5 animals). r is the Spearman’s rank correlations coefficient.

Correlation of PK-PD indices (Cmax/MIC, AUC/MIC and T) with in vivo efficacy in mice for bicyclic 4-nitroimidazole analogs in total plasma concentration (A), free plasma concentration (B) and total lung concentration (C).

Each data point represents Δ Mean log lung CFU reduction compared to untreated controls (mean value ± SEM from 5 animals). r is the Spearman’s rank correlations coefficient.

Discussion

Understanding pharmacokinetic-pharmacodynamic (PK-PD) relationships in the early drug discovery process is essential to minimize the attrition rate during the pre-clinical and clinical development phases. In murine models of TB, PK-PD relationships have been established for several standard TB drugs, such as rifampicin [35], isoniazid [36], fluoroquinolones (FQ) [37] and TMC207 [38]. Based on the PK-PD findings with rifampicin, further clinical studies are still in progress to optimize the clinical dose [39]–[42]. PA-824, a bicyclic 4-nitroimidazole has demonstrated bactericidal activity in both preclinical and clinical settings [7], [8], [10]. Extensive medicinal chemistry efforts to improve aqueous solubility, metabolic stability, in vitro potency and in vivo efficacy have independently generated several series of NI analogs [5], [14]–[23]. All the bicyclic 4-nitroimidazole analogs analyzed in this study showed cofactor F420 dependent bio-activation (Table 1) suggesting the mechanism of action of these compounds similar to PA-824 [24], [43]. Comprehensive in vivo efficacy studies are generally resource/time intensive and are particularly true for TB. Thus, prioritizing potential lead compounds for in vivo efficacy studies would be useful based on PK parameters. This study is a retrospective analysis of in vivo efficacy with PK for bicyclic 4-nitroimidazole analogs to identify the PK parameters and PK-PD indices that correlate with the in vivo potency. The results of this analysis could potentially be exploited to prioritize new analogs for efficacy studies. Mtb mainly resides in lung granulomatous structures and hence it is important for a drug to be available at the site of the infection for it to be active. The volume of distribution is a primary PK parameter defined by the physico-chemical properties of the compound that indicates the extent of compound distribution in the body. Azithromycin, a macrolide antibiotic, with very high Vss (33 L/kg) is known to have higher lung concentration than serum (AUC lung/serum = 21) and it correlates well with in vivo activity against respiratory pathogens [44], [45]. Likewise, moxifloxacin displays a high volume of distribution (Vss = 2 to 5 L/kg) resulting in pronounced penetration into tissues (AUC L/P ratio of 3.3) [46], [47] possibly leading to its potent in vivo efficacy against TB [48]. Recently, moxifloxacin has been shown to penetrate and accumulate in granulomatous lesions in TB infected rabbit lungs [49]. TMC207, a diarylquinoline analog, extensively distributes to lungs (AUC L/P ratio of 22) and is efficacious against Mtb [50]. In this study, NI analogs having moderate-to-high volume of distribution (Vss = 1.6 to 4.2 L/kg) and L/P ratio of >2 showed good efficacy in a murine TB model (Δlog CFUs ranging from 1.23 to 1.56) (Table 2, 3 and 4, Figure 4). Interestingly, NI-622 and NI-644 that showed lower lung to plasma ratio displayed only a marginal efficacy (Δlog CFUs ranging from 0.48 to 0.89). Although, NI-135 and NI-136 showed higher lung to plasma ratio (3.6 to 4.6), their absolute lung concentrations were 2.5 to 7.5 fold lower than PA-824. However, both these compounds displayed 10 times better in vitro potency resulting in comparable in vivo efficacy to PA-824. Overall, the relationship between in vivo efficacy of bicyclic 4-nitroimidazoles displayed positive correlation with Vss (r = 0.45) (Figure 4). Based on these observations, the Vss and lung distribution could give an initial indication about a compound’s potential for in vivo efficacy and thus these two parameters could be used for initial prioritization of compounds during early drug discovery.
Figure 4

Correlation of volume of distribution with in vivo efficacy in mice for bicyclic 4-nitroimidazole analogs.

Each data point represents Δ Mean log lung CFU reduction compared to untreated controls (mean value ± SEM from 5 animals). r is the Spearman’s rank correlations coefficient.

Correlation of volume of distribution with in vivo efficacy in mice for bicyclic 4-nitroimidazole analogs.

Each data point represents Δ Mean log lung CFU reduction compared to untreated controls (mean value ± SEM from 5 animals). r is the Spearman’s rank correlations coefficient. A thorough dose fractionation study of PA-824 in a murine model showed that the primary PK-PD driver for in vivo efficacy is the duration during which the free concentration are above MIC (fT) in plasma [11]. In this study lung PK parameters have not been measured. Further, fT in plasma of 22%, 48% and 77% is required for it to show bacteriostatic, 1-log10 and 1.59 log10 kill respectively. In general, the PK-PD parameter driving efficacy is conserved within a given class of compounds [51], for example, the efficacy of all FQ analogs is driven by AUC/MIC [37], [52], [53], while the efficacy of β-lactams correlates with T [53]–[55]. These studies are done with thorough dose fractionation of single compound with multiple doses and dosing regimen. During lead optimization program prioritization of promising compounds that show good in vivo efficacy is important to reduce the overall turnaround time. In this retrospective analysis with 7 different bicyclic 4-nitroimidazole analogs, we attempted to correlate in vivo efficacy at 25 mg/kg with PK parameters. On the contrary to what has been observed by Ahmad et al., in this study, with seven bicyclic 4-nitroimidazole analogs having varied lung distribution, in vivo efficacy showed weak correlation with free T in plasma. However, the total T in lungs showed positive correlation with in vivo efficacy (r = 0.88) likely due to their preferential distribution into lungs for some analogs. For all the compounds analyzed, the total lung T ranged between 64–100% resulting in 0.9–1.56 log lung CFU reduction, hence efficacy studies at lower doses (resulting in T less than 65%) might be necessary to see a better correlation. Overall, in this study a diverse set of bicyclic 4-nitroimidazoles with Vss ranging from 0.7 L/kg to 4.2 L/kg, lung to plasma ratio ranging from 0.5 to 4.6 showed positive correlation with lung T than with any other parameters. The results presented in this study must be interpreted with a couple of limitations in mind. First, single-dose PK parameters determined in healthy mice were assumed to be similar to multiple-dose PK parameters in infected animals and were correlated with efficacy data. This assumption is supported by published preclinical data that has shown the absence of plasma accumulation of PA-824 in mice dosed for 2 months [56]. Further, in clinical studies with PA-824, the PK parameters from a single dose phase I study were similar to a multiple dose phase II study in patients [9], [10], [12], [13]. Another limitation of this study is that total concentrations in lungs rather than the free lung concentrations were used for the PK-PD analysis. It is well accepted that for a given compound unbound drug concentrations in plasma are equivalent to unbound tissue concentrations when active transport is not involved in the drug distribution [57], [58]. Further, it is the unbound concentration of a compound at its target site driving the pharmacological effect [58]–[63]. Nevertheless, whole-tissue concentrations can be of some value in early drug discovery providing a first assessment of partition into the lungs [61]. Techniques like microdialysis in lungs can be applied to assess unbound tissue concentration [64], [65]. In TB patients, Mtb mainly resides in diverse and heterogeneous lesions in lungs. In general, interpretation of PD activity of anti-TB compounds is complicated by differential lung pathophysiology. PK in intrapulmonary compartments like the epithelial lining fluid and alveolar macrophages have also been studied in humans for standard TB drugs like rifampicin [66], isoniazid [67], ethambutol [68], pyrazinamide [69], rifapentine [70], moxifloxacin [71], ofloxacin [72] and linezolid [73]. The concentration in these sites could be the key factor governing the efficacy of anti-TB drugs. However, measurement of compound concentration in lungs by microdialysis, epithelial lining fluid and alveolar macrophages have limitations in sampling, methodology and interpretation of results [61], [74]; and such studies have not been explored in preclinical settings for TB. The total lung concentration may not be equal to the concentration in Mtb lesions, thus warranting lesion PK analysis to improve the predictive power for efficacy. Recently PK in lung lesions of mycobacterium-infected rabbits has been investigated for isoniazid, rifampicin, pyrazinamide and moxifloxacin [49], [75]. Although lesion PK can offer better insights in understanding PK-PD relationships, it is not easily applicable to early drug discovery especially with mouse efficacy model as it doesn’t display spectrum of lesions observed in TB patients or in higher animal models. In addition, similar studies with bicyclic 4-nitroimidazoles may be challenging as they undergo enzymatic transformation in Mtb to multiple stable and unstable metabolites [43]. Our findings show that the efficacy of all bicyclic 4-nitroimidazole analogs is most likely driven by PK parameters in lungs. A simple efficacy surrogate would be useful during the lead optimization to prioritize candidates for lengthy efficacy studies. For this class, efficacy correlated better with concentration in lungs rather than in plasma, consistent with Vss and differential lung: plasma distributions. The results of this analysis potentially be exploited to prioritize new analogs for efficacy studies based on in vitro potency, volume of distribution and lung concentration. Figure S1: Plasma concentration time profile for representative bicyclic 4-nitroimidazole analogs following a single 25 mg/kg oral dose in mice. Table S1: Correlation of PK parameters with in vivo efficacy in mice for bicyclic 4-nitroimidazole analogs. Table S2: Correlation of PK-PD indices with in vivo efficacy in mice for bicyclic 4-nitroimidazole analogs. (DOCX) Click here for additional data file.
  71 in total

1.  Concentrations of moxifloxacin in serum and pulmonary compartments following a single 400 mg oral dose in patients undergoing fibre-optic bronchoscopy.

Authors:  A Soman; D Honeybourne; J Andrews; G Jevons; R Wise
Journal:  J Antimicrob Chemother       Date:  1999-12       Impact factor: 5.790

2.  Structure-activity relationships for amide-, carbamate-, and urea-linked analogues of the tuberculosis drug (6S)-2-nitro-6-{[4-(trifluoromethoxy)benzyl]oxy}-6,7-dihydro-5H-imidazo[2,1-b][1,3]oxazine (PA-824).

Authors:  Adrian Blaser; Brian D Palmer; Hamish S Sutherland; Iveta Kmentova; Scott G Franzblau; Baojie Wan; Yuehong Wang; Zhenkun Ma; Andrew M Thompson; William A Denny
Journal:  J Med Chem       Date:  2011-12-29       Impact factor: 7.446

Review 3.  Microdialysis for in vivo pharmacokinetic/pharmacodynamic characterization of anti-infective drugs.

Authors:  Martin Brunner; Hartmut Derendorf; Markus Müller
Journal:  Curr Opin Pharmacol       Date:  2005-10       Impact factor: 5.547

4.  Structure-activity relationships of antitubercular nitroimidazoles. 3. Exploration of the linker and lipophilic tail of ((s)-2-nitro-6,7-dihydro-5H-imidazo[2,1-b][1,3]oxazin-6-yl)-(4-trifluoromethoxybenzyl)amine (6-amino PA-824).

Authors:  Joseph Cherian; Inhee Choi; Amit Nayyar; Ujjini H Manjunatha; Tathagata Mukherjee; Yong Sok Lee; Helena I Boshoff; Ramandeep Singh; Young Hwan Ha; Michael Goodwin; Suresh B Lakshminarayana; Pornwaratt Niyomrattanakit; Jan Jiricek; Sindhu Ravindran; Thomas Dick; Thomas H Keller; Veronique Dartois; Clifton E Barry
Journal:  J Med Chem       Date:  2011-07-26       Impact factor: 7.446

Review 5.  Challenges and opportunities in developing novel drugs for TB.

Authors:  Takushi Kaneko; Christopher Cooper; Khisimuzi Mdluli
Journal:  Future Med Chem       Date:  2011-09       Impact factor: 3.808

6.  Pharmacokinetic evaluation of the penetration of antituberculosis agents in rabbit pulmonary lesions.

Authors:  Maria C Kjellsson; Laura E Via; Anne Goh; Danielle Weiner; Kang Min Low; Steven Kern; Goonaseelan Pillai; Clifton E Barry; Véronique Dartois
Journal:  Antimicrob Agents Chemother       Date:  2011-10-10       Impact factor: 5.191

7.  Overton's rule helps to estimate the penetration of anti-infectives into patients' cerebrospinal fluid.

Authors:  Marija Djukic; Martin Munz; Fritz Sörgel; Ulrike Holzgrabe; Helmut Eiffert; Roland Nau
Journal:  Antimicrob Agents Chemother       Date:  2011-11-21       Impact factor: 5.191

8.  Pharmacokinetics and pharmacodynamics of TMC207 and its N-desmethyl metabolite in a murine model of tuberculosis.

Authors:  Marie-Claude Rouan; Nacer Lounis; Tom Gevers; Lieve Dillen; Ron Gilissen; Araz Raoof; Koen Andries
Journal:  Antimicrob Agents Chemother       Date:  2011-12-12       Impact factor: 5.191

Review 9.  Why Do We Use 600 mg of Rifampicin in Tuberculosis Treatment?

Authors:  Jakko van Ingen; Rob E Aarnoutse; Peter R Donald; Andreas H Diacon; Rodney Dawson; Georgette Plemper van Balen; Stephen H Gillespie; Martin J Boeree
Journal:  Clin Infect Dis       Date:  2011-05       Impact factor: 9.079

10.  Synthesis and structure-activity relationships of varied ether linker analogues of the antitubercular drug (6S)-2-nitro-6-{[4-(trifluoromethoxy)benzyl]oxy}-6,7-dihydro-5h-imidazo[2,1-b][1,3]oxazine (PA-824).

Authors:  Andrew M Thompson; Hamish S Sutherland; Brian D Palmer; Iveta Kmentova; Adrian Blaser; Scott G Franzblau; Baojie Wan; Yuehong Wang; Zhenkun Ma; William A Denny
Journal:  J Med Chem       Date:  2011-09-01       Impact factor: 7.446

View more
  12 in total

1.  Synthesis and evaluation of pretomanid (PA-824) oxazolidinone hybrids.

Authors:  David F Bruhn; Michael S Scherman; Aman P Singh; Lei Yang; Jiuyu Liu; Anne J Lenaerts; Richard E Lee
Journal:  Bioorg Med Chem Lett       Date:  2015-12-07       Impact factor: 2.823

2.  Synthesis and Biological Evaluation of Polar Functionalities Containing Nitrodihydroimidazooxazoles as Anti-TB Agents.

Authors:  Kushalava Reddy Yempalla; Gurunadham Munagala; Samsher Singh; Gurleen Kour; Shweta Sharma; Reena Chib; Sunil Kumar; Priya Wazir; G D Singh; Sushil Raina; Sonali S Bharate; Inshad Ali Khan; Ram A Vishwakarma; Parvinder Pal Singh
Journal:  ACS Med Chem Lett       Date:  2015-09-11       Impact factor: 4.345

3.  Modeling and Simulation of Pretomanid Pharmacodynamics in Pulmonary Tuberculosis Patients.

Authors:  Michael A Lyons
Journal:  Antimicrob Agents Chemother       Date:  2019-09-30       Impact factor: 5.191

4.  Mutations in genes for the F420 biosynthetic pathway and a nitroreductase enzyme are the primary resistance determinants in spontaneous in vitro-selected PA-824-resistant mutants of Mycobacterium tuberculosis.

Authors:  Hana L Haver; Adeline Chua; Pramila Ghode; Suresh B Lakshminarayana; Amit Singhal; Barun Mathema; René Wintjens; Pablo Bifani
Journal:  Antimicrob Agents Chemother       Date:  2015-06-22       Impact factor: 5.191

5.  Modeling and Simulation of Pretomanid Pharmacokinetics in Pulmonary Tuberculosis Patients.

Authors:  Michael A Lyons
Journal:  Antimicrob Agents Chemother       Date:  2018-06-26       Impact factor: 5.191

6.  Identification of a new series of benzothiazinone derivatives with excellent antitubercular activity and improved pharmacokinetic profiles.

Authors:  Lu Xiong; Chao Gao; Yao-Jie Shi; Xin Tao; Juan Rong; Kun-Lin Liu; Cui-Ting Peng; Ning-Yu Wang; Qian Lei; Yi-Wen Zhang; Luo-Ting Yu; Yu-Quan Wei
Journal:  RSC Adv       Date:  2018-03-21       Impact factor: 4.036

Review 7.  New Drugs for the Treatment of Tuberculosis.

Authors:  Elisa H Ignatius; Kelly E Dooley
Journal:  Clin Chest Med       Date:  2019-12       Impact factor: 4.967

Review 8.  Nitroaromatic Antibiotics as Nitrogen Oxide Sources.

Authors:  Allison M Rice; Yueming Long; S Bruce King
Journal:  Biomolecules       Date:  2021-02-12

9.  Delamanid or pretomanid? A Solomonic judgement!

Authors:  Saskia E Mudde; Anna M Upton; Anne Lenaerts; Hannelore I Bax; Jurriaan E M De Steenwinkel
Journal:  J Antimicrob Chemother       Date:  2022-03-31       Impact factor: 5.790

10.  Assessment of a pretomanid analogue library for African trypanosomiasis: Hit-to-lead studies on 6-substituted 2-nitro-6,7-dihydro-5H-imidazo[2,1-b][1,3]thiazine 8-oxides.

Authors:  Andrew M Thompson; Andrew J Marshall; Louis Maes; Nigel Yarlett; Cyrus J Bacchi; Eric Gaukel; Stephen A Wring; Delphine Launay; Stephanie Braillard; Eric Chatelain; Charles E Mowbray; William A Denny
Journal:  Bioorg Med Chem Lett       Date:  2017-10-27       Impact factor: 2.823

View more

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