Quantitative systems pharmacology (QSP) modeling is an emerging computational medicine approach with growing applications and significance in modern drug development. QSP models are generally formulated based on multiscale disease mechanisms and drug-target interactions, which makes them capable of integrating multimodal data from the preclinical and clinical space. This also enables them to generate quantitative characterization of the dynamic disease progression as well as high-throughput predictions of drug-induced efficacy and toxicity signals. Therefore, QSP modeling and model-based virtual clinical trials have been widely implemented to guide drug development, in scenarios such as target identification and assessment, clinical trial design, evaluation of combination therapy and biomarkers, and personalized medicine. In US and Europe, QSP modeling has been developing rapidly in the past 10 years and is now an integral part of the model-informed drug development paradigm; however, in China it is still a nascent field. Here we will present a comprehensive review of the recent advancements of QSP and its impact in modern drug development through a number of case studies. This review will provide guidance for the future drug development efforts and the growth of QSP practice in China.

Polygonatum franchetii Hua is a medicinal plant endemic to China from Polygonatum Mill. The chloroplast genomes of two P. franchetii individuals sampled from two different habitats were sequenced by using the DNBSEQ-T7 high-throughput sequencing platform. After assembly and annotation, the two complete chloroplast genomes were characterized, and then comparative and phylogenetic analyses were performed with other published chloroplast genome sequences from Polygonatum. The whole chloroplast genomes of the two P. franchetii individuals were 155 942 and 155 962 bp in length, with a large single copy region (LSC, 84 670 and 84 722 bp), a small single copy region (SSC, 18 564 and 18 566 bp) and a pair of reverse repeats (IRa/IRb, 26 354 and 26 337 bp), respectively. Both of them contained 113 genes, including 79 protein-coding genes (PCGs), 30 transfer RNA (tRNA) genes, and 4 ribosomal RNA (rRNA) genes. Comparative analyses showed that the genome length, the guanine and cytosine (GC) content, genes content and order were highly conserved between the two P. franchetii individuals and among different Polygonatum species. The detected repeat sequences, including dispersed repeats, tandem repeats and simple sequence repeats (SSRs), were also relatively similar in types and positions, though showing a slightly difference in number. No significant expansion or contraction of the inverted repeat regions was found. Sequences variation between the two P. franchetii individuals was lower than that among different Polygonatum species. Besides, coding sequences (CDS) showed less divergence than noncoding sequences, and sequence divergence of IRs regions was lower than that of the LSC and SSC regions, both intraspecifically and interspecifically. Eight sequences with high nucleotide diversity among different species were screened, all of which were found located in the LSC and SSC regions. Phylogenetic inference showed that all Polygonatum species clustered into a monophyletic clade with a 100% bootstrap value, within which, species in section Verticillata formed a distinct group, section Sibirica and section Polygonatum were sister groups. The two P. franchetii individuals grouped together and showed the closest phylogenetic affinity to P. stenophyllum Maxim., belonging to the section Verticillata. The chloroplast genome of P. franchetii and its phylogenetic position in Polygonatum were comprehensively investigated and clearly elucidated in this study, the results may lay a foundation for the resource development and utilization of P. franchetii, as well as further molecular identification and phylogenetic studies of medicinal Polygonatum species.

Anthocyanidin reductase (ANR) is one of the key enzyme in the flavonoid biosynthetic pathway, and its catalytic activity is important for the synthesis of plant anthocyanin. In this study, specific primers were designed according to the transcriptome data of Lonicera japonica Thunb., and the CDS, gDNA and promoter sequences of ANR genes from Lonicera japonica Thunb. and Lonicera japonica Thunb. var. chinensis (Wats.) Bak. were cloned. The results showed that the CDS sequences of LjANR and rLjANR were 1 002 bp, the gDNA sequences were 2 017 and 2 026 bp respectively, and the promoter sequences were 1 170 and 1 164 bp respectively. LjANR and rLjANR both contain 6 exons and 5 introns, which have the same length of exons and large differences in introns. The promoter sequences both contain a large number of light response, hormone response and abiotic stress response elements. Bioinformatics analysis showed that both LjANR and rLjANR encoded 333 amino acids and were predicted to be stable hydrophobic proteins without transmembrane segments and signal peptides. The secondary structures of LjANR and rLjANR were predicted to be mainly consisted of α-helix and random coil. Sequence alignment and phylogenetic analysis showed that LjANR and rLjANR had high homology with Actinidia chinensis var. chinensis, Camellia sinensis and Camellia oleifera, and were closely related to them. The expression levels of LjANR and rLjANR were the highest in flower buds and the lowest in roots. The expression patterns at different flowering stages were similar, with higher expression levels in S1 and S2 stages and then gradually decreased until reaching the lowest level in S4 stage, after a slow increase in S5 stage, the expression levels decreased again. The expression levels of ANR genes in the two varieties showed significant differences in roots, S2 and S5 stages, while the differences in stems, flower buds, S1, S3 and S6 stages were extremely significant. The prokaryotic expression vector pET-32a-LjANR was constructed for protein expression. The target protein was successfully expressed of about 59 kD. This study lays a foundation for further study on the function of ANR gene and provides theoretical guidance for breeding new varieties of Lonicera japonica Thunb.

Hemagglutinin and neuraminidase, two important glycoproteins on the surface of influenza virus, play a considerable role in the entry and release stage of the viral life cycle, respectively. With in-depth investigation of influenza virus glycoproteins and the continuous innovation of drug discovery strategies, a new generation of glycoproteins inhibitors have been continuously discovered. From the point of view of medicinal chemistry, this review summarizes the current advances in seeking small-molecule inhibitors targeting influenza virus glycoproteins, hoping to provide valuable guidance for future development of novel antiviral drugs.

In this study, the mechanism of Xiaoyan Lidan formula (XYLDF) against 3, 5-diethoxycarbonyl-1, 4-dihydro-2, 4, 6-collidine (DDC)-induced chronic intrahepatic cholestasis (CIHC) in mice was investigated based on metabolomics, molecular docking and pharmacological methods. In the pharmacodynamics study, a dosage of 5 g·kg^-1 (clinical equivalent) XYLDF was administered in DDC-induced mice, then the effect of XYLDF against CIHC was evaluated by measuring the levels of alanine aminotransferase (ALT), aspartate aminotransferase (AST), alkaline phosphatase (AKP) as well as total bilirubin (TBIL) in serum and observing liver histopathological changes. All experiments were approved by the Ethical Committee Experimental Animal Center of Guangzhou University of Chinese Medicine (ZYD-2021-001). The serum metabolites of mice in each group were detected and identified based on ultra-performance liquid chromatography quadrupole time-of-flight tandem mass spectrometry, and the relevant biological pathways and molecular key targets were further enriched. Molecular docking technology was used to further evaluate the binding activity of the main active ingredients of XYLDF with potential targets. Subsequently, the in vitro experiment was conducted for the validation of the vital target. The results showed that compared with the model group, XYLDF significantly decreased the levels of ALT, AST, AKP and TBIL in the serum of CIHC mice, as well as alleviated inflammatory infiltration and hepatocyte necrosis in liver tissue. According to the metabonomic study, a total of 35 differential metabolites was identified as biomarkers associated with cholestasis, 12 of which were significantly recovered by XYLDF treatment. These biomarkers were involved in the pathways of primary bile acid biosynthesis and linoleic metabolism, which are closely related to the mechanism of XYLDF against CIHC. Protein-protein interaction network indicated that cytochrome P450 3A4 (CYP3A4) and cytochrome P450 1A1 (CYP1A1) are significant potential targets with good binding properties with six major active ingredients of XYLDF. Furthermore, it was found that 4-methoxy-5-hydroxycanthin-6-one, dehydroandrographolide and isodocarpin, three of the main active components in XYLDF, markedly induced the expression of CYP3A4 mRNA in vitro. This study revealed that XYLDF mainly mediates the biosynthesis of bile acids in CIHC mice to improve liver tissue lesions and bile efflux disorders, among which, CYP3A4 is the key target in the protection of XYLDF against CIHC. This research provides a reference for further elucidation of the pharmacological mechanism of XYLDF.

In this study, we synthesized six tetrazine-dipyrromethene boron difluoride (BODIPY) probes and achieved a remarkable up to 14-fold increase in singlet oxygen yield via tetrazine bioorthogonal click-to-release reactions. We systematically investigated the photodynamic activity of these probes, revealing crucial structure-activity relationships. Additionally, we evaluated the stability and release kinetics of these probes and identified P5 and P6 as ideal candidates for photodynamic therapy in live cells. This innovative strategy opens new avenues for fine-tuning the photodynamic properties of BODIPY dyes, thereby expanding their utility in cancer therapy.

We performed an extensively targeting metabolomic detecting using ultra-performance liquid chromatography-mass spectrometry (UPLC-MS/MS) to compare the secondary metabolites in Dang shen [Codonopsis pilosula (Franch.) Nannf.] from Shanxi and Gansu provinces. The findings showed that 161 secondary metabolites in 6 groups (phenolic acids, flavonoids, lignans and coumarins, alkaloids, terpenoids, others) were found from Dang shen in Changzhi city of Shanxi province and Dingxi city of Gansu province. There were 98 secondary metabolites which is differed significantly. In comparison to Dingxi city, 33 different secondary metabolites of Dang shen in Changzhi city had a greater relative content, whereas relative content of 65 different metabolites in Dingxi city was higher. Metabolic pathway enrichment analysis revealed that phenolic acids and flavonoids were significantly different in the secondary metabolites of Dang shen from different producing places. This may be one of the reasons for the difference in the quality of Dang shen in Shanxi and Gansu provinces. This work compared and analyzed the secondary metabolites of Dang shen from Dingxi city in Gansu province and Changzhi city in Shanxi province for the first time, which lays the foundation for further study on the quality of Dang shen.

This study primarily concentrated on scientific problems of poor taste caused by unclear critical quality attributes of oral preparations manufactured by Chinese materia medica, successfully established an identification method for taste critical quality attribute and a taste improvement method combining electronic tongue with human senses, and determined the optimal taste formula, to improve patients' oral medication compliance. The study received ethical approval from the Review Committee of the Beijing University of Chinese Medicine. The results showed that the proportion of bitterness of Xiaoer Qingrening Granule was 61.8%, and its bitterness grade was 3.70, it was determined that bitterness is the critical quality attribute that caused the poor taste of Xiaoer Qingrening Granule. Additionally, the optimal taste formula per milliliter of Xiaoer Qingrening sugar-free intermediate was determined with allowable daily intake, solubility, and sweetness as the limiting conditions, which was 40 mg hydroxypropyl β-cyclodextrin, 180 mg trehalose, and 1.5 mg acesulfame potassium. Compared with the Xiaoer Qingrening Granule, the sensory evaluation score of the optimal taste formula was increased by 37.5 points. In conclusion, this study achieved the taste improvement of Xiaoer Qingrening Granule and formed a set of taste improvement strategies including the identification of taste critical quality attribute, the selection of the type and dosage of corrigent, and the optimization of taste formula, which provided a thought reference for the taste improvement of other oral preparations and a new perspective for quality control of intelligent manufacturing of traditional Chinese medicines.

Currently, the resistance of first-line anti-tuberculosis drugs has made the prevention and treatment of tuberculosis increasingly difficult, posing a serious threat to global public health. Several studies have shown that efflux pumps are one of the important causes for bacteria to develop multi-drug resistance and extremely-drug resistance, and efflux pump inhibitors can inhibit the efflux of antibacterial drugs, thereby reducing bacterial drug resistance. Numerous natural products and synthetic compounds have been reported to possess efflux pump inhibitory activity, but they have not been applied in clinical settings because of their toxicity, pharmacokinetic properties, etc. Therefore, we summarized the efflux pump inhibitory activity, antimicrobial activity, and structure-activity relationships of reported efflux pump inhibitors against Mycobacterium tuberculosis in recent years, providing references for the development of new efflux pump inhibitors with better activity and lower toxicity.

Based on the dual needs of analgesia and anti-inflammation in trauma treatment, this study uses acetaminophen and moxifloxacin hydrochloride as active pharmaceutical ingredients and develops a composite bilayer tablet with a dual-phase drug release system by using binder jet 3D printing technology. Due to the complexity of the 3D printing process, there is an interaction between the various parameters. Through the optimization of the process, the relationship between the key process parameters can be determined more intuitively. In this study, the process of extended-release tablets was optimized to maintain the mechanical properties of the tablets while realizing the regulation of release. The full-factor experimental design of three central points 2³ was used to analyze the factors that significantly affect the quality attributes of extended-release tablets and the interaction between factors. The optimal extended-release process parameters were obtained by the response optimizer: the inkjet quantity of the printing ink was 10 (about 13.8 pL), the powder thickness was 180 μm, and the running speed was 360 mm·s^-1. The in vitro of release of 3D printed composite bilayer tablets showed that the in vitro of release of 3D printed tablets and commercially available tablets conformed to the Ritger-Peppas release model. The results of porosity showed that the immediate-release layer of the preparation has many pores and large pore size, and the dissolution of the immediate release layer within 15 min was greater than 85%. The internal pore size of the extended release layer is large, but it can still release slowly for up to 8 h, the mechanism may be related to the extended release of HPMC gelation. On the basis of verifying the rationality of the design goal of 3D printed composite bilayer tablets, this study also provides a theoretical basis for the preparation of 3D printing complex preparations.