[Shen Shiwan : A translator for medicine in the Period of the Republic of China].

Shen Shiwan was a translator and physician for both Chinese medicine and western medicine in the period of the Republic of China. This paper examined the life and the main academic contributions of Shen Shiwan. It was found that Shen's main contribution were translating medical works, founding journals and opening the door of Chinese medicine to the world. Additionally, he introduced western medicine, such as pathology, fertility and nutrition to Chinese medical professionals. He also introduced the medical schools of Japanese traditional medicine (Chinese medicine in Japan, Han Yi) into China. Shen's contribution in medicine played an important role for the medical professionals in China in understanding the development of western medicine and Chinese medicine in Japan.

Mitochondrial Stress in Metabolic Inflammation: Modest Benefits and Full Losses.

Energy intake and metabolic balance are the pillars of health preservation. Overnutrition causes nonspecific, persistently low inflammatory state known as metabolic inflammation. This condition contributes to the pathophysiology of various metabolic disorders, such as atherosclerosis, obesity, diabetes mellitus, and metabolic syndrome. The mitochondria maintain the balance of energy metabolism. Excessive energy stress can lead to mitochondrial dysfunction, which promotes metabolic inflammation. The inflammatory environment further impairs mitochondrial function. Accordingly, excellent organism design keeps the body metabolically healthy in the context of mitochondrial dysfunction, and moderate mitochondrial stress can have a beneficial effect. This review summarises the research progress on the multifaceted characterisation of mitochondrial dysfunction and its role in metabolic inflammation.

Insights Into the Role of Mitochondria in Vascular Calcification.

Vascular calcification (VC) is a growing burden in aging societies worldwide, and with a significant increase in all-cause mortality and atherosclerotic plaque rupture, it is frequently found in patients with aging, diabetes, atherosclerosis, or chronic kidney disease. However, the mechanism of VC is still not yet fully understood, and there are still no effective therapies for VC. Regarding energy metabolism factories, mitochondria play a crucial role in maintaining vascular physiology. Discoveries in past decades signifying the role of mitochondrial homeostasis in normal physiology and pathological conditions led to tremendous advances in the field of VC. Therapies targeting basic mitochondrial processes, such as energy metabolism, damage in mitochondrial DNA, or free-radical generation, hold great promise. The remarkably unexplored field of the mitochondrial process has the potential to shed light on several VC-related diseases. This review focuses on current knowledge of mitochondrial dysfunction, dynamics anomalies, oxidative stress, and how it may relate to VC onset and progression and discusses the main challenges and prerequisites for their therapeutic applications.

Magic and mystery of microRNA-32.

MicroRNAs (miRNAs) are a group of endogenous, small (∼22 nts in length) noncoding RNA molecules that function specifically by base pairing with the mRNA of genes and regulate gene expression at the post-transcriptional level. Alterations in miR-32 expression have been found in numerous diseases and shown to play a vital role in cell proliferation, apoptosis, oncogenesis, invasion, metastasis and drug resistance. MiR-32 has been documented as an oncomiR in the majority of related studies but has been also verified as a tumour suppressor miRNA in conflicting reports. Moreover, it has a crucial role in metabolic and cardiovascular disorders. This review provides an in-depth look into the most recent finding regarding miR-32, which is involved in the expression, regulation and functions in different diseases, especially tumours. Additionally, this review outlines novel findings suggesting that miR-32 may be useful as a noninvasive biomarker and as a targeted therapeutic in several diseases.

Highly expressed lncRNA FOXD3-AS1 promotes non-small cell lung cancer progression via regulating miR-127-3p/mediator complex subunit 28 axis.

OBJECTIVE: The present study aimed to determine the expression of long non-coding RNA (lncRNA) FOXD3 antisense RNA 1 (FOXD3-AS1) in lung cancer tissues and to explore its underlying mechanisms in mediating non-small cell lung cancer (NSCLC) progression. MATERIALS AND METHODS: Gene expression levels were determined by quantitative real-time PCR; lung cancer cell proliferation and invasion were determined by in vitro functional assays; protein levels were determined by Western blot assay; xenograft nude mice model was used to evaluate the in vivo tumor growth of lung cancer cells; Luciferase reporter assay determined the interactions among FOXD3-AS1, miR-127-3p, and mediator complex subunit 28 (MED28). RESULTS: Data mining and analysis of the clinical sample showed that FOXD3-AS1 expression was significantly up-regulated in lung cancer tissues. In vitro functional assays demonstrated that FOXD3-AS1 overexpression promoted NSCLC cell proliferation and invasion, while FOXD3-AS1 knockdown exerted tumor-suppressive effects on NSCLC cells. Moreover, FOXD3-AS1 interacted with miR-127-3p by acting as a competing endogenous RNA to suppress miR-127-3p expression, while miR-127-3p repressed MED28 expression by targeting MED28 3' untranslated region in NSCLC cells. Mechanistically, the oncogenic effects of FOXD3-AS1 overexpression were significantly attenuated by miR-127-3p overexpression and MED28 knockdown in NSCLC cells. In the xenograft mice model, FOXD3-AS1 knockdown suppressed in vivo tumor growth of A549 cells, and also up-regulated miR-127-3p expression and repressed MED28 expression in the xenograft tumors. In the clinical aspect, the downregulation of miR-127-3p and up-regulation of MED28 were respectively detected in lung cancer tissues. CONCLUSIONS: Our findings provided new evidence that the FOXD3-AS1 regulated NSCLC progression via targeting the miR-127-3p/MED28 axis.

[Clinical efficacy of anterior decompression and stability reconstruction in patients with cervical hyperextension injury].

Objective: To evaluate the clinical efficacy of anterior decompression and stability reconstruction in patients with cervical hyperextension injury. Methods: Postoperative data from 60 patients with cervical hyperextension injury between April 2009 and December 2016 were analyzed retrospectively; the patients included 50 males and 10 females, aged 21-87 years [average, (57±13) years]. All patients had various degrees of spinal cord injury, and were treated with anterior cervical decompression, fusion, and internal fixation.The preoperative and postoperative neurological function were compared to evaluate the clinical efficacy of the treatment.The t test was applied when preoperative and postoperative data were compared. Results: Of the 60 patients, 5 underwent anterior cervical discectomy and fusion, 26 underwent anterior cervical corpectomy and fusion, and 29 received treatment with the hybrid technique.The average follow-up was (5.1±2.1) years (range, 1.6-9.1 years). The American Spinal Injury Association (ASIA) scores, abbreviated injury scale (AIS grades), and Japanese Orthopedic Association Scores (JOA scores) at the 1 week post operation and final follow-up were significantly better than those obtained preoperatively (all P<0.01). The JOA scores before operation and at the final follow-up was 10(7, 11) and 16(14, 17), respectively.Based on the recovery rate of JOA scores, the rate of cure was 28.3% (17 cases), the rate of significant efficiency was 60.0% (36 cases), the rate of efficiency was 8.3% (5 cases), the rate of inefficiency was 3.3% (2 cases). Among the patients showing cure, 5 were satisfied with the life function, however, the extent of their injury was still of Grade D owing to the incomplete recovery of muscle force.Further, when the 60 patients were divided into A and B groups according to whether they had congenital spinal stenosis or not, respectively, and no significant difference was found in ASIA scores, AIS grades and JOA scores between the groups (all P<0.01). Conclusion: Anterior cervical decompression, fusion and internal fixation is safe and effective for cervical hyperextension injury without continuous-type ossification of the posterior longitudinal ligament.

Disrupting G6PD-mediated Redox homeostasis enhances chemosensitivity in colorectal cancer.

Glucose-6-phosphate dehydrogenase (G6PD) is a key enzyme that generates NADPH to maintain reduced glutathione (GSH), which scavenges reactive oxygen species (ROS) to protect cancer cell from oxidative damage. In this study, we mainly investigate the potential roles of G6PD in colorectal cancer (CRC) development and chemoresistance. We discover that G6PD is overexpressed in CRC cells and patient specimens. High expression of G6PD predicts poor prognosis and correlated with poor outcome of oxaliplatin-based first-line chemotherapy in patients with CRC. Suppressing G6PD decreases NADPH production, lowers GSH levels, impairs the ability to scavenge ROS levels, and enhances oxaliplatin-induced apoptosis in CRC via ROS-mediated damage in vitro. In vivo experiments further shows that silencing G6PD with lentivirus or non-viral gene delivery vector enhances oxaliplatin anti-tumor effects in cell based xenografts and PDX models. In summary, our finding indicated that disrupting G6PD-mediated NADPH homeostasis enhances oxaliplatin-induced apoptosis in CRC through redox modulation. Thus, this study indicates that G6PD is a potential prognostic biomarker and a promising target for CRC therapy.

[Clinical outcomes of single-level lumbar spondylolisthesis by minimally invasive transforaminal lumbar interbody fusion with bilateral tubular channels].

Objective: To evaluate the clinical effectiveness of minimally invasive transforaminal lumbar interbody fusion (MIS-TLIF) for single-level lumbar spondylolisthesis treatment with bilateral Spotlight tubular channels. Methods: A total of 21 patients with lumbar spondylolisthesis whom underwent MIS-TLIF via bilateral Spotlight tubular channels were retrospectively analyzed from October 2014 to November 2015. The 21 patients included 11 males and 10 females ranged from 35 to 82 years (average aged 60.7 years). In term of spondylolisthesis category, there were 18 cases of degenerative spondylolisthesis and 3 cases of isthmic spondylolisthesis. With respect to spondylolisthesis degree, 17 cases were grade Ⅰ° and 4 cases were grade Ⅱ°. Besides, 17 cases at L(4-5) and 4 cases at L(5)-S(1)were categorized by spondylolisthesis levels. Operation duration, blood loss, postoperative drainage and intraoperative exposure time were recorded, functional improvement was defined as an improvement in the Oswestry Disability Index (ODI), Visual Analog Scale (VAS) was also employed at pre and post-operation (3 months and the last follow-up), to evaluate low back and leg pain. Furthermore, to evaluate the recovery of the intervertebral foramen and of lumbar sagittal curvature, average height of intervertebral space, Cobb angles of lumbar vertebrae and operative segments, spondylolisthesis index were measured. At the last follow-up, intervertebral fusion was assessed using Siepe evaluation criteria and the clinical outcome was assessed using the MacNab scale. Radiographic and functional outcomes were compared pre- and post-operation using the paired T test to determine the effectiveness of MIS-TLIF. Statistical significance was defined as P<0.05. Results: All patients underwent a successful MIS-TLIF surgery. The operation time (235.2±30.2) mins, intraoperative blood loss (238.1±130.3) ml, postoperative drainage (95.7±57.1) ml and intraoperative radiation exposure (47.1±8.8) were recorded. Different significance between 3 months post-operative follow-up and pre-operation was exhibited (P<0.01) in respects of lumbar VAS (t=11.1, P<0.01) and leg VAS (t=17.8, P<0.01). Moreover, final follow-up compared with pre-operation, and final follow-up compared with 3 months post-operative follow-up, VAS scores were also statistical difference (P<0.01). At the final follow-up, there were significant differences compared with pre-operation in ODI scores (t=30.1, P<0.01). Comparison between 3 months post-operative follow-up and pre-operation, statistical distinctions were demonstrated (P<0.05) in terms of mean height of intervertebral space (t=-10.9, P<0.01), the Cobb angles of lumbar vertebrae (t=-2.4, P<0.05), operative segments Cobb angles (t=-5.2, P<0.01) and Lumbar spondylolisthesis incidence (t=17.1, P<0.01). In addition, there was statistical difference between final follow-up and pre-operation (P<0.05) as well. For instance, mean height of intervertebral space (t=-10.5, P<0.01), the Cobb angles of lumbar vertebrae (t=-2.7, P<0.05), operative segments Cobb angles (t=-4.2, P<0.01) and Lumbar spondylolisthesis incidence (t=18.6, P<0.01) were involved. All spondylolisthesis vertebrae were restored completely. Lastly, at the last follow-up, 12 cases of grade 1 and 7 cases of grade 2 fusion were present as determined by the Siepe evaluation criteria. McNab scale assessment classified 17 patients having excellent clinical outcome, 3 patients in good and 1 patient having a better clinical outcome. Conclusion: MIS-TLIF with bilateral Spotlight tubular channels is a safe and effective approach for single segment lumbar spondylolisthesis.

Comparative pharmacokinetics of diaveridine in pigs and chickens following single intravenous and oral administration.

Comparative pharmacokinetic profiles of diaveridine following single intravenous and oral dose of 10 mg/kg body weight in healthy pigs and chickens were investigated, respectively. Concentrations of diaveridine in plasma samples were determined using a validated high-performance liquid chromatography-ultraviolet (HPLC-UV) method. The concentration-time data were subjected to noncompartmental kinetic analysis by WinNonlin program. The corresponding pharmacokinetic parameters in pigs or chickens after single intravenous administration were as follows, respectively: t1/2ß (elimination half-life) 0.74 ± 0.28 and 3.44 ± 1.07 h; Vd (apparent volume of distribution) 2.70 ± 0.99 and 3.86 ± 0.92 L/kg; ClB (body clearance) 2.59 ± 0.62 and 0.80 ± 0.14 L/h/kg; and AUC0-∞ (area under the blood concentration vs. time curve) 4.11 ± 1.13 and 12.87 ± 2.60 µg∙h/mL. The corresponding pharmacokinetic parameters in pigs or chickens after oral administration were as follows, respectively: t1/2ß 1.78 ± 0.41 and 2.91 ± 0.57 h; Cmax (maximum concentration) 0.43 ± 0.24 and 1.45 ± 0.57 µg/mL; Tmax (time to reach Cmax ) 1.04 ± 0.67 and 3.25 ± 0.71 h; and AUC0-∞ 1.33 ± 0.55 and 9.28 ± 2.69 µg∙h/mL. The oral bioavailability (F) of diaveridine in pigs or chickens was determined to be 34.6% and 72.2%, respectively. There were significant differences between the pharmacokinetics profiles in these two species.

[Survey of the first case of Zika virus disease in the mainland of China].

OBJECTIVE: To analyze epidemiological characteristics of the first laboratory confirmed case of Zika virus disease in the mainland China, and provide evidence for the prevention and control of the spread of Zika virus disease. METHODS: Epidemiological survey was conducted for the first suspect case of Zika virus disease in China, and medical observation was conducted for the close contacts. The nucleic acid detection of Zika virus was conducted with Real-time RT-PCR by using blood and urine samples collected from the cases. RESULTS: The first case of Zika virus disease was confirmed. The case was cured and discharged on 18(th) day after the onset. The blood sample on 10(th) day after the onset and urine samples on 11-13(th) days after the onset were all positive for Zika virus. The case had the history of mosquito bites in the endemic country before the onset. The close contacts showed no Zika virus disease-like symptoms during the medical observation period. CONCLUSION: The case reported by Jiangxi province was the first imported case of Zika virus disease in the mainland of China. The infection was related with the mosquito bites in Venezuela where Zika virus disease in endemic.

Estimating danofloxacin withdrawal time in broiler chickens based on physiologically based pharmacokinetics modeling.

In this study, a physiologically based pharmacokinetics (PBPK) model was firstly developed for danofloxacin in healthy broiler chickens after a single oral administration at 5 mg/kg bw. Then, the model extrapolation from healthy chickens to those infected with Pasteurella multocidaones was performed. The healthy model was validated through a comparison of predicted and previously published concentrations, which indicated that the healthy PBPK model had good predictive ability in plasma, lung, muscle, liver, and kidney, especially at the later sampling time points. Multiple dosing of administration was incorporated into the healthy and infected models. In addition, a Monte Carlo simulation (MCS) included 1000 iterations was further incorporated into both models to predict the withdrawal times of danofloxacin in healthy and infected chickens, which were estimated to be 3 and 2 days, respectively.

Estimating marbofloxacin withdrawal time in broiler chickens using a population physiologically based pharmacokinetics model.

Residue depletion of marbofloxacin in broiler chicken after oral administration at 5 mg/kg/day for three consecutive days was studied in this study. The areas under the concentration-time curve from 0 h to ∞ (AUC0-∞ s) of marbofloxacin in tissues and plasma were used to calculate tissue/plasma partition coefficients (PX s). Based on PX s and the other parameters derived from published studies, a flow-limited physiologically based pharmacokinetics (PBPK) model was developed to predict marbofloxacin concentrations, which were then compared with those derived from the residue depletion study so as to validate this model. Considering individual difference in drug disposition, a Monte Carlo simulation included 1000 iterations was further incorporated into the validated model to generate a population PBPK model and to estimate the marbofloxacin residue withdrawal times in edible tissues. The withdrawal periods were compared to those derived from linear regression analysis. The PBPK model presented here successfully predicted the measured concentrations in all tissues. The withdrawal times in all edible tissues derived from the population PBPK model were longer than those from linear regression analysis, and based on the residues in kidney, a withdrawal time of 4 days was estimated for marbofloxacin after oral administration at 5 mg/kg/day for three consecutive days. It was shown that population PBPK model could be used to accurately predict marbofloxacin residue withdrawal time in edible tissues in broiler chickens.

Pharmacokinetics of doxycycline in tilapia (Oreochromis aureus × Oreochromis niloticus) after intravenous and oral administration.

The pharmacokinetics of doxycycline was studied in plasma after a single dose (20 mg/kg) of intravenous or oral administration to tilapia (Oreochromis aureus × Oreochromis niloticus) reared in fresh water at 24 °C. Plasma samples were collected from six fish per sampling point. Doxycycline concentrations were determined by high-performance liquid chromatography with a 0.005 µg/mL limit of detection, then were subjected to noncompartmental analysis. Following oral administration, the double-peak phenomenon was observed, and the first (Cmax1 ) and second (Cmax2) peaks were 1.99 ± 0.43 µg/mL at 2.0 h and 2.27 ± 0.38 µg/mL at 24.0 h, respectively. After the intravenous injection, a Cmax2 (12.12 ± 1.97 µg/mL) was also observed, and initial concentration of 45.76 µg/mL, apparent elimination rate constant (λz) of 0.018 per h, apparent elimination half-life (t1/2λz) of 39.0 h, systemic total body clearance (Cl) of 41.28 mL/h/kg, volume of distribution (Vz) of 2323.21 mL/kg, and volume of distribution at steady-state (Vss) of 1356.69 mL/kg were determined, respectively. While after oral administration, the λz, t1/2λz, and bioavailability of doxycycline were 0.009 per h, 77.2 h, and 23.41%, respectively. It was shown that doxycycline was relatively slowly and incompletely absorbed, extensively distributed, and slowly eliminated in tilapia, in addition, doxycycline might undergo enterohepatic recycling in tilapia.

Metabolic regulation of cancer cell side population by glucose through activation of the Akt pathway.

Side population (SP) cells within tumors are a small fraction of cancer cells with stem-like properties that can be identified by flow cytometry analysis based on their high ability to export certain compounds such as Hoechst 33342 and chemotherapeutic agents. The existence of stem-like SP cells in tumors is considered as a key factor contributing to drug resistance, and presents a major challenge in cancer treatment. Although it has been recognized for some time that tumor tissue niches may significantly affect cancer stem cells (CSCs), the role of key nutrients such as glucose in the microenvironment in affecting stem-like cancer cells and their metabolism largely remains elusive. Here we report that SP cells isolated from human cancer cells exhibit higher glycolytic activity compared to non-SP cells. Glucose in the culture environment exerts a profound effect on SP cells as evidenced by its ability to induce a significant increase in the percentage of SP cells in the overall cancer cell population, and glucose starvation causes a rapid depletion of SP cells. Mechanistically, glucose upregulates the SP fraction through ATP-mediated suppression of AMPK and activation of the Akt pathway, leading to elevated expression of the ATP-dependent efflux pump ABCG2. Importantly, inhibition of glycolysis by 3-BrOP significantly reduces SP cells in vitro and impairs their ability to form tumors in vivo. Our data suggest that glucose is an essential regulator of SP cells mediated by the Akt pathway, and targeting glycolysis may eliminate the drug-resistant SP cells with potentially significant benefits in cancer treatment.

The correlation between vertebral wedge-shaped changes in X-ray imaging at supine and standing positions and the efficacy of operative treatment of thoracolumbar spinal fracture in the elderly.

STUDY DESIGN: By analyzing a large number of surgical patients, we identified the roles of wedge-shaped changes in related surgeries. OBJECTIVES: To illustrate the relevance of vertebral wedge-shaped changes in X-ray imaging at supine and standing positions in patients with percutaneous kyphoplasty as well as the postoperative effect. SETTING: All patient data were collected from a hospital in China. METHODS: Between June 2006 and May 2010, 77 surgical patients (9 men and 68 women) with wedge-shaped compression fractures were retrospectively analyzed. Patients were divided into group A (ΔWR2.5%) and group B (ΔWR<2.5%) according to the dynamic changes in the percentage of vertebral body wedge-shaped variable ratio (WR) at supine and standing positions. The intensity of back pain in different positions pre- and postoperatively was evaluated with a visual analog pain scale (VAS). RESULTS: The WRs in both standing and supine positions were significantly reduced by kyphoplasty in both groups A and B. In agreement with the improvement in WRs, the VAS was significantly decreased in three positions for patients in group A and in turning over and standing position for patients in group B. With respect to ΔWR changes, group B revealed significantly lower values compared with group A preoperatively (P<0.001), but there was no significant difference between groups A and B postoperatively and at 1-month follow-up (P=0.179 and P=0.558, respectively). CONCLUSIONS: Improvement in symptoms after kyphoplasty is better in patients with wedge-shaped changes in supine and standing positions, and the efficacy of height restoration of the spine would be better in unstable vertebrae by balloon dilatation.

Estimating tulathromycin withdrawal time in pigs using a physiologically based pharmacokinetics model.

A physiologically based pharmacokinetics model was developed to predict tulathromycin concentrations in edible swine tissues. Physiological parameters included volumes of and plasma flows through different tissues which were obtained from the literatures. The tissue/plasma partition coefficient was calculated according to the area method, and the model was validated through a comparison of predicted and observed concentrations. Withdrawal times in different tissues were predicted. The physiologically based pharmacokinetics model presented here provided accurate predictions of the observed concentrations in all tissues. The results showed that the injection site had the longest withdrawal time (21 days), followed by skin together with fat (19 days) and then kidney (10 days), lung (6 days), liver (4 days) and muscle (1 day). A withdrawal time of 21 days was finally predicted for tulathromycin in swine after a single intramuscular injection at 2.5 mg/kg body weight.

A physiologically based pharmacokinetics model for florfenicol in crucian carp and oral-to-intramuscular extrapolation.

In this study, an oral physiologically based pharmacokinetics (PBPK) model was developed for florfenicol in crucian carp (Carassius auratus). Subsequently, oral-to-intramuscular extrapolation was performed and the two models were used to predict florfenicol concentrations in the edible tissues of crucian carp. The oral model gave good predictions in most tissues, except for kidney and liver in which the florfenicol concentrations were underestimated at the later time points. In contrast, using the intramuscular model, the concentrations in the kidney were overestimated at the later time points. Both models had the best predictive ability in the main edible tissue, the muscle. The oral model also accurately predicted the florfenicol concentrations in the muscle after multiple doses. The present study demonstrated the feasibility of predicting florfenicol concentrations in the edible tissues of crucian carp using a route-to-route extrapolation method.

Potential role of AQP1 and VEGF in the development of malignant pleural effusion in mice.

Recent studies have revealed the important role of aquaporin-1 (AQP1) in tumor cell migration and angiogenesis. However, the function of AQP1 in malignant pleural effusion has not been well characterized. We established a mouse model to examine the role of AQP1 and vascular endothelial growth factor (VEGF) in the development of malignant pleural effusion. We showed that elevated expressions of AQP1 mRNA and VEGF protein were associated with increased volume of malignant pleural effusion. These results suggest that AQP1 and VEGF play important roles in the development of malignant pleural effusion in mice, which may help us find new strategies for the prevention and treatment of malignant pleural effusion.

Use of a Monte Carlo analysis within a physiologically based pharmacokinetic model to predict doxycycline residue withdrawal time in edible tissues in swine.

The pharmacokinetics of doxycycline were studied following a single intravenous (I.V.) and intramuscular (I.M.) injection of 10 mg/kg into eight healthy pigs. The steady-state tissue/plasma partition coefficients were obtained via a 3-h constant rate infusion (CRI) in four pigs. Based on the results of in vivo studies and the parameters derived from published work, a physiologically based pharmacokinetic (PBPK) model was developed to predict the drug concentration in edible tissues. The predicted values were then compared with those derived from a previous study. To account for individual differences in the processes of drug metabolism and/or diffusion, a Monte Carlo (MC) run of 1000 simulations was incorporated into the PBPK model to predict the doxycycline residue withdrawal times in edible tissues in swine. The withdrawal periods were compared with those derived from linear regression analysis. The PBPK model presented here provided accurate predictions of the observed concentrations in all tissues except for the injection site. The withdrawal times in all edible tissues derived from the MC analysis were longer than those from linear regression analysis. Based on the residues in the injection site and muscle tissue, the MC analysis predicted a withdrawal time of 33 days. Here, we illustrate that MC analysis can be incorporated into the PBPK model to accurately predict doxycycline residue withdrawal time in edible tissues in swine.