Literature DB >> 36081821

Formulation and evaluation of SGLT2 inhibitory effect of a polyherbal mixture inspired from Ayurvedic system of medicine.

Ankit Kumar1, Anoop Singh Negi2, Ashutosh Chauhan3, Ravindra Semwal1, Rajnish Kumar2, Ruchi Badoni Semwal4, Randhir Singh5,6, Tushar Joshi7, Subhash Chandra8, Sunil Kumar Joshi9, Deepak Kumar Semwal10.   

Abstract

Background and aim: The ingredients viz., Artemisia roxburghiana, Cissampelos pareira, Stephania glabra, Drimia indica, Roylea cinerea, Tinospora sinensis and Curcuma longa of the present formulation are used to treat diabetes in the Indian traditional medical system. Adopting the concept of multiple herbal mixtures for better therapeutic effects from the ancient Ayurvedic text Sarangdhar Samhita, the present study aimed to develop a polyherbal formulation (PHF) of seven herbs and to evaluate its sodium-glucose cotransporter protein-2 (SGLT2) inhibitory effect on type 2 diabetic rats. Experimental procedure: Streptozotocin (STZ) (60 mg/kg) and nicotinamide (NAM) (120 mg/kg) were intraperitoneally administered to induce type 2 diabetes in Wistar rats. The animals were divided into 5 groups viz. normal control, diabetic control, positive control (dapagliflozin at 0.1 mg/kg) and two test groups (PHF at 250 and 500 mg/kg). Various parameters including blood glucose, serum glutamic pyruvic transaminase (SGPT), serum glutamic-oxaloacetic transaminase (SGOT), bilirubin, triglycerides and creatinine were measured. Results and conclusion: The treatment with PHF (250 and 500 mg/kg) showed a significant (p < 0.05) decrease in blood glucose levels by 56.37% and 58.17%, respectively. The levels of SGOT, SGPT and bilirubin were significantly reduced in PHF-fed diabetic rats. Histopathological examination revealed no major changes in the treated groups as compared to the normal control. The molecular docking study showed strong binding of β-sitosterol, insulanoline, warifteine, dehydrocorydalmine, taraxerol acetate, lupeol, corydalmine and luteolin to SGLT2 protein. The present study concludes that PHF has promising antidiabetic activity via inhibiting SGLT2 protein without showing any adverse effects.
© 2022 Center for Food and Biomolecules, National Taiwan University. Production and hosting by Elsevier Taiwan LLC.

Entities:  

Keywords:  Diabetes mellitus; Histopathology; Liver function test; Molecular docking; Polyherbal formulation; SGLT2 protein

Year:  2022        PMID: 36081821      PMCID: PMC9446025          DOI: 10.1016/j.jtcme.2022.03.003

Source DB:  PubMed          Journal:  J Tradit Complement Med        ISSN: 2225-4110


Introduction

Diabetes mellitus (DM), a metabolic chronic disorder, is one of the major health issues that reached its alarming stage and affected around 500 million people worldwide irrespective of their genders. This number is projected to reach 700 million by the year 2045. Out of the total diabetes cases, almost 90% belong to type 2 DM which is mainly managed by different oral hypoglycemic drugs including metformin, glipizide, acarbose and miglitol. However, the rest 10% cases belong to type 1 DM which needs insulin therapy to manage the glucose level. The global expenditure on diabetes and its related complications was estimated to be $760 billion in the adults' population of which almost 70% of expenditure belongs to the aged population of 50 years and above. Almost 90% of the existing traditional medicines worldwide are mainly obtained from medicinal plants and are commonly used in their crude forms. The most advantage of herbal treatment is its high therapeutic value with fewer adverse effects. Importantly, most of the herbal remedies including spices can be used without consulting any practitioner because their use is based on traditional knowledge that transfers from one generation to another. It is evident that the remedial values of medicinal plants are due to the presence of their secondary metabolites such as alkaloids, flavonoids, steroids and terpenoids. In addition, primary metabolites like carbohydrates, proteins, vitamins and lipids also played a significant role in human health mainly in body-building and energy production. The ethnobotanical knowledge had reported more than 800 medicinal plants particularly useful in DM. Most of these plants are also found in India and still in traditional practice in rural areas. For the past two decades, the medicinal plants are extensively studied for their efficacy against DM for their scientific validation and also to find out an alternate for the existing antidiabetic drugs which are known for their serious side effects. Based on the recent findings, these plants are capable to reduced blood glucose level via different mechanisms including improvement in insulin secretion and glucose utilisation, activation of insulin receptors, regeneration of pancreatic β-cells and inhibition of α-glucosidase, α-amylase and sodium-glucose transport protein 2 (SGLT2)., In addition to using as crude, purified molecules like berberine, curcumin, mangiferin, stevioside and capsaicin obtained from these plants have also been investigated for their antidiabetic effect and found more effective than their source extracts. Inhibition of SGLT2, a protein responsible for 80–90% glucose reabsorption in the kidney, is presently one of the key targets in the treatment of type 2 DM and diabetes-associated complications like diabetic kidney disease.12, 13, 14 However, possible side effects like urinary tract infections, foot amputations and kidney damage by existing SGLT2 inhibitors viz. canagliflozin, dapagliflozin and empagliflozin are a matter of serious concern., Hence, herbal medicine with high efficacy and the least adverse effect can be developed as a substitute for existing SGLT2 inhibitory drugs to treat type 2 diabetic patients. In this direction, a formulation was developed with seven traditional antidiabetic herbs viz. Artemisia roxburghiana (aerial parts), Cissampelos pareira (roots), Stephania glabra (tubers), Drimia indica (tubers), Roylea cinerea (aerial parts), Tinospora sinensis (stems) and Curcuma longa (rhizomes) to evaluate its SGLT2 inhibitory effect. Artemisia roxburghiana Wall. ex Besser (Compositae) is used as a traditional medicine in India for diabetes, rheumatic arthritis and malarial fever., Earlier scientific studies revealed its antiprotozoal, anthelmintic, anti-inflammatory and antidiabetic potential. Various secondary metabolites such as taraxerol acetate, lupeol, friedelan-3β-ol, friedelin, betulinic acid, betulin, apigenin, artemisinin, apigenin-7,4-dimethyl ether, α-selinene, α-copaene, α-gurjunene, α-eudesmol, curcumene, caryophyllene oxide, β-selinene, germacrene D, β-eudesmol, δ-cadinene, bicyclogermacrene, artemisinin, scopoletin and quercetin have been reported from this plant. Cissampelos pareira L. (Menispermaceae) is a popular traditional remedy for ulcers, wounds, arthritis, asthma, cholera, diarrhoea, cancer, inflammation, snakebite, malaria, rabies and diabetes. Its antidiabetic, antimalarial, antiparasitic, antiulcer, antidiarrhoeal and antivenom activities have been scientifically validated in different experimental models. Chemical analysis of its roots showed the presence of bulbocapnine, corydine, dicentrine, norrufescine, insulanoline, isochondrodendrine, warifteine, pelosine, cissampareine, insularine, cissampentin, hayatine, hayatinin, dehydrodicentrine, cycleanine, reserpine, sepeerine, berberine, cycleanine, corytuberine, pareitropone, cissamine, laudanosine and milonine. Stephania glabra (Roxb.) Miers (Menispermaceae) is also used traditionally in asthma, tuberculosis, dysentery, diabetes, cancer, fever and inflammation in many Asian countries., Similar to C. pareira, it is also rich in alkaloidal constituents like roemerine, reticuline, palmatrubine, dehydrocorydalmine, corydalmine, jatrorrhizine, stepharanine, capaurine, magnoflorine, stepholidine, n-desmethylcycleanine, tuduranine, tetrahydropalmatine, columbamine, laurifoline, stepharine, corynoxidine and palmatine. It is evident from the literature that tubers of this plant showed strong antidiabetic, antipsychotic, antibacterial and antihypertensive effects. The bulbs of Drimia indica (Roxb.) Jessop (Asparagaceae) are popular for their Ayurvedic properties in cardiac protection, GI disorders, asthma, leprosy, skin problems and diabetes., Cichoriol A, urgineol A-F, anhydroscilliphaeosidin, proscillaridin A, scillaren A, scilliphaeoside, scillarenin, β-sitosterol, luteolin, stigmasterol, campesterol, kaempferol and quercetin are the common constituents reported from this plant. Pharmacological studies on the bulbs revealed anticancer, antidiabetic, antimicrobial, anti-inflammatory and bronchodilator action on different models. Roylea cinerea (D. Don) Baillon (Lamiaceae) is traditionally used for treating diabetes, malaria, stomachache and skin diseases. The aerial parts of this plant exhibited different pharmacological properties including antidiabetic and anticancer. Its phytochemical profile showed the presence of 4-methoxybenzo[b]azet-2(1H)-one, 3β-hydroxy-35-(cyclohexyl-5′-propan-7′-one)-33-ethyl-34-methyl-bacteriohop-16-ene, stigmasterol, β-sitosterol, β-amyrin and cetyl alcohol., Tinospora sinensis (Lour.) Merr. (Menispermaceae) is well-known Ayurvedic medicine for Jaundice, fever, rheumatism, gonorrhoea and diabetes. Although it showed diverse pharmacological activities including antidiabetic and antipyretic, its role as an anti-severe acute respiratory syndrome coronavirus 2 (anti-SARS-CoV2) agent as well as a prophylaxis for coronavirus disease (COVID-19) has gained high popularity in India and hence considered as an essential Ayurvedic drug for COVID-19. Columbine, tinosporaside, jatrorrhizine, palmatine, berberine, tembeterine, tinocordifolioside, choline, tinospora acid, tinospora, tinosporin and tinosporide. Curcuma longa L. (Zingiberaceae) is traditionally used as a dietary spice and medicinal herb for different infections, wounds, gastrointestinal (GI) troubles and diabetes., Several studies reported its multiple pharmacological properties such as anti-inflammatory, antioxidant, antitumor, antibacterial, anticoagulant and antidiabetic activities. It contains curcumin, demethoxycurcumin, bisdemethoxycurcumin, turmerones, furanodiene, ferulic acid, coumaric acid, myristicin and cyclocurcumin as bioactive constituents. Since no established treatment is yet available for DM, research is continuously going on to find out a permanent solution mainly for type 2 DM. On the other hand, various earlier researches revealed that a single herb is sometimes insufficient to achieve the required medicinal effect. However, the combination of more than one herb, called polyherbalism, produces enhanced medicinal properties as well as less toxicity, perhaps due to the synergistic effect, and this remains the actual challenge in the development of polyherbal formulation. The herbs used in the present study are used in traditional medicine either as a single herb or in the form of a formulation.,39, 40, 41 Hence, the present research aimed to develop a polyherbal formulation inspired from Ayurveda to treat diabetes via inhibiting SGLT2 protein.

Materials and methods

Collection and authentication of plant material

The plant samples were collected from the different localities of the Garhwal and Kumaon regions of Uttarakhand Himalaya. The authentication of the plants was done by the Botanical Survey of India, Dehradun and R&D Centre, Uttarakhand Ayurved University, Dehradun. The accession/voucher numbers 118189, 118595, 118596, 118597, UAU211, UAU229 and UAU253 were assigned to Artemisia roxburghiana Wall. ex Besser, Drimia indica (Roxb.) Jessop, Cissampelos pareira L., Stephania glabra (Roxb.) Miers, Tinospora sinensis (Lour.) Merr., Roylea cinerea (D.Don) Baill. and Curcuma longa L., respectively.

Preparation of extracts

The freshly collected plant parts were thoroughly washed with plenty of tap water to remove impurities and put in the hot air oven at 50 °C for 48 h with occasional turns upside down. The dried and coarsely powdered material was soaked with distilled water and absolute ethanol (1:1) for 72 h (three times, each for 24 h) for complete extraction. The solvent was removed using a vacuum rotary evaporator (Equitron, India) with bath temperature at 50 °C to obtain a dry extract which was then stored at 4 °C until further use.

GC-MS/MS metabolomics analysis of extracts

GC-MS analysis of the extract was performed on Agilent 7890B gas chromatograph (Agilent Technologies, CA, USA) coupled with an Agilent 5977B mass detector. The sample was injected into GC-MS by an automatic sampler (CN1700443 3 series). HP-5 MS column (5% phenyl methyl polysiloxane; 30 m × 0.25 mm i.d. x 0.25 μm) and helium as a carrier gas were used for metabolites separation. For GC-MS analysis, the extract was dried in an Eppendorf concentrator and resuspended in 100% methanol and further dried in the concentrator, and then directly derivatized with 70 μL of N-methyl-N-(trimethylsilyl) trifluoroacetamide for 60 min at 37 °C in a 1.5 mL centrifuge tube followed by centrifugation at 13000 rpm for 10 min. The resulting supernatant was taken out and directly analyzed by the GC-MS system. The injection volume was 1 μL with splitless mode. The injector temperature was set at 280 °C. The oven temperature was initially set at 80 °C for 2 min, then ramped to 220 °C at a rate of 10 °C/min without any hold and further increased at 310 °C at the rate of 20 °C/min held for 10 min and with a solvent delay of 5 min. The column flow rate was 1 mL/min. The conditions for the operation of the mass spectrometer were set as follows: ion source temperature 230 °C, MS Quad temperature 150 °C, electron energy (70eV) and scanning range of m/z, 25–1000 amu. Metabolites were identified by matching the mass spectra of target metabolite (3:1 signal to noise ratio) with the NIST-17 mass spectral library. Metabolite identity was reported only when the matching value of the mass spectra comparison was more than 80%, and an increase in the area of the corresponding peak was observed when spiking the sample with the corresponding pure standard. The identified metabolites were further compared with the available literature on the plant.

Preparation of polyherbal formulation

The polyherbal formulation (PHF) was prepared by mixing all hydroalcoholic extracts in a definite ratio by following Bhavprakash Nighantu, an Ayurvedic classical text. The PHF contained 7.42% A. roxburghiana (aerial parts), 22.26% C. pareira (roots), 7.42% S. glabra (tuber), 0.15% D. indica (tuber), 7.42% R. cinerea (aerial parts), 29.69% T. sinensis (stems) and 25.64% Curcuma longa (rhizomes).

Experimental animals

Wistar rats with 80–120 g weight were selected for the study irrespective of their sex. The animal studies were conducted at M.M. College of Pharmacy, Ambala with ethical approval No. MMCP/IAEC/55/2019. All rats were acclimatized to the laboratory condition for one week before starting the study. The general behaviour, body weight, and feed-water intake of the rats were observed during the acclimatization period. The animals were maintained as per Control and Supervision of Experiments on Animals (CPCSEA) guidelines, and experimentation was done as per the approved protocol. The rats were kept on a permitted diet and water ad lib during the study.

Acute toxicity study

The rats were divided into 2 groups of 6 rats per group irrespective of their sex to assess the acute toxicity as per Organisation for Economic Co-operation and Development (OECD) guidelines. Group 1 served as a control and provided a normal diet whereas group 2 is given a single oral dose of PHF at 2000 mg/kg, body weight and also allowed a normal diet. All rats were monitored continuously every hour for the first day and then every day for fourteen days to assess the different parameters including skin, fur, eyes, respiratory pattern, salivation, diarrhoea, urination, tremors, ptosis, relaxation, gait, posture, lethargy, sleep, coma and food/water intake pattern together with other behavioural changes. However, the body weight of the rats was measured after every 4 days.

Induction of diabetes

The diabetes was induced in overnight fasted rats by a single intraperitoneal injection of 60 mg/kg of STZ (Sisco Research Laboratories, India), freshly prepared in 0.1 M citrate buffer (pH 4.5), 15 min after the intraperitoneal administration of 120 mg/kg of NAM (Sisco Research Laboratories, India), which was prepared in saline solution. The diabetes was confirmed by the elevated blood glucose level determined at 72 h after the STZ/NAM administration. The rats with blood glucose levels of more than 250 mg/dL were considered diabetic and used for the present study.

Experimental design

The rats were divided into five groups (n = 6) in which the first group comprised of non-diabetic rats whereas the rest of the groups included diabetic rats. The normal control group (G1NC) and diabetic control group (G2DC) did not receive any treatment but were allowed free access to distilled water. The third group (G3S) received standard drug dapagliflozin orally at the dose of 0.1 mg/kg/day in distilled water whereas the fourth group (G4T1) and fifth group (G5T2) received PHF at the oral doses of 250 and 500 mg/kg/day, respectively in distilled water for 28 days. The doses of PHF were determined based on acute toxicity study.

Biochemical analysis

The blood samples were collected from retro-orbital plexus using ether anaesthesia. Blood collected into the Clot Activator tubes or Fluoride tubes (for glucose test) was stored at 4 °C for 3 h and then centrifuged at 3000 revolutions per minute for 10 min using MiniSpin Plus Microcentrifuge (Eppendorf, Germany) to separate the serum. The isolated serum samples were used to measure the levels of fasting blood glucose (FBG), SGPT, SGOT, total bilirubin, cholesterol, triglyceride, creatinine, uric acid, urea and lactate dehydrogenase (LDH). All the parameters were assessed routinely on the 0th, 7th, 14th, 21st and 28th day using Standard Assay Erba Diagnostic Kits (Erba Semi-Auto Biochemistry Analyzer, UK). In addition, body weight of rats was also measured regularly using digital weighing balance.

Isolation of organs and their histopathology

At the end of the 28th day, the rats were painlessly sacrificed by cervical dislocation. The liver, kidney and pancreas of each rat were carefully removed for histopathological examination. The organs were washed with normal saline and preserved in a 5% formalin solution (Rankem, India) at room temperature until further use. The histopathological examination of isolated organs was performed by Dr. Lal Pathlabs Ltd., New Delhi, India using a microscope at 10x, 20x, 40x and 100x magnifications.

Molecular modelling for SGLT2 protein

Due to the unavailability of three dimensional (3-D) structure of SGLT2 of Homo sapiens, the sequence of SGLT2 was retrieved from NCBI (accession number- NP_003032.1), and the 3D structure was predicted using the “Easy Modeller” homology modelling module. The generated 3D structures of SGLT2 were subjected to further optimization, which includes minimization, side chain as well as loop refinement. The optimized structure was submitted into the PDBsum server of protein data bank (PDB) to validate the modelled structure.

Active site confirmation

After modelling, the computed atlas of surface topography of proteins (CASTp) webserver was used to find out the active side in protein and locate delineating, and measuring geometric and topological properties of the modelled SGLT2 structure. The CASTp 3.0 predicted many active amino acid residues which may present in the active site and are responsible for protein-ligand interaction.

Molecular docking and visualization

Docking was performed to obtain a population of possible orientations and conformations for the ligand at the binding site by using InstaDock open-source software. This software performs the prediction of the bound conformations based on the binding affinity. The grid centre for docking analysis was set to X = −25.462, Y = -43.18, and Z = 61.307, and the dimensions of the grid box were set as 84 × 110 × 99 Å having a spacing of 0.375 Å between the grids points. The virtual screening of compounds was conducted by rigid molecular docking in the active site of SGLT2 protein. A total of 170 compounds for molecular docking study were selected from our previous reports on the individual herbs of the present PHF and other antidiabetic agents.,,,,, Molecular interactions between protein-ligand complexes, including hydrogen bonds and other bonds, were analyzed and depicted using discovery studio software.

Statistical analysis

All the values are shown as means ± standard error of the mean (SEM). Statistical analyses were performed using one-way ANOVA followed by Tukey's test and analyzed using the GraphPad Prism 8 software. The p values equal to or less than 0.05 were considered statistically significant.

Results

Identification of metabolites

Based on the GC-MS/MS analysis, most intense peaks with highest peak width were characterized from each methyl derivatized extract. The results revealed 2-pyrrolidinone; L-arabitol; α-lyxopyranose; palmitic acid; bisphenol A, bis(trimethylsilyl) ether and acethydrazide, N2-[4-(thiitan-3-yloxy)benzylideno]- as main constituents of the extract of C. pareira. A. roxburghiana extract was found to contain protocatechoic acid; caffeic acid; arachidonic acid and azelaic acid. Similarly, 2-methylpropanoic acid; 2-hexenedioic acid, bis(trimethylsilyl) ester; d-tagatofuranose, pentakis(trimethylsilyl) ether; acetin, bis-1,3-trimethylsilyl ether and levoglucosan were characterized from D. indica extract. In case of T. sinensis methyl derivatized extract, 2-butenedioic acid; 5α-androstan-3β-ol, O-methyl-; ((11,11-dimethyl-8-methylenebicyclo[7.2.0]undecan-4-yl)methoxy)trimethylsilane; 2,5-dichlorohydroquinone; ethanesulfonic acid, 2-[(trimethylsilyl)oxy]-, trimethylsilyl ester; syringaresinol and β-epinephrine were found as major constituents. However, the extract of S. glabra comprised of cyclopropane, 1-(1-hydroxyethyl)-1-(diethylphosphonyl)-2-methylene-; 5α-androstan-3β-ol, O-methyl-; pyridine-3-sulfonamide, 2,6-dichloro-4-methyl-N-cyclohexyl-; 6H-dibenzo[a,g]quinolizine, 5,8,13,13a-tetrahydro-2,3,9,10-tetramethoxy-; 11-hydroxyetiocholanolone (3α,5α,11β) and 2-{2-oxo-4-[(trimethylsilyl)oxy]pyrrolidin-1-yl}-N-(trimethylsilyl)acetamide as main constituents.

Acute toxicity

All the animals treated with a single oral dose of PHF at 2000 mg/kg showed normal behavioural, motor, and neuronal functions. No mortality was observed during the experimental period up to 14 days. The monitoring of skin, fur, eyes, respiratory pattern, autonomic nervous system (ANS) characteristics (i.e., salivation, diarrhoea and urination) and central nervous system (CNS) characteristics (i.e., tremors, ptosis, relaxation, gait and posture) in the treated rats remained unaffected. Moreover, the water and feed intake pattern of the rats was regular and consistent during the study period. The changes in body weight of the treated group did not record any substantial changes when compared with the control. Hence, based on the acute toxicity outcome, PHF can be considered safe up to 2000 mg/kg body weight.

Effect on body weight

The weekly measurement of body weight showed no significant changes in diabetic control and standard group. However, the formulation with a higher dose showed a slight increase in the body weight of diabetic rats when compared with the control group. The average body weight of the animals treated with 500 mg/kg was recorded to be 100 ± 8.003 g on the 28th day which was increased from 87 ± 7.278 on a basal day whereas, at a dose of 250 mg/kg, the body weight increased to 100 ± 3.921 from 94 ± 4.416.

Effect on blood glucose level

After confirmation of hyperglycemia (fasting glucose level >250 mg/dL), the rats were treated with PHF and standard drug. The results (Fig. 1) showed that blood glucose levels of all treated groups were significantly (p < 0.05) decreased on the 7th day of treatment and remains effective till the 28th day when compared to the normal control group. The PHF showed 56.37% and 58.17% reduction in glucose levels after 28 days at the doses of 250 and 500 mg/kg, respectively whereas 53.90% reduction was recorded with 0.1 mg/kg of dapagliflozin.
Fig. 1

Effect of PHF on blood glucose level of rats. G1NC (normal control group), G2DC (diabetic control group), standard group (G3S), test group 1 (G4T1) and test group 2 (G5T2); (∗) shows significance level with p value of <0.05 when compared to G2DC.

Effect of PHF on blood glucose level of rats. G1NC (normal control group), G2DC (diabetic control group), standard group (G3S), test group 1 (G4T1) and test group 2 (G5T2); (∗) shows significance level with p value of <0.05 when compared to G2DC.

Effect on liver function

The STZ-NAM-induced hyperglycemia significantly increased the SGPT and SGOT levels in rats from the 7th to the 28th day. However, there were no significant changes occurred in the levels of SGPT and SGOT in the standard and PHF-treated groups for up to 14 days. Interestingly, 14th day onwards, a significant (p < 0.05) decrease was noted till the end of the experiment (Table 1). Similarly, total bilirubin levels were increased in the diabetic control group during the entire experimental period while the changes measured in the treated groups were found insignificant till 14 days and significantly increased thereafter. Hence, it can be said that the PHF has long-term significant effects on the levels of SGPT, SGPT and total bilirubin in STZ-NAM-induced diabetic rats.
Table 1

Effect of PHF on liver function in experimental rats.

DayG1NCG2DCG3SG4T1G5T2
SGPT (U/L)
Basal (0) Day28.40 ± 4.6932.10 ± 5.3534.69 ± 3.2532.25 ± 6.1637.28 ± 0.49
7th Day31.54 ± 5.7741.77 ± 9.7943.30 ± 7.2736.86 ± 5.2947.14 ± 8.92
14th Day31.45 ± 5.0576.02 ± 16.5487.51 ± 12.3054.16 ± 7.8058.75 ± 7.68
21st Day33.78 ± 4.74119.53 ± 8.4460.76 ± 5.13∗47.28 ± 15.81∗53.51 ± 13.84∗
28th Day32.49 ± 3.3698.64 ± 11.8040.28 ± 15.03∗40.39 ± 6.84∗48.31 ± 7.45∗
SGOT (U/L)
Basal (0) Day89.10 ± 5.6092.82 ± 10.97103.47 ± 12.56101.48 ± 14.8689.09 ± 4.63
7th Day86.95 ± 7.84129.61 ± 6.72119.05 ± 33.85116.36 ± 8.33112.38 ± 3.40
14th Day92.82 ± 10.37169.17 ± 8.87140.74 ± 8.26137.20 ± 7.92113.84 ± 10.44∗
21st Day86.03 ± 12.70147.52 ± 15.57138.70 ± 13.21123.69 ± 13.05101.70 ± 5.89∗
28th Day81.33 ± 10.25181.20 ± 7.29119.72 ± 9.46∗115.34 ± 14.68∗118.35 ± 8.37∗
Total bilirubin (mg/dL)
Basal (0) Day0.66 ± 0.140.79 ± 0.160.62 ± 0.100.83 ± 0.090.88 ± 0.10
7th Day0.87 ± 0.350.94 ± 0.160.86 ± 0.220.91 ± 0.300.92 ± 0.13
14th Day0.68 ± 0.100.82 ± 0.220.73 ± 0.180.89 ± 0.090.84 ± 0.09
21st Day0.88 ± 0.060.94 ± 0.090.53 ± 0.07∗0.71 ± 0.08∗0.77 ± 0.23∗
28th Day0.87 ± 0.081.03 ± 0.220.69 ± 0.46∗0.79 ± 0.30∗0.88 ± 0.32∗

G1NC (normal control group), G2DC (diabetic control group), standard group (G3S), test group 1 (G4T1) and test group 2 (G5T2); (∗) Shows significance level with p value of <0.05 when compared to G2DC.

Effect of PHF on liver function in experimental rats. G1NC (normal control group), G2DC (diabetic control group), standard group (G3S), test group 1 (G4T1) and test group 2 (G5T2); (∗) Shows significance level with p value of <0.05 when compared to G2DC.

Effect on lipid profile

The levels of cholesterol and triglycerides in the serum of experimental rats were measured to evaluate the effect of PHF on lipid profile. Their levels were found higher in the diabetic control group when compared to the normal control group. The results (Table 2) revealed a significant (p < 0.05) change after two weeks in the levels of cholesterol and triglycerides in diabetic rats treated with both PHF and dapagliflozin when compared to the diabetic control group.
Table 2

Effect of PHF on lipid profile of experimental rats.

DayG1NCG2DCG3SG4T1G5T2
Cholesterol (mg/dL)
Basal (0) Day36.64 ± 8.8340.52 ± 4.7945.15 ± 3.4639.81 ± 6.9038.46 ± 9.06
7th Day41.56 ± 4.7743.06 ± 1.6846.76 ± 5.4050.5 ± 10.9843.79 ± 4.98
14th Day41.37 ± 5.3183.28 ± 9.6559.85 ± 11.49∗44.13 ± 11.69∗45.39 ± 7.64∗
21st Day37.51 ± 4.8174.81 ± 15.9959.86 ± 11.9844.7 ± 5.27∗48.00 ± 6.92∗
28th Day36.79 ± 12.2380.05 ± 4.7256.67 ± 7.29∗40.9 ± 3.59∗45.02 ± 11.02∗
Triglycerides (mg/dL)
Basal (0) Day37.61 ± 5.3237.13 ± 2.8934.81 ± 8.9440.61 ± 7.0047.11 ± 57.55
7th Day37.81 ± 14.8642.56 ± 9.2448.07 ± 13.5843.74 ± 7.3942.87 ± 3.27
14th Day44.58 ± 18.5145.84 ± 10.9542.30 ± 17.5549.62 ± 22.6841.13 ± 17.62
21st Day40.68 ± 15.5880.24 ± 22.1248.47 ± 21.23∗42.55 ± 30.34∗48.20 ± 12.40∗
28th Day47.92 ± 7.8086.86 ± 13.1542.06 ± 16.00∗44.47 ± 12.43∗40.74 ± 9.17∗

G1NC (normal control group), G2DC (diabetic control group), standard group (G3S), test group 1 (G4T1) and test group 2 (G5T2); (∗) Shows significance level with p value of <0.05 when compared to G2DC.

Effect of PHF on lipid profile of experimental rats. G1NC (normal control group), G2DC (diabetic control group), standard group (G3S), test group 1 (G4T1) and test group 2 (G5T2); (∗) Shows significance level with p value of <0.05 when compared to G2DC.

Effect on renal function

The present study showed that dapagliflozin and PHF did not produce any significant changes in the levels of creatinine, uric acid and urea in serum of diabetic rats. A slight increase in the levels of urea was although recorded in all treated groups but it was comparatively lower than that of the diabetic control group (Table 3).
Table 3

Effect of PHF on renal function of experimental rats.

DayG1NCG2DCG3SG4T1G5T2
Creatinine (mg/dL)
Basal (0) Day0.69 ± 0.020.68 ± 0.020.70 ± 0.050.73 ± 0.030.67 ± 0.06
7th Day0.74 ± 0.140.73 ± 0.100.73 ± 0.190.76 ± 0.100.69 ± 0.13
14th Day0.66 ± 0.050.75 ± 0.160.70 ± 0.040.74 ± 0.050.65 ± 0.08
21st Day0.84 ± 0.170.89 ± 0.190.80 ± 0.130.82 ± 0.070.74 ± 0.12
28th Day0.70 ± 0.020.85 ± 0.070.70 ± 0.020.68 ± 0.050.68 ± 0.09
Uric acid (mg/dL)
Basal (0) Day2.26 ± 0.382.62 ± 0.222.21 ± 0.652.46 ± 1.122.44 ± 0.48
7th Day3.58 ± 1.903.70 ± 0.792.81 ± 1.062.90 ± 1.833.45 ± 0.72
14th Day2.89 ± 0.783.38 ± 0.502.91 ± 0.942.23 ± 0.492.93 ± 0.44
21st Day2.87 ± 0.293.46 ± 0.872.48 ± 1.083.21 ± 1.092.91 ± 0.48
28th Day2.19 ± 0.522.58 ± 0.612.48 ± 0.392.18 ± 0.532.74 ± 1.10
Urea (mg/dL)
Basal (0) Day82.89 ± 14.2870.40 ± 17.9985.49 ± 12.6781.87 ± 19.2280.27 ± 13.20
7th Day86.12 ± 11.8480.96 ± 3.45111.40 ± 3.6791.87 ± 8.9782.33 ± 6.49
14th Day83.47 ± 6.3680.10 ± 8.8090.43 ± 3.5497.93 ± 1.9888.37 ± 6.71
21st Day82.47 ± 3.1888.87 ± 4.0196.71 ± 8.5697.81 ± 6.7185.80 ± 5.69
28th Day85.87 ± 18.68118.58 ± 18.1299.32 ± 11.1197.55 ± 14.5292.34 ± 7.60

G1NC (normal control group), G2DC (diabetic control group), standard group (G3S), test group 1 (G4T1) and test group 2 (G5T2); (∗) Shows significance level with p value of <0.05 when compared to G2DC.

Effect of PHF on renal function of experimental rats. G1NC (normal control group), G2DC (diabetic control group), standard group (G3S), test group 1 (G4T1) and test group 2 (G5T2); (∗) Shows significance level with p value of <0.05 when compared to G2DC.

Effect on LDH level

The diabetic control group showed a significant increase in their LDH level on the 28th day when compared to the normal control group. The treated groups also showed an increase in the LDH level but this change was found insignificant (Fig. 2). There was no significant difference in the effect of PHF (at 250 and 500 mg/kg) and dapagliflozin (0.1 mg/kg).
Fig. 2

Effect of PHF on LDH level in experimental rats. G1NC (normal control group), G2DC (diabetic control group), standard group (G3S), test group 1 (G4T1) and test group 2 (G5T2); the results were found insignificant when compared to G2DC.

Effect of PHF on LDH level in experimental rats. G1NC (normal control group), G2DC (diabetic control group), standard group (G3S), test group 1 (G4T1) and test group 2 (G5T2); the results were found insignificant when compared to G2DC.

Histopathology of isolated organs

Histopathology of liver

Microscopic examination of sections revealed the normal histological structure of hepatic lobules in both treated and untreated animals including the normal control group. The hepatocytes were arranged in cords radiating from the central vein and separated by blood sinusoids in all the groups. Mild sinusoidal congestion, as well as congestion of central and portal vein, was seen in G2DC and G5T2 groups similar to the G1NC group, whereas mild congestion of sinusoids was seen in the G4T1 group. Hepatocytes were polyhedral in shape with slightly vacuolated granular cytoplasm and vesicular nuclei in all the groups similar to the normal control group. Focal area of mild degenerative changes with swollen hepatocytes and some cells with karyolitic nuclei were present in G5T2 similar to the G1NC group, whereas the focal area of mild degenerative changes with mild steatosis and some cells with karyolitic nuclei were present in G4T1. Blood sinusoids separating the hepatic cords lined by endothelial cells and kupffer cells were seen in the G3S group only. Bi-nucleated hepatocytes with no infiltration or granuloma were also seen in all the groups (Fig. 3).
Fig. 3

Histopathology of right kidney (KR), left kidney (KL), pancrease (P) and liver (L) of experimental rats. (↓): Mild sinusoidal congestion; (□): mild degenerative changes; (⌂): moderate degenerative changes; (○) high degenerative changes; (◯) depletion of beta cells; (Δ): desquamation of epithelial linings; G1NC (normal control group), G2DC (diabetic control group), standard group (G3S), test group 1 (G4T1) and test group 2 (G5T2); the results were found insignificant when compared to G2DC.

Histopathology of right kidney (KR), left kidney (KL), pancrease (P) and liver (L) of experimental rats. (↓): Mild sinusoidal congestion; (□): mild degenerative changes; (⌂): moderate degenerative changes; (○) high degenerative changes; (◯) depletion of beta cells; (Δ): desquamation of epithelial linings; G1NC (normal control group), G2DC (diabetic control group), standard group (G3S), test group 1 (G4T1) and test group 2 (G5T2); the results were found insignificant when compared to G2DC.

Histopathology of pancreas

The section of G1NC showed normal architecture of the pancreas with very mild degenerative changes. The exocrine pancreas is composed of closely packed acinar cells and arranged into small lobules. Pancreatic lobules are separated by intact intra-lobular and interlobular connective tissue septa. The islet cells were seen interspersed between the acinar cells. Islets appeared lightly stained than the surrounding acinar cells. However, the histopathology of the pancreas of G2DC revealed moderate to high degenerative changes in both exocrine and endocrine components. Acinar cells were swollen and small vacuoles were observed in acinar cells. Interlobular ducts were lined with flattened epithelium. Some lobules showed depletion of islet β-cells. On the other hand, all the treated groups showed mild degenerative changes in both exocrine and endocrine components similar to G1NC. In addition, a few acinar cells swollen with small vacuoles were observed. Interlobular ducts were lined with flattened epithelium. A small number of lobules in G3S, G4T1 and G5T2 showed depletion of islet β-cells (Fig. 3).

Histopathology of right kidney

Microscopic examination of sections of the right kidney revealed mild haemorrhages present in between intertubular spaces in all the groups. Mildly congested blood vessels were also seen in all the groups. There are mild to moderate degenerative changes of epithelial linings of tubules seen in G2NC whereas all treated groups showed mild degenerative changes similar to the normal control group. Some tubules show desquamation of epithelial linings in G2DC and G3S groups. However, focal interstitial inflammatory mononuclear cell infiltrate were seen in G2DC but not in G3S. Similar to G1NC, no granuloma was seen in the rats treated with standard and PHF (Fig. 3).

Histopathology of left kidney

Microscopic examination of the left kidney showed most histological features similar to the right kidney. G5T2 showed focal mild intertubular haemorrhages. An infiltration of polymorph cells was present in surrounding fibro-fatty tissues in G5T2. Small cystic areas were only present in the G4T1 group. In the G3S group, some of the tubules showed desquamation of epithelial linings and mild expansion of the glomerular cavity similar to the normal control group. Additionally, no infiltration of inflammatory cells or granuloma was seen in the G3S group. G1NC and G2DC showed focal interstitial inflammatory mononuclear cell infiltrates. Focal tubular basophilia, nuclear crowding, and thickened basement membranes were also present in the G2DC group. In the treatment group G4T1, tubules showed necrosis of epithelial linings (Fig. 3).

Homology modelling and active site confirmation

The SGLT2 protein sequence of Homo sapiens was retrieved from NCBI, and a basic local alignment search tool-protein (BLASTp) was used to identify homologous templates from the PDB database. BLASTp showed similarities of the SGLT2 query with the crystal structure of secretogranin-2 (SGC2) of Vibrio parahaemolyticus (PDB ID 2XQ2 A, 3DH4 and 2XQ2 B), and these structures share 69%, 64% and 63% identity respectively query cover. After that, Easy Modeller software was used for the prediction of the 3D structure of SGLT2 using four structures as templates. The 3-D model of protein also showed 90.5% of the residues in the allowed region of the Ramachandran plot. It demonstrates the stereochemical reliability of the homology template of SGLT2 protein. The active site (Fig. 4) residue numbers of SGLT2 of Homo sapiens are 32, 33, 36, 37, 40, 41, 43, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 65, 66, 69, 70, 72, 73,74, 75, 76, 77, 78, 79, 80, 83, 84, 95, 96, 98, 99, 102, 123, 124, 125, 126,127,129, 130, 136, 139, 142, 143, 145, 146, 149, 150, 153, 154, 157, 158, 161, 184, 187, 188, 189, 190, 191, 192, 193, 194, 195, 196, 197, 198, 200, 201, 202, 203, 204, 205, 208, 209, 283, 286, 287, 289, 290, 291, 294, 295, 296, 297, 300, 301, 304, 320, 321, 324, 329, 387, 388, 389, 390, 391, 392, 393, 395, 396, 397, 398, 399, 400, 402, 403, 404, 406, 407, 411, 412, 414, 415, 416, 417, 418, 419, 420, 422, 423, 424, 426, 427, 430, 452, 453, 456, 460, 463, 464, 467, 471, 552, 553, 659, 662, 663, 666, 667, 670, 671, 672. Molecular docking was performed with SGLT2 of Homo sapiens, specifying the active site using the amino acids area.
Fig. 4

Active sites of SGLT2 protein.

Active sites of SGLT2 protein.

Molecular docking

Phytochemicals of plants were screened by molecular docking against the SGLT2 receptor. Out of 170 compounds, β-sitosterol and insulanoline showed the highest binding energy. The top eight active compounds with their docking scores, natural sources, 2D and 3D molecular docking interactions with SGLT2 protein are shown in Fig. 5.
Fig. 5

2D and 3D molecular docking interactions of bioactive compounds with SGLT2 protein. DS (docking score - kcal/mol), BS (botanical source).

2D and 3D molecular docking interactions of bioactive compounds with SGLT2 protein. DS (docking score - kcal/mol), BS (botanical source).

Discussion

Based on the global research outcomes to date, it is understood that DM is incurable, although it can be managed with the help of selected medications including modern and traditional medicine. The treatment of DM with herbal medicine is considered a safer approach which has been practised in different traditional systems of medicine since ancient times. Moreover, the treatment with multiple herbs showed better effects with lesser toxicity than a single herb, perhaps due to the synergistic effect of different constituents present in these herbs., In the past couple of years, tremendous work has been done on the polyherbal formulations and many of them were found effective in controlling hyperglycemia in different experimental models.56, 57, 58 The evaluation of toxicity of a plant sample is a key step for conducting pharmacological studies in higher models. Such studies give initial information about the safety profile of a drug candidate and also help in deciding the doses for in vivo studies. In the present study, PHF was evaluated for its acute toxicity in rats and based on the outcome of this experiment, the doses for antidiabetic activity were determined. Since the highest doses at 2000 mg/kg did not produce any toxic effect in rats up to 14 days, the oral LD50 of PHF can be considered to be higher than 2000 mg/kg. For the past couple of years, a combination of STZ and NAM is most accepted to develop type 2 DM because STZ causes β-cell damage while NAM partially protects the insulin-secreting cells against STZ and hence, protects them from permanent damage. The present study found STZ-NAM-induced diabetic model as the most suitable to evaluate the SGLT2 inhibitory effect in rats. The administration of STZ and NAM produces a significant elevation in the glucose levels of experimental rats when compared to the normal group. The reduction in blood glucose level is the first priority in developing an antidiabetic drug because hyperglycemia can cause various complications including retinopathy, nephropathy and neuropathy. In the present study, the treatment of PHF (250 and 500 mg/kg) significantly (p < 0.05) decreased the elevated blood glucose levels in STZ-NAM induced diabetes rats within one week and this effect was also recorded during the entire experimental period when compared with the diabetic control group. The study also found that the effect of PHF was comparable to that of the positive control group at the end of the 28th day. Moreover, body weight is an important parameter to be monitored regularly in the event of diabetes. However, in the present study, no significant changes were recorded in body weight of experimental rats. The increased levels of lipids, mainly cholesterol and triglycerides, in diabetic patients are risk factors for diabetes-associated cardiovascular diseases. It is evident that the risk of developing hypertension is higher in DM patients than in others. Moreover, chronic diabetes plays a key role in developing vascular diseases by increasing the formation of advanced glycation endproducts (AGEs), activation of the receptor for advanced glycation endproducts (RAGE), oxidative stress and inflammation. The results obtained from the present study revealed a significant change in the levels of total cholesterol and triglycerides both in test groups and positive control. The maximum reduction was recorded on the 28th day of the experiment. Similarly, the increased levels of SGPT, SGOT and bilirubin due to liver injury are common conditions mainly in chronic diabetes. This condition can also affect the metabolism of lipids, carbohydrates and proteins and hence causes fatty liver. Moreover, developing liver cirrhosis, hepatocellular carcinoma and oxidative stress were also recorded in the patients with chronic DM. In the present study, the PHF significantly decreased the levels of SGPT and SGOT enzymes after the 14th day. Moreover, total bilirubin levels were increased in the diabetic control group during the entire experimental period while the changes measured in the treated groups were found significant when compared to G2DC. Hence, it can be said that the PHF has a significant effect on the levels of SGPT, SGPT and total bilirubin in STZ-NAM-induced diabetic rats and as a result, it might play a potential role in the hepatoprotective activity. On the other hand, elevated levels of creatinine, urea and uric acid in the blood are indications of kidney problems. The present study did not find a significant change in the levels of these components when compared to the control group. The histopathological examination of isolated organs was done to find any internal toxic sign in the treated animals. The sample can be considered safe if its histology matches with normal animals. This study showed that PHF does not cause any harmful effect on the experimental rats when compared to diabetic rats. The results of PHF for liver and kidney function were supported by histopathological analysis. Moreover, the pancreas was also found unaffected after the 28-days of in vivo study. Overall, the formulation was found effective and safe for further use. Dapagliflozin, a synthetic antidiabetic drug, has the potential to reduce blood glucose specifically by inhibiting SGLT2 protein in the renal tubule to block glucose reabsorption. In addition to reducing HbA1c level, this drug also helps in reducing body weight. However, previous studies indicate that the SGLT2 inhibitors may produce adverse effects on the skeleton which results in the weakening of bones in diabetic patients. The present study used dapagliflozin as a standard drug to compare the antidiabetic effect of PHF. The results showed that the values of most of the parameters studied for PHF were comparable to the standard drug. Hence, in view of the potential glucose-lowering effect and safety profile, the PHF was also evaluated for its SGLT2 inhibitory effect via in silico method using major bioactives of the herbs present in PHF. Nowadays, molecular docking studies are most popular among the scientific community for the early inclusion of pharmacokinetics consideration in the process of novel drug discovery with the help of advanced software. The present study identified eight bioactive phytochemicals including β-sitosterol and insulanoline which showed strong binding to the SGLT2 protein with a docking score of more than nine. Since most of the constituents of PHF showed in silico SGLT2 inhibitory activity, it can be suggested that the glucose-lowering effect of PHF involved the SGLT2 inhibition pathway. Moreover, the individual molecules having high docking scores can be further studied using higher models.

Conclusion

This study concluded that the polyherbal formulation, comprised of 7 herbs viz. A. roxburghiana, C. pareira, S. glabra, D. indica, R. cinerea, T. sinensis and C. longa, showed significant glucose-lowering effect via inhibiting SGLT2 protein. The levels of SGOT, SGPT, bilirubin, cholesterol and triglycerides were significantly decreased in the PHF-fed rats. The toxicity study, as well as histopathological examination, revealed the nontoxic nature of the formulation. Hence, in view of the antidiabetic effects in animal models, the PHF can be developed as a novel antihyperglycemic agent. However, further studies are needed to consider it for clinical use.

Animal ethical approval

The animal studies were conducted at M.M. College of Pharmacy, Ambala with ethical approval No. MMCP/IAEC/55/2019. The animals were maintained as per Control and Supervision of Experiments on Animals (CPCSEA) guidelines, and experimentation was done as per the approved protocol.

Authors’ contribution

DKS: conceptualized and supervised the work, edited and reviewed the manuscript; AK, ASN, RK and RSD: developed the formulation and conducted the animal experiments; TJ and SC: conducted the molecular docking studies; AK and RS: prepared the first draft of the manuscript; RS and AC: statistically analyzed the data: RBS and SKJ: edited and reviewed the manuscript.

Declaration of competing interest

The authors declared no conflict of interest.
  48 in total

1.  Antidiabetic plants used by Sikkim and Darjeeling Himalayan tribes, India.

Authors:  D R Chhetri; P Parajuli; G C Subba
Journal:  J Ethnopharmacol       Date:  2005-04-11       Impact factor: 4.360

Review 2.  Efficacy and safety of Stephania glabra: an alkaloid-rich traditional medicinal plant.

Authors:  Deepak Kumar Semwal; Ruchi Badoni Semwal
Journal:  Nat Prod Res       Date:  2014-09-04       Impact factor: 2.861

Review 3.  Streptozotocin-nicotinamide-induced diabetes in the rat. Characteristics of the experimental model.

Authors:  Tomasz Szkudelski
Journal:  Exp Biol Med (Maywood)       Date:  2012-05-22

Review 4.  Protective and therapeutic effects of natural products against diabetes mellitus via regenerating pancreatic β-cells and restoring their dysfunction.

Authors:  Deepak Kumar Semwal; Ankit Kumar; Sonali Aswal; Ashutosh Chauhan; Ruchi Badoni Semwal
Journal:  Phytother Res       Date:  2020-09-28       Impact factor: 5.878

Review 5.  Sodium-Glucose Cotransporter 2 Inhibition and Diabetic Kidney Disease.

Authors:  Radica Z Alicic; Joshua J Neumiller; Emily J Johnson; Brad Dieter; Katherine R Tuttle
Journal:  Diabetes       Date:  2019-02       Impact factor: 9.461

6.  Development of an antidiabetic formulation (ADJ6) and its inhibitory activity against α-amylase and α-glucosidase.

Authors:  Anand Duraiswamy; Devanand Shanmugasundaram; Changam Sheela Sasikumar; Sanjay M Cherian; Kotturathu Mammen Cherian
Journal:  J Tradit Complement Med       Date:  2015-02-28

Review 7.  Phytomedicines Used for Diabetes Mellitus in Ghana: A Systematic Search and Review of Preclinical and Clinical Evidence.

Authors:  Michael Buenor Adinortey; Rosemary Agbeko; Daniel Boison; William Ekloh; Lydia Enyonam Kuatsienu; Emmanuel Ekow Biney; Obed O Affum; Jeffery Kwarteng; Alexander Kwadwo Nyarko
Journal:  Evid Based Complement Alternat Med       Date:  2019-04-18       Impact factor: 2.629

8.  Pharmacological evaluation of "sugar remedy," a polyherbal formulation, on streptozotocin-induced diabetic mellitus in rats.

Authors:  Sandeep Singhal; Arvind Singh Rathore; Vikram Lohar; Rakesh Dave; Jeetesh Dave
Journal:  J Tradit Complement Med       Date:  2014-07

9.  Antiproliferative Effects of Roylea cinerea (D. Don) Baillon Leaves in Immortalized L6 Rat Skeletal Muscle Cell Line: Role of Reactive Oxygen Species Mediated Pathway.

Authors:  Astha Bhatia; Harpal Singh Buttar; Rohit Arora; Balbir Singh; Amritpal Singh; Sarabjit Kaur; Saroj Arora
Journal:  Front Pharmacol       Date:  2020-03-13       Impact factor: 5.810

Review 10.  Potential Safety Issues with Use of Sodium-Glucose Cotransporter 2 Inhibitors, Particularly in People with Type 2 Diabetes and Chronic Kidney Disease.

Authors:  Tamara Y Milder; Sophie L Stocker; Richard O Day; Jerry R Greenfield
Journal:  Drug Saf       Date:  2020-10-23       Impact factor: 5.606

View more

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