CN107292126B - 一种中药对复杂性疾病所致“失和”网络整合调节作用的定量评价方法 - Google Patents
一种中药对复杂性疾病所致“失和”网络整合调节作用的定量评价方法 Download PDFInfo
- Publication number
- CN107292126B CN107292126B CN201710308991.3A CN201710308991A CN107292126B CN 107292126 B CN107292126 B CN 107292126B CN 201710308991 A CN201710308991 A CN 201710308991A CN 107292126 B CN107292126 B CN 107292126B
- Authority
- CN
- China
- Prior art keywords
- network
- traditional chinese
- chinese medicine
- node
- regulation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000003814 drug Substances 0.000 title claims abstract description 131
- 230000033228 biological regulation Effects 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 62
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 title claims abstract description 52
- 201000010099 disease Diseases 0.000 title claims abstract description 51
- 230000000694 effects Effects 0.000 title claims abstract description 28
- 238000011158 quantitative evaluation Methods 0.000 title claims abstract description 8
- 230000010354 integration Effects 0.000 title abstract description 6
- 108090000623 proteins and genes Proteins 0.000 claims abstract description 91
- 102000004169 proteins and genes Human genes 0.000 claims abstract description 58
- 229940126680 traditional chinese medicines Drugs 0.000 claims abstract description 33
- 238000003032 molecular docking Methods 0.000 claims abstract description 28
- 230000003334 potential effect Effects 0.000 claims abstract description 18
- 230000001105 regulatory effect Effects 0.000 claims abstract description 8
- 239000000523 sample Substances 0.000 claims description 37
- 206010000891 acute myocardial infarction Diseases 0.000 claims description 31
- 230000037361 pathway Effects 0.000 claims description 21
- 150000001875 compounds Chemical class 0.000 claims description 14
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 13
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 13
- 235000008434 ginseng Nutrition 0.000 claims description 13
- 230000014509 gene expression Effects 0.000 claims description 12
- 238000011156 evaluation Methods 0.000 claims description 10
- 230000009471 action Effects 0.000 claims description 9
- 230000019491 signal transduction Effects 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 239000000203 mixture Substances 0.000 claims description 8
- 230000003993 interaction Effects 0.000 claims description 7
- 239000003446 ligand Substances 0.000 claims description 7
- 208000024172 Cardiovascular disease Diseases 0.000 claims description 4
- 210000000941 bile Anatomy 0.000 claims description 4
- 230000003828 downregulation Effects 0.000 claims description 4
- 238000010201 enrichment analysis Methods 0.000 claims description 4
- 230000003827 upregulation Effects 0.000 claims description 4
- 241000208340 Araliaceae Species 0.000 claims 1
- 238000011160 research Methods 0.000 abstract description 25
- 239000000126 substance Substances 0.000 abstract description 13
- 230000000052 comparative effect Effects 0.000 abstract description 4
- 235000013350 formula milk Nutrition 0.000 description 70
- 241000700159 Rattus Species 0.000 description 52
- 239000000306 component Substances 0.000 description 47
- 241000252212 Danio rerio Species 0.000 description 28
- 230000033115 angiogenesis Effects 0.000 description 20
- 238000002474 experimental method Methods 0.000 description 19
- 230000002107 myocardial effect Effects 0.000 description 19
- 206010019280 Heart failures Diseases 0.000 description 15
- 244000131316 Panax pseudoginseng Species 0.000 description 15
- 230000037149 energy metabolism Effects 0.000 description 14
- 230000008569 process Effects 0.000 description 14
- 210000001519 tissue Anatomy 0.000 description 13
- 239000000463 material Substances 0.000 description 12
- 230000025102 vascular smooth muscle contraction Effects 0.000 description 12
- 210000004369 blood Anatomy 0.000 description 11
- 239000008280 blood Substances 0.000 description 11
- 210000004204 blood vessel Anatomy 0.000 description 11
- 229940079593 drug Drugs 0.000 description 11
- 238000011529 RT qPCR Methods 0.000 description 10
- 230000015572 biosynthetic process Effects 0.000 description 10
- 239000003270 steroid hormone Substances 0.000 description 10
- IAZDPXIOMUYVGZ-UHFFFAOYSA-N Dimethylsulphoxide Chemical compound CS(C)=O IAZDPXIOMUYVGZ-UHFFFAOYSA-N 0.000 description 9
- 208000010125 myocardial infarction Diseases 0.000 description 9
- 239000000843 powder Substances 0.000 description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 9
- 210000004165 myocardium Anatomy 0.000 description 8
- 238000003786 synthesis reaction Methods 0.000 description 8
- 241001465754 Metazoa Species 0.000 description 7
- 230000017531 blood circulation Effects 0.000 description 7
- 210000000936 intestine Anatomy 0.000 description 7
- 230000003285 pharmacodynamic effect Effects 0.000 description 7
- 238000004088 simulation Methods 0.000 description 7
- 238000012360 testing method Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 230000000875 corresponding effect Effects 0.000 description 6
- 238000001514 detection method Methods 0.000 description 6
- 230000008014 freezing Effects 0.000 description 6
- 238000007710 freezing Methods 0.000 description 6
- 230000000144 pharmacologic effect Effects 0.000 description 6
- 238000002360 preparation method Methods 0.000 description 6
- LTMHDMANZUZIPE-AMTYYWEZSA-N Digoxin Natural products O([C@H]1[C@H](C)O[C@H](O[C@@H]2C[C@@H]3[C@@](C)([C@@H]4[C@H]([C@]5(O)[C@](C)([C@H](O)C4)[C@H](C4=CC(=O)OC4)CC5)CC3)CC2)C[C@@H]1O)[C@H]1O[C@H](C)[C@@H](O[C@H]2O[C@@H](C)[C@H](O)[C@@H](O)C2)[C@@H](O)C1 LTMHDMANZUZIPE-AMTYYWEZSA-N 0.000 description 5
- 241000304195 Salvia miltiorrhiza Species 0.000 description 5
- 235000011135 Salvia miltiorrhiza Nutrition 0.000 description 5
- 239000011575 calcium Substances 0.000 description 5
- 208000029078 coronary artery disease Diseases 0.000 description 5
- LTMHDMANZUZIPE-PUGKRICDSA-N digoxin Chemical compound C1[C@H](O)[C@H](O)[C@@H](C)O[C@H]1O[C@@H]1[C@@H](C)O[C@@H](O[C@@H]2[C@H](O[C@@H](O[C@@H]3C[C@@H]4[C@]([C@@H]5[C@H]([C@]6(CC[C@@H]([C@@]6(C)[C@H](O)C5)C=5COC(=O)C=5)O)CC4)(C)CC3)C[C@@H]2O)C)C[C@@H]1O LTMHDMANZUZIPE-PUGKRICDSA-N 0.000 description 5
- 229960005156 digoxin Drugs 0.000 description 5
- LTMHDMANZUZIPE-UHFFFAOYSA-N digoxine Natural products C1C(O)C(O)C(C)OC1OC1C(C)OC(OC2C(OC(OC3CC4C(C5C(C6(CCC(C6(C)C(O)C5)C=5COC(=O)C=5)O)CC4)(C)CC3)CC2O)C)CC1O LTMHDMANZUZIPE-UHFFFAOYSA-N 0.000 description 5
- 230000000857 drug effect Effects 0.000 description 5
- 235000008216 herbs Nutrition 0.000 description 5
- 239000007788 liquid Substances 0.000 description 5
- 239000012452 mother liquor Substances 0.000 description 5
- 208000031225 myocardial ischemia Diseases 0.000 description 5
- 238000011084 recovery Methods 0.000 description 5
- 150000003431 steroids Chemical class 0.000 description 5
- 210000003462 vein Anatomy 0.000 description 5
- 238000012795 verification Methods 0.000 description 5
- 241000251468 Actinopterygii Species 0.000 description 4
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 4
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 4
- BHPQYMZQTOCNFJ-UHFFFAOYSA-N Calcium cation Chemical compound [Ca+2] BHPQYMZQTOCNFJ-UHFFFAOYSA-N 0.000 description 4
- 102100031181 Glyceraldehyde-3-phosphate dehydrogenase Human genes 0.000 description 4
- 101001012157 Homo sapiens Receptor tyrosine-protein kinase erbB-2 Proteins 0.000 description 4
- 206010061216 Infarction Diseases 0.000 description 4
- 241001098657 Pterois Species 0.000 description 4
- 102100030086 Receptor tyrosine-protein kinase erbB-2 Human genes 0.000 description 4
- RYMZZMVNJRMUDD-UHFFFAOYSA-N SJ000286063 Natural products C12C(OC(=O)C(C)(C)CC)CC(C)C=C2C=CC(C)C1CCC1CC(O)CC(=O)O1 RYMZZMVNJRMUDD-UHFFFAOYSA-N 0.000 description 4
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 4
- 241000282898 Sus scrofa Species 0.000 description 4
- 230000004071 biological effect Effects 0.000 description 4
- 239000012503 blood component Substances 0.000 description 4
- 229910052791 calcium Inorganic materials 0.000 description 4
- 229910001424 calcium ion Inorganic materials 0.000 description 4
- 230000000747 cardiac effect Effects 0.000 description 4
- RNFNDJAIBTYOQL-UHFFFAOYSA-N chloral hydrate Chemical compound OC(O)C(Cl)(Cl)Cl RNFNDJAIBTYOQL-UHFFFAOYSA-N 0.000 description 4
- 229960002327 chloral hydrate Drugs 0.000 description 4
- 239000002299 complementary DNA Substances 0.000 description 4
- 210000002808 connective tissue Anatomy 0.000 description 4
- 210000004351 coronary vessel Anatomy 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 108020004445 glyceraldehyde-3-phosphate dehydrogenase Proteins 0.000 description 4
- 230000007574 infarction Effects 0.000 description 4
- NOESYZHRGYRDHS-UHFFFAOYSA-N insulin Chemical group N1C(=O)C(NC(=O)C(CCC(N)=O)NC(=O)C(CCC(O)=O)NC(=O)C(C(C)C)NC(=O)C(NC(=O)CN)C(C)CC)CSSCC(C(NC(CO)C(=O)NC(CC(C)C)C(=O)NC(CC=2C=CC(O)=CC=2)C(=O)NC(CCC(N)=O)C(=O)NC(CC(C)C)C(=O)NC(CCC(O)=O)C(=O)NC(CC(N)=O)C(=O)NC(CC=2C=CC(O)=CC=2)C(=O)NC(CSSCC(NC(=O)C(C(C)C)NC(=O)C(CC(C)C)NC(=O)C(CC=2C=CC(O)=CC=2)NC(=O)C(CC(C)C)NC(=O)C(C)NC(=O)C(CCC(O)=O)NC(=O)C(C(C)C)NC(=O)C(CC(C)C)NC(=O)C(CC=2NC=NC=2)NC(=O)C(CO)NC(=O)CNC2=O)C(=O)NCC(=O)NC(CCC(O)=O)C(=O)NC(CCCNC(N)=N)C(=O)NCC(=O)NC(CC=3C=CC=CC=3)C(=O)NC(CC=3C=CC=CC=3)C(=O)NC(CC=3C=CC(O)=CC=3)C(=O)NC(C(C)O)C(=O)N3C(CCC3)C(=O)NC(CCCCN)C(=O)NC(C)C(O)=O)C(=O)NC(CC(N)=O)C(O)=O)=O)NC(=O)C(C(C)CC)NC(=O)C(CO)NC(=O)C(C(C)O)NC(=O)C1CSSCC2NC(=O)C(CC(C)C)NC(=O)C(NC(=O)C(CCC(N)=O)NC(=O)C(CC(N)=O)NC(=O)C(NC(=O)C(N)CC=1C=CC=CC=1)C(C)C)CC1=CN=CN1 NOESYZHRGYRDHS-UHFFFAOYSA-N 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 4
- 210000003205 muscle Anatomy 0.000 description 4
- 238000001558 permutation test Methods 0.000 description 4
- 239000002504 physiological saline solution Substances 0.000 description 4
- 230000001737 promoting effect Effects 0.000 description 4
- 238000010839 reverse transcription Methods 0.000 description 4
- 210000005241 right ventricle Anatomy 0.000 description 4
- RYMZZMVNJRMUDD-HGQWONQESA-N simvastatin Chemical compound C([C@H]1[C@@H](C)C=CC2=C[C@H](C)C[C@@H]([C@H]12)OC(=O)C(C)(C)CC)C[C@@H]1C[C@@H](O)CC(=O)O1 RYMZZMVNJRMUDD-HGQWONQESA-N 0.000 description 4
- 229960002855 simvastatin Drugs 0.000 description 4
- 235000014347 soups Nutrition 0.000 description 4
- 239000006228 supernatant Substances 0.000 description 4
- 210000000115 thoracic cavity Anatomy 0.000 description 4
- 230000002792 vascular Effects 0.000 description 4
- 238000005303 weighing Methods 0.000 description 4
- WEVYAHXRMPXWCK-UHFFFAOYSA-N Acetonitrile Chemical compound CC#N WEVYAHXRMPXWCK-UHFFFAOYSA-N 0.000 description 3
- 206010002091 Anaesthesia Diseases 0.000 description 3
- SMDOOINVMJSDPS-UHFFFAOYSA-N Astragaloside Natural products C1=C(O)C(OC)=CC(C2=C(C(=O)C3=C(O)C=C(O)C=C3O2)OC2C(C(OC3C(C(O)C(O)C(CO)O3)O)C(O)C(CO)O2)O)=C1 SMDOOINVMJSDPS-UHFFFAOYSA-N 0.000 description 3
- 241000202726 Bupleurum Species 0.000 description 3
- 101100227322 Caenorhabditis elegans fli-1 gene Proteins 0.000 description 3
- 244000037364 Cinnamomum aromaticum Species 0.000 description 3
- 235000014489 Cinnamomum aromaticum Nutrition 0.000 description 3
- 241000380450 Danaus melanippus Species 0.000 description 3
- OKKJLVBELUTLKV-UHFFFAOYSA-N Methanol Chemical compound OC OKKJLVBELUTLKV-UHFFFAOYSA-N 0.000 description 3
- 239000012901 Milli-Q water Substances 0.000 description 3
- 101100281205 Mus musculus Fli1 gene Proteins 0.000 description 3
- 244000236658 Paeonia lactiflora Species 0.000 description 3
- 235000008598 Paeonia lactiflora Nutrition 0.000 description 3
- 235000003181 Panax pseudoginseng Nutrition 0.000 description 3
- 244000197580 Poria cocos Species 0.000 description 3
- 235000008599 Poria cocos Nutrition 0.000 description 3
- 241000282894 Sus scrofa domesticus Species 0.000 description 3
- 241001455273 Tetrapoda Species 0.000 description 3
- DOQPXTMNIUCOSY-UHFFFAOYSA-N [4-cyano-4-(3,4-dimethoxyphenyl)-5-methylhexyl]-[2-(3,4-dimethoxyphenyl)ethyl]-methylazanium;chloride Chemical compound [H+].[Cl-].C1=C(OC)C(OC)=CC=C1CCN(C)CCCC(C#N)(C(C)C)C1=CC=C(OC)C(OC)=C1 DOQPXTMNIUCOSY-UHFFFAOYSA-N 0.000 description 3
- 230000037005 anaesthesia Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 3
- QMNWISYXSJWHRY-XWJCTJPOSA-N astragaloside Chemical compound O1[C@H](C(C)(O)C)CC[C@]1(C)[C@@H]1[C@@]2(C)CC[C@]34C[C@]4(CC[C@H](O[C@H]4[C@@H]([C@@H](O)[C@H](O)CO4)O)C4(C)C)C4[C@@H](O[C@H]4[C@@H]([C@@H](O)[C@H](O)[C@@H](CO)O4)O)CC3[C@]2(C)C[C@@H]1O QMNWISYXSJWHRY-XWJCTJPOSA-N 0.000 description 3
- 230000037396 body weight Effects 0.000 description 3
- 230000028956 calcium-mediated signaling Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 3
- 238000010828 elution Methods 0.000 description 3
- 230000036541 health Effects 0.000 description 3
- 230000004155 insulin signaling pathway Effects 0.000 description 3
- 238000012353 t test Methods 0.000 description 3
- 238000013518 transcription Methods 0.000 description 3
- 230000035897 transcription Effects 0.000 description 3
- 229960000881 verapamil hydrochloride Drugs 0.000 description 3
- 239000012224 working solution Substances 0.000 description 3
- 101150068220 CACNB2 gene Proteins 0.000 description 2
- 241000721047 Danaus plexippus Species 0.000 description 2
- RTZKZFJDLAIYFH-UHFFFAOYSA-N Diethyl ether Chemical compound CCOCC RTZKZFJDLAIYFH-UHFFFAOYSA-N 0.000 description 2
- 102000004877 Insulin Human genes 0.000 description 2
- 108090001061 Insulin Proteins 0.000 description 2
- 206010022489 Insulin Resistance Diseases 0.000 description 2
- 101150033052 MAS5 gene Proteins 0.000 description 2
- 238000009004 PCR Kit Methods 0.000 description 2
- 229930182555 Penicillin Natural products 0.000 description 2
- JGSARLDLIJGVTE-MBNYWOFBSA-N Penicillin G Chemical group N([C@H]1[C@H]2SC([C@@H](N2C1=O)C(O)=O)(C)C)C(=O)CC1=CC=CC=C1 JGSARLDLIJGVTE-MBNYWOFBSA-N 0.000 description 2
- 238000010802 RNA extraction kit Methods 0.000 description 2
- 239000013614 RNA sample Substances 0.000 description 2
- 101100344462 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) YDJ1 gene Proteins 0.000 description 2
- 101150081767 TXN2 gene Proteins 0.000 description 2
- 102000004987 Troponin T Human genes 0.000 description 2
- 108090001108 Troponin T Proteins 0.000 description 2
- 210000000683 abdominal cavity Anatomy 0.000 description 2
- 230000004913 activation Effects 0.000 description 2
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 2
- 229910052782 aluminium Inorganic materials 0.000 description 2
- 210000000702 aorta abdominal Anatomy 0.000 description 2
- 230000006793 arrhythmia Effects 0.000 description 2
- 206010003119 arrhythmia Diseases 0.000 description 2
- 238000003705 background correction Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000000740 bleeding effect Effects 0.000 description 2
- 238000009395 breeding Methods 0.000 description 2
- 230000001488 breeding effect Effects 0.000 description 2
- 230000024683 calcium ion homeostasis Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 210000000038 chest Anatomy 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 230000008828 contractile function Effects 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 2
- 239000013078 crystal Substances 0.000 description 2
- ZOOGRGPOEVQQDX-KHLHZJAASA-N cyclic guanosine monophosphate Chemical compound C([C@H]1O2)O[P@](O)(=O)O[C@@H]1[C@H](O)[C@H]2N1C(N=C(NC2=O)N)=C2N=C1 ZOOGRGPOEVQQDX-KHLHZJAASA-N 0.000 description 2
- 238000012899 de-mixing Methods 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 210000002889 endothelial cell Anatomy 0.000 description 2
- 230000004720 fertilization Effects 0.000 description 2
- 238000005206 flow analysis Methods 0.000 description 2
- 239000011888 foil Substances 0.000 description 2
- 235000011389 fruit/vegetable juice Nutrition 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 230000002496 gastric effect Effects 0.000 description 2
- 238000003304 gavage Methods 0.000 description 2
- 229930182494 ginsenoside Natural products 0.000 description 2
- 210000005003 heart tissue Anatomy 0.000 description 2
- 238000009396 hybridization Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 229940125396 insulin Drugs 0.000 description 2
- 239000007928 intraperitoneal injection Substances 0.000 description 2
- 229960004393 lidocaine hydrochloride Drugs 0.000 description 2
- YECIFGHRMFEPJK-UHFFFAOYSA-N lidocaine hydrochloride monohydrate Chemical compound O.[Cl-].CC[NH+](CC)CC(=O)NC1=C(C)C=CC=C1C YECIFGHRMFEPJK-UHFFFAOYSA-N 0.000 description 2
- 238000004895 liquid chromatography mass spectrometry Methods 0.000 description 2
- 230000013011 mating Effects 0.000 description 2
- 208000030159 metabolic disease Diseases 0.000 description 2
- 230000004060 metabolic process Effects 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 231100000062 no-observed-adverse-effect level Toxicity 0.000 description 2
- 238000004806 packaging method and process Methods 0.000 description 2
- 230000036407 pain Effects 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 229940049954 penicillin Drugs 0.000 description 2
- 238000003752 polymerase chain reaction Methods 0.000 description 2
- 239000013641 positive control Substances 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 210000002966 serum Anatomy 0.000 description 2
- 210000000329 smooth muscle myocyte Anatomy 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000007447 staining method Methods 0.000 description 2
- 239000011550 stock solution Substances 0.000 description 2
- 208000011580 syndromic disease Diseases 0.000 description 2
- 230000009261 transgenic effect Effects 0.000 description 2
- 208000001072 type 2 diabetes mellitus Diseases 0.000 description 2
- 229910021642 ultra pure water Inorganic materials 0.000 description 2
- 239000012498 ultrapure water Substances 0.000 description 2
- SGTNSNPWRIOYBX-UHFFFAOYSA-N 2-(3,4-dimethoxyphenyl)-5-{[2-(3,4-dimethoxyphenyl)ethyl](methyl)amino}-2-(propan-2-yl)pentanenitrile Chemical compound C1=C(OC)C(OC)=CC=C1CCN(C)CCCC(C#N)(C(C)C)C1=CC=C(OC)C(OC)=C1 SGTNSNPWRIOYBX-UHFFFAOYSA-N 0.000 description 1
- 241000282465 Canis Species 0.000 description 1
- 206010052337 Diastolic dysfunction Diseases 0.000 description 1
- 102000001301 EGF receptor Human genes 0.000 description 1
- 108060006698 EGF receptor Proteins 0.000 description 1
- 206010048554 Endothelial dysfunction Diseases 0.000 description 1
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 description 1
- 108010078321 Guanylate Cyclase Proteins 0.000 description 1
- 102000014469 Guanylate cyclase Human genes 0.000 description 1
- 239000009636 Huang Qi Substances 0.000 description 1
- 206010020565 Hyperaemia Diseases 0.000 description 1
- 206010021143 Hypoxia Diseases 0.000 description 1
- 206010061218 Inflammation Diseases 0.000 description 1
- 238000012351 Integrated analysis Methods 0.000 description 1
- 102000048238 Neuregulin-1 Human genes 0.000 description 1
- 108090000556 Neuregulin-1 Proteins 0.000 description 1
- 235000006484 Paeonia officinalis Nutrition 0.000 description 1
- 241001106477 Paeoniaceae Species 0.000 description 1
- 241000180649 Panax notoginseng Species 0.000 description 1
- 235000003143 Panax notoginseng Nutrition 0.000 description 1
- 208000037273 Pathologic Processes Diseases 0.000 description 1
- 235000001188 Peltandra virginica Nutrition 0.000 description 1
- 241001619461 Poria <basidiomycete fungus> Species 0.000 description 1
- 206010063837 Reperfusion injury Diseases 0.000 description 1
- YMGFTDKNIWPMGF-QHCPKHFHSA-N Salvianolic acid A Natural products OC(=O)[C@H](Cc1ccc(O)c(O)c1)OC(=O)C=Cc2ccc(O)c(O)c2C=Cc3ccc(O)c(O)c3 YMGFTDKNIWPMGF-QHCPKHFHSA-N 0.000 description 1
- YMGFTDKNIWPMGF-UCPJVGPRSA-N Salvianolic acid A Chemical compound C([C@H](C(=O)O)OC(=O)\C=C\C=1C(=C(O)C(O)=CC=1)\C=C\C=1C=C(O)C(O)=CC=1)C1=CC=C(O)C(O)=C1 YMGFTDKNIWPMGF-UCPJVGPRSA-N 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 208000005392 Spasm Diseases 0.000 description 1
- 102000005789 Vascular Endothelial Growth Factors Human genes 0.000 description 1
- 108010019530 Vascular Endothelial Growth Factors Proteins 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- XBJFCYDKBDVADW-UHFFFAOYSA-N acetonitrile;formic acid Chemical compound CC#N.OC=O XBJFCYDKBDVADW-UHFFFAOYSA-N 0.000 description 1
- 239000002671 adjuvant Substances 0.000 description 1
- 230000036592 analgesia Effects 0.000 description 1
- 238000010171 animal model Methods 0.000 description 1
- 230000003527 anti-angiogenesis Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- QMNWISYXSJWHRY-YLNUDOOFSA-N astragaloside IV Chemical compound O1[C@H](C(C)(O)C)CC[C@]1(C)[C@@H]1[C@@]2(C)CC[C@]34C[C@]4(CC[C@H](O[C@H]4[C@@H]([C@@H](O)[C@H](O)CO4)O)C4(C)C)[C@H]4[C@@H](O[C@H]4[C@@H]([C@@H](O)[C@H](O)[C@@H](CO)O4)O)C[C@H]3[C@]2(C)C[C@@H]1O QMNWISYXSJWHRY-YLNUDOOFSA-N 0.000 description 1
- QMNWISYXSJWHRY-BCBPIKMJSA-N astragaloside IV Natural products CC(C)(O)[C@@H]1CC[C@@](C)(O1)[C@H]2[C@@H](O)C[C@@]3(C)[C@@H]4C[C@H](O[C@@H]5O[C@H](CO)[C@H](O)[C@@H](O)[C@H]5O)[C@H]6C(C)(C)[C@H](CC[C@@]67C[C@@]47CC[C@]23C)O[C@@H]8OC[C@@H](O)[C@H](O)[C@H]8O QMNWISYXSJWHRY-BCBPIKMJSA-N 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 210000000748 cardiovascular system Anatomy 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 230000004663 cell proliferation Effects 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 229940126678 chinese medicines Drugs 0.000 description 1
- 238000013375 chromatographic separation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- PFKIBRPYVNVMRU-UHFFFAOYSA-N cyclosieversioside F Natural products CC(C)(O)C1COC(C)(C1)C2C(O)CC3(C)C4CC(OC5OC(CO)C(O)C(O)C5O)C6C(C)(C)C(CCC67CC47CCC23C)OC8OCC(O)C(O)C8O PFKIBRPYVNVMRU-UHFFFAOYSA-N 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000002074 deregulated effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 235000014113 dietary fatty acids Nutrition 0.000 description 1
- 238000006471 dimerization reaction Methods 0.000 description 1
- 208000035475 disorder Diseases 0.000 description 1
- 238000004090 dissolution Methods 0.000 description 1
- 230000002222 downregulating effect Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000000132 electrospray ionisation Methods 0.000 description 1
- 230000008694 endothelial dysfunction Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 229930195729 fatty acid Natural products 0.000 description 1
- 239000000194 fatty acid Substances 0.000 description 1
- 150000004665 fatty acids Chemical class 0.000 description 1
- HQVFCQRVQFYGRJ-UHFFFAOYSA-N formic acid;hydrate Chemical compound O.OC=O HQVFCQRVQFYGRJ-UHFFFAOYSA-N 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000003209 gene knockout Methods 0.000 description 1
- 229940089161 ginsenoside Drugs 0.000 description 1
- 239000008103 glucose Substances 0.000 description 1
- 230000012010 growth Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000004128 high performance liquid chromatography Methods 0.000 description 1
- 230000007954 hypoxia Effects 0.000 description 1
- 238000011534 incubation Methods 0.000 description 1
- 239000000411 inducer Substances 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 230000004054 inflammatory process Effects 0.000 description 1
- 230000003834 intracellular effect Effects 0.000 description 1
- 238000005040 ion trap Methods 0.000 description 1
- 150000002500 ions Chemical class 0.000 description 1
- 208000028867 ischemia Diseases 0.000 description 1
- 230000000302 ischemic effect Effects 0.000 description 1
- 238000004811 liquid chromatography Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004949 mass spectrometry Methods 0.000 description 1
- 210000003584 mesangial cell Anatomy 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000010413 mother solution Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000036284 oxygen consumption Effects 0.000 description 1
- 230000009054 pathological process Effects 0.000 description 1
- 230000035755 proliferation Effects 0.000 description 1
- 230000005588 protonation Effects 0.000 description 1
- 239000008412 qishen yiqi Substances 0.000 description 1
- 230000007115 recruitment Effects 0.000 description 1
- 229930183842 salvianolic acid Natural products 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000009877 shengmai Substances 0.000 description 1
- 210000002460 smooth muscle Anatomy 0.000 description 1
- 238000010186 staining Methods 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 230000002195 synergetic effect Effects 0.000 description 1
- 231100000816 toxic dose Toxicity 0.000 description 1
- 230000009724 venous congestion Effects 0.000 description 1
- 229960001722 verapamil Drugs 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B50/00—ICT programming tools or database systems specially adapted for bioinformatics
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Biotechnology (AREA)
- Evolutionary Biology (AREA)
- Biophysics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Bioethics (AREA)
- Databases & Information Systems (AREA)
- Genetics & Genomics (AREA)
- Molecular Biology (AREA)
- Medicines Containing Plant Substances (AREA)
Abstract
本发明公开了一种中药对复杂性疾病所致“失和”网络整合调节作用的定量评价方法,该方法包括:(1)获取中药对复杂性疾病的潜在作用靶点,利用所述潜在作用靶点构建所述“失和”生物分子网络;(2)通过综合考虑中药中入血成分与靶蛋白的对接得分值以及成分在中药中的相对含量,计算获得中药对复杂性疾病所致的“失和”生物分子网络的调控力度值;(3)比较不同中药之间的调控力度值,评价各中药对所述“失和”生物分子网络的整合调节作用。本发明方法通过辨识中药化学物质基础与复杂性疾病所致的“失和”生物分子网络的相关性,开展中药之间内在异同的比较研究,探索诠释中药的配伍科学内涵,丰富和拓展了中药配伍规律的研究思路。
Description
技术领域
本发明涉及中药药理学领域,尤其涉及一种中药对复杂性疾病所致“失和”网络整合调节作用的定量评价方法。
背景技术
中医药是中华民族数千年来在生活实践和与疾病斗争的过程中积累的瑰宝,是在几千年临床实践中逐步发展而成的医学科学,直至今日仍在国民健康事业中具有举足轻重的作用,同时也为全人类的健康做出了重要贡献。
现代研究表明,中药中多种药效成分通过协调互济作用从整体水平调节机体状态,具有“多成分、多靶点、多途径”的整合调节作用,在治疗复杂性疾病方面具有不可替代的优势。
中药通过其化学物质组成与机体生命系统的相互作用发挥药效作用,因此,中药化学物质成分辨识、药效药理作用研究以及成分与生物效应之间的相关性辨识等成为了研究中药配伍科学内涵的关键环节。
目前,国内外对中药配伍规律的研究主要集中在其具体内涵的理论研究和应用现代实验手段揭示其配伍原理等方面。其中,文献研究着重于文献的整理和知识挖掘,主要有基于频次统计和关联规则、基于复杂系统熵聚类、基于人工神经网络等方法;而实验研究主要以中药化学成分辨识、中药生物学效应或化学成分与生物效应之间的交互作用为切入点来探讨中药的配伍规律及其科学内涵。
例如:(1)丁选胜等(丁选胜,安娜芹,姚丹,赵一鸣.人参白虎汤不同配伍条件下人参皂苷Rg1和Re变化的研究.中草药2007,No.421(05):702-704.)通过HPLC法测定不同配比的人参白虎汤中人参皂苷Re和Rg1的含量变化情况,探讨人参白虎汤的配伍规律,结果发现人参配伍方剂中其他药材后,一定程度增加了Re和Rg1的溶出量,证明了人参白虎汤配伍比例的合理性。
(2)商洪才等(商洪才,张伯礼,高秀梅,王怡,李玉红,康立源,杜嵘,江海涛.丹参三七药对不同配比药效学比较研究.辽宁中医杂志2002(05):297-299.)采用冠状动脉左前降支结扎法构建犬急性心肌缺血模型,观察不同比例丹参-三七配伍后的药效,结果表明丹参、三七配伍后药效作用强于各单味药材的药效,且以10:6和10:3配伍的丹参-三七药效最突出。
(3)刘慧兰等(刘慧兰,欧阳建军.桂枝、柴胡与白芍分别配伍的相关药效学研究.湖南中医药大学学报2007,No.127(03):31-33.)研究了桂枝、柴胡与白芍分别配伍后的作用,结果表明柴胡和白芍配伍后具有解痉和镇痛作用;桂枝和白芍配伍后具有抗炎和镇痛作用,且在镇痛方面具有协同效应。
(4)刘明平等(刘明平,林铿,黄兆胜,吴庆光.通脉复方有效组分抑制肾小球系膜细胞增殖的活性筛选及配伍研究.中国实验方剂学杂志2011,v.17(08):172-175.)通过考察通脉复方(黄芪和三七)组成中有效组分对肾小球系膜细胞增殖的影响,发现总多糖和总皂苷是通脉复方中的主要有效成分,且4:1为最佳的剂量配比。
(5)吕琳星等(吕琳星,范雪梅,梁琼麟,王文凯,王义明,罗国安.基因芯片用于组分中药新双龙方的配伍机制研究.高等学校化学学报2012,v.33(11):2397-2404.)采用基因芯片技术研究新双龙方剂及其有效组分(人参总皂苷和丹参总酚酸)治疗急性心肌梗塞大鼠的作用机制,研究结果发现新双龙方组能比有效组分组更积极的调控相关的靶基因,且新双龙方能富集到更多的通路,从而证实了新双龙方的配伍优势。
随着网络药理学的兴起,基于中药网络药理学的方剂配伍规律研究也拉开了帷幕,并取得了一定的进展。
吴磊宏等(Wu L,Wang Y,Li Z,Zhang B,Cheng Y,Fan X.Identifying roles of“Jun-Chen-Zuo-Shi”component herbs of QiShenYiQi formula in treating acutemyocardial ischemia by network pharmacology.Chinese medicine 2014,9(1):1-9.)通过基于网络药理学的方法探讨芪参益气方治疗急性心肌缺血的药理作用,从分子网络和通路水平解释了方剂中药材“君臣佐使”的配伍原则。任钧国等(任钧国,马晓斌,林成仁,李鸿海,王敏,李军梅,王杨慧,刘建勋.基于网络方法的中药方剂药理作用配伍规律分析.中国中药杂志2010,v.35(18):2469-2474.)以加味生脉散为研究对象,构建基于网络的中药方剂药理作用配伍网络,研究了其抗大鼠急性心肌缺血-再灌注损伤作用的配伍规律。Li等(Li X,Xu X,Wang J,Yu H,Wang X,Yang H,Xu H,Tang S,Li Y,Yang L.A system-levelinvestigation into the mechanisms of Chinese Traditional Medicine:CompoundDanshen Formula for cardiovascular disease treatment.PLoS One 2012,7(9):e43918.)提出了基于化学成分、化学基因组学和药理学数据的网络方法系统的探究了复方丹参方治疗心血管疾病的药理作用机制,并阐明了方剂中丹参是君药、三七是臣药、冰片是佐药的配伍原则。
虽然已有大量关于中药配伍规律的研究,并在一批中药的配伍规律研究上取得了良好进展,但相关研究多以还原论思维模式开展,从单一的药效指标或物质基础角度分析中药配伍规律,这与中药的整体性原则存在分歧,且往往与病证脱节、缺乏传统理论依据,与临床实际应用仍有距离。导致这一局面的主要原因是,中药配伍规律研究仍存在以下瓶颈问题尚未突破:1)中药数量庞大且绝大部分由古人所留下,其配伍规律及科学内涵的描述和记载与现代医学体系分属两套不同的描述符号系统,难以被理解与接受。2)大部分的中药,尤其是方剂,均由两味以上的药材组成,化学成分复杂,且成分间存在相互作用,加大了中药配伍规律研究的难度。3)中药通过作用于机体的特定证候而发挥作用,而生命机体高度复杂,且中医证候现代化研究相对匮乏,严重阻碍了中药配伍规律及其科学内涵的深入研究。
因此,面对如此庞大的中药数量,面临着中药化学组成和机体生命系统的双重复杂性,亟需拓宽思路,引入现代技术手段对宏观理念和微观事实进行有机衔接,以期揭示中药配伍的科学内涵。
发明内容
本发明提供了一种中药对复杂性疾病所致“失和”网络整合调节作用的定量评价方法,该方法结合中药中入血成分与靶蛋白的对接得分值和成分在中药中的相对含量,有效辨识中药化学组成与机体生物分子网络之间的相关性,定量描述中药对其相关“失和”生物分子网络的调控力度(即中药对复杂性疾病所致“失和”网络的整合调节作用)。
一种中药对复杂性疾病所致“失和”网络整合调节作用的定量评价方法,包括以下步骤:
(1)获取中药对复杂性疾病的潜在作用靶点,利用所述潜在作用靶点构建所述“失和”生物分子网络;
(2)计算获得中药对所述“失和”生物分子网络的调控力度值;
a)利用所述中药中入血成分与所述“失和”生物分子网络中靶蛋白的对接得分值,计算所述“失和”生物分子网络中各连接的权重;
式(1)中,Scomp表示中药中入血成分与网络中蛋白进行分子对接的得分值;Sreference表示网络中蛋白质复合结构的原始配体与蛋白进行分子对接的最低得分值;EW表示网络中的任何一个节点到下游直接节点的连接的权重;
b)将所述权重代入式(2)和式(3)中,计算获得网络效能和网络流;
式(2)中,dwij为网络中节点i与节点j之间的最短距离,即节点i到节点j之间所经过的边的权重之和最小;N为网络中所有节点的集合,NE为网络效能,即网络中任何一个节点i到下游节点j的最短路径的倒数和;
式(3)中,dwim为网络中节点i与节点之间的最短距离,即节点i到节点m之间所经过的边的权重之和最小;N为网络中所有节点的集合,NF为网络流,即网络中任何一个节点i到网络中叶节点m的最短路径的倒数和;LeafNode表示叶节点,是指网络中出度为0的节点;
c)将所述网络效能和网络流代入式(4)中,计算获得中药中入血成分对所述“失和”生物分子网络的调控能力;
式(4)中,当所述“失和”生物分子网络中所有连接的权重均为初始值1时,计算得到NEmax和NFmax;NEcomp表示中药中单个入血成分对所述“失和”生物分子网络的网络效能,NFcomp表示中药中单个入血成分对所述“失和”生物分子网络的网络流;RE表示中药中入血成分对所述“失和”生物分子网络的调控能力;
d)将所述“失和”生物分子网络的调控能力代入式(5)中,计算得到中药对复杂性疾病所致“失和”生物分子网络的调控力度值;
式(5)中,k表示中药中的入血成分,G表示中药中所有入血成分的集合,Ak表示成分k在中药所有入血成分中所占的峰面积百分比,FRE表示中药对复杂性疾病所致“失和”生物分子网络的调控力度值;
(3)比较不同中药之间的调控力度值,评价各中药对所述“失和”生物分子网络的调控能力。
本发明方法中的分子对接是应用Molecular Operating Environment(MOE,v2014.09)的分子对接模块完成蛋白质和化合物的分子对接工作;其主要流程包括蛋白结构矫正(Structure Preparation)、质子化(Protonate3D)、配体构象优化和应用打分函数进行评分等。
本发明结合入血成分在LC-MS分析中的相对含量设计了用于评价中药对网络调控能力的方法。网络中的任何一个节点(i)到下游直接节点(j)的连接权重根据公式(1)计算得到;设定网络中所有连接的初始权重为1,最大权重为200。每一个网络中蛋白质复合结构中的原始配体与该蛋白进行分子对接模拟,选取最低得分值作为该网络的Sreference,Scomp为化合物与网络中蛋白进行分子对接的得分值。根据网络中连接权重的不同,按照公式(2)计算网络效能,其中dwij为网络中任意两节点间具有最小值的最短距离,N为网络中所有节点的集合。在有向网络中,下游蛋白对网络功能具有更重要的调控作用,因此引入网络流权衡网络中不同节点的重要程度。
进一步地,所述中药为四逆汤类方。
更优选,所述中药为四逆汤、四逆加人参汤、茯苓四逆汤、通脉四逆汤或通脉四逆加猪胆汁汤。
进一步地,所述复杂性疾病为心脑血管疾病。
更优选,所述复杂性疾病为急性心肌梗塞。
进一步地,步骤(1)中,获取所述潜在作用靶点的方法,包括以下步骤:
A)获得中药作用于复杂性疾病的芯片数据;
B)采用cMAP法得到中药作用复杂性疾病的潜在相关靶蛋白;
C)采用EoR法得到中药作用复杂性疾病的潜在相关基因;
D)整合所述潜在相关靶蛋白和潜在相关基因,得到中药作用复杂性疾病的潜在作用靶点。
进一步地,所述cMAP法,包括以下步骤:
(a)计算探针的FC值,当探针的FC绝对值大于1.2,且p<0.05时,定义为显著差异表达的探针;
(b)根据探针改变的幅度大小,挑选中药芯片数据中上调或下调幅度最大的显著差异表达的探针与cMAP数据库中的芯片数据进行相似度匹配,找到与中药相关性最强的化合物;
(c)根据STITCH4.0数据库中化合物-蛋白质相互作用的数据信息,得到潜在作用蛋白靶点。
FC(Fold Change)值,即FC(SNs/Model)值,是指中药给药组(SNs)与模型对照组(Model)基因转录水平表达强度的比值。
进一步地,所述EoR法,包括以下步骤:
A)根据芯片数据计算节点的回调状态;
式(6)中,RL′表示回调状态;LogE为各组的基因表达水平取log2对数后的均值;Edrug为中药给药组的芯片数据;Esham为假手术组的芯片数据;Emodel为模型对照组的芯片数据。
B)计算节点的回复调控效率
EoR=100%-|100%-RL′| (7)
式(7)中,EoR表示为回复调控效率,是用于表征节点回调效率的指标;RL′表示回调状态;EoR值大于0的节点为所述潜在相关基因。
进一步地,所述“失和”生物分子网络的构建,包括以下步骤:
A)应用DAVID Bioinformatics Resources 6.7软件对靶点列表进行通路富集分析,找到所述复杂性疾病相关的生物学信号通路,定义p<0.05为显著富集通路;
B)以KEGG数据库中显著富集的通路为基础,构建所述“失和”生物分子网络。
进一步地,所述通路为Calcium signaling pathway(钙离子信号通路)、ErbBsignaling pathway(ErbB信号通路)、Vascular smooth muscle contraction(血管平滑肌收缩通路)、Insulin signaling pathway(胰岛素信号通路)和Steroid hormonebiosynthesis(类固醇激素的生物合成)。
所述“失和”生物分子网络为钙离子稳态相关“失和”网络(Net.Calcium)、血管新生相关“失和”网络(Net.Angiogenesis)、血管平滑肌收缩相关“失和”网络(Net.VSMContraction)、能量代谢相关“失和”网络(Net.EnergyMetabolism)或类固醇激素合成相关“失和”网络(Net.Steroid)。
与现有技术相比,本发明具有以下有益效果:
(1)本发明方法结合中药中入血成分与靶蛋白的对接得分值和成分在中药中的相对含量,有效辨识中药化学组成与机体生物分子网络之间的相关性,定量描述中药对其相关“失和”生物分子网络的调控力度。
(2)本发明以一系列中药(如类方)为载体,在系统剖析各中药化学组成及生物学效应的基础上,通过比较各中药之间的内在异同,探讨诠释中药的配伍科学内涵,丰富和拓展了中药配伍规律的研究思路。
附图说明
图1为本发明中药对复杂性疾病所致的“失和”生物分子网络调控能力的评价方法的流程图。
图2为四逆汤类方给药对大鼠心肌梗死范围的影响;
(A)TTC染色法评价心肌梗死面积的统计结果。其中model、SIN和RSIN,n=6;FSIN、TSIN和ZSIN,n=7;**代表与模型对照组相比p<0.01;
(B)心肌梗死大鼠的代表性心脏切片图,正常心肌呈红色,梗死心肌呈白色;其中Sham为假手术组大鼠的心脏切片图。
图3为本发明方法获得的四逆汤类方相关的急性心肌梗塞所致的“失和”生物分子网络;
其中,(1)能量代谢相关“失和”网络(Net.EnergyMetabolism);(2)血管新生相关“失和”网络(Net.Angiogenesis);(3)血管平滑肌收缩相关“失和”网络(Net.VSMContraction);(4)钙离子稳态相关“失和”网络(Net.Calcium);(5)类固醇激素合成相关“失和”网络(Net.Steroid)。
图4为四逆汤类方对急性心肌梗塞所致的相关“失和”生物分子网络的调控差异;
其中,SIN为四逆汤组;RSIN为四逆加人参汤组;FSIN为茯苓四逆汤组;TSIN为通脉四逆汤组;ZSIN通脉四逆加猪胆汁汤组;
Net.EnergyMetabolism为能量代谢相关“失和”网络;Net.Angiogenesis为血管新生相关“失和”网络;Net.VSMContraction为血管平滑肌收缩相关“失和”网络;Net.Calcium为钙离子稳态相关“失和”网络;Net.Steroid为类固醇激素合成相关“失和”网络。
图5为四逆汤类方对血管缺失模型斑马鱼的影响;
其中,正常对照组:Control;模型对照组:Model;阳性药对照组(30μM):Ast;(A)四逆汤类方对血管缺失斑马鱼肠下血管面积影响的表型图;虚线区域为肠下血管面积;(B)四逆汤类方对血管缺失斑马鱼肠下血管分支数影响的表型图;白色箭头所指为肠下血管分支数;(C)四逆汤类方对血管缺失模型斑马鱼肠下血管面积的保护率。
图6为四逆汤类方对心衰斑马鱼的影响;
其中,正常对照组:Control;模型对照组:Model;阳性药对照组(0.1μg/mL):Digoxin;(A)四逆汤类方对心衰斑马鱼心脏静脉淤血改善表型图;白色箭头所指为静脉淤血处;(B)四逆汤类方对心衰斑马鱼心脏静脉淤血的保护率;(C)四逆汤类方对心衰斑马鱼血流速度的保护率。
具体实施方式
一种中药对复杂性疾病所致“失和”网络整合调节作用的定量评价方法,具体步骤如下:
(1)获取中药对复杂性疾病的潜在作用靶点,利用所述潜在作用靶点构建所述“失和”生物分子网络;
A)获得中药作用于复杂性疾病的芯片数据。
以大鼠为研究对象,构建大鼠的急性心肌梗塞模型,对患有急性心肌梗塞的大鼠实施灌胃给药(即中药)进行疾病的治疗,取大鼠的心肌组织用于基因芯片检测,获得芯片数据。
①实验组别及剂量设置
根据《伤寒论》中人的日给药量换算得到动物实验中大鼠的日给药量,四逆汤类方粉末,包括四逆汤、四逆加人参汤、茯苓四逆汤、通脉四逆汤和通脉四逆加猪胆汁汤均自制。深黄色粉末,常温干燥保存,临用前方剂称重,加适量Milli-Q水溶解、混匀并定容至方剂的给药体积。
设置实验的组别,分为假手术组Sham(构建急性心肌梗塞模型时,大鼠只穿线不结扎;不进行给药,给予相应量的Milli-Q水)、模型对照组Model(构建急性心肌梗塞模型时,大鼠结扎;不进行给药,给予相应量的Milli-Q水)、中药给药组SNs(构建急性心肌梗塞模型时,大鼠结扎;进行给药)。
②大鼠急性心肌梗塞模型的制备
SD大鼠,禁食过夜,自由饮水。实验当天腹腔注射360mg/kg水合氯醛麻醉,每只肌注0.1mL盐酸利多卡因抗心律失常。开胸前连接ALC-V8呼吸机,开胸后迅速于大鼠左侧第三和第四肋间小心掏出心脏,并用6/0线结扎冠状动脉左前降支。结扎后小心地将心脏放回胸腔,排尽胸腔内空气,用4/0线缝合肌肉组织和皮肤。每只大鼠肌肉注射青霉素16万单位。手术后第二天,通过眼眶取血方法对手术后存活的大鼠进行采血,血液室温静置30min后,4℃下4000rpm/min离心15min,取上清测定大鼠血清中心肌肌钙蛋白T(cTnT)的浓度。根据cTnT浓度判断大鼠急性心肌梗塞模型的成功与否。
TTC用心肌组织取材:SD大鼠称重,腹腔注射360mg/kg水合氯醛麻醉后,腹主动脉取血后,用20mL预冷的生理盐水对心脏进行灌流后,剪下心脏于预冷的生理盐水中清洗。剔除结缔组织和右心室等后取结扎线以下的部位置于预冷的2mL冻存管中,再-80℃贮存备用。
基因芯片检测用心肌组织取材:所用手术器械使用之前首先用含有0.1%DEPC水处理,再进行高压灭菌并烘干。用20mL预冷的含0.1%DEPC生理盐水对心肌进行灌流后,剪下心脏置于预冷的含0.1%DEPC生理盐水中清洗,并剔除结缔组织和右心室等。取结扎线以下梗死和正常心肌交界区的心肌组织两份,每份约150mg,用铝箔纸包装后分别置于预冷的2mL冻存管中,迅速于液氮中冻存,之后再-80℃贮存备用。
③基因芯片实验及数据预处理
采用的芯片为Affymetrix Rat Genome 230 2.0 Array,该芯片平台为大鼠全基因组芯片。应用mirVanaTM RNA Isolation Kit(Applied Biosystem p/n AM1556)试剂盒提取心肌样品的总RNA。利用NanoDrop ND-2100对样品总RNA进行定量并采用AgilentBioanalyzer 2100检测RNA完整性。
质检合格后,采用Affymetrix表达谱芯片配套试剂盒,参照标准操作流程进行样本标记、芯片的杂交以及洗脱等。经Affymetrix Scanner 3000扫描得到原始图像并利用AGCC软件(Affymetrix GeneChip Command Console version4.0)提取得到原始数据(.CEL格式)。原始数据导入ArrayTrack软件中,应用MAS5算法进行背景校正,再对所有数据进行中位数1000的标准化处理,方便后续芯片之间的比较研究。
④基因芯片实验结果验证
对同一批RNA样本,采用实时定量PCR(Real-time Quantitative PCR)实验验证芯片数据结果的可靠性。Txn2、Pkp2、Cxcl12、Cacnb2和GAPDH的引物由上海祥音生物科技有限公司合成,其中GAPDH为内参基因。
实时定量PCR验证实验过程分为3个步骤。
1)RNA逆转录合成cDNA。采用AMBION动物试剂盒抽提大鼠心脏组织中RNA。样品中RNA浓度采用NanoDrop 2000 Spectrophotometer进行测定,再按照逆转录试剂盒的标准操作流程逆转录合成cDNA。
2)样品实时定量PCR测定。引物经适量DEPC水溶解后得到100μM的引物储备液,再根据QuantiFast SYBR Green PCR Kit试剂盒的标准操作流程进行样品的实时定量PCR检测。
3)应用2-ΔΔCt方法计算基因转录水平表达的改变强度(Fold Change,FC),最后比较实时定量PCR和基因芯片的结果。
B)采用cMAP法得到中药作用复杂性疾病的潜在相关靶蛋白。
将上述获取的芯片数据代入cMAP(Connectivity Map)数据库中进行分析,获得潜在的相关靶蛋白。
(a)计算探针的FC(SNs/Model)值,当探针的FC(SNs/Model)绝对值大于1.2,且p<0.05时,定义为显著差异表达的探针;
(b)根据探针改变的幅度大小,挑选中药芯片数据中上调或下调幅度最大的显著差异表达的探针与cMAP数据库中的芯片数据进行相似度匹配,找到与中药相关性最强的化合物;
(c)根据STITCH4.0数据库中化合物-蛋白质相互作用的数据信息,得到中药潜在作用蛋白靶点。
C)采用EoR法得到中药作用复杂性疾病的潜在相关基因。
根据芯片数据计算节点的回调状态;
式(6)中,RL′表示回调状态;LogE为各实验组的基因表达水平取log2对数后的均值;Edrug为中药给药组的芯片数据;Esham为假手术组的芯片数据;Emodel为模型对照组的芯片数据。
再计算节点的回复调控效率(Efficiency of Recovery regulation,EoR)
EoR=100%-|100%-RL′| (7)
式(7)中,EoR表示为回复调控效率(Efficiency of Recovery regulation),是用于表征节点(芯片数据)回调效率的指标;RL′表示回调状态。
根据计算,EoR大于0的节点对复杂性疾病所致“失和”生物分子网络具有复衡贡献,而EoR小于等于0的节点对其贡献值为0。故,定义EoR值大于0的节点为获得的潜在相关基因。
D)整合相关靶蛋白和相关基因,得到中药作用复杂性疾病的靶蛋白列表,再通过靶蛋白列表构建中药相关的复杂性疾病所致的“失和”生物分子网络。
(a)对潜在相关靶基因和潜在相关靶蛋白之间的对应关系进行合并,去重;
(b)应用DAVID Bioinformatics Resources 6.7软件对靶蛋白列表进行通路富集分析,找到所述复杂性疾病相关的生物学信号通路,定义p<0.05为显著富集通路;
(c)以KEGG数据库中所述显著富集通路为基础,构建所述“失和”生物分子网络。
(2)计算获得中药对复杂性疾病所致“失和”生物分子网络的调控力度值。
a)利用中药中入血成分与“失和”生物分子网络中靶蛋白的对接得分值,计算“失和”生物分子网络中各连接的权重;
式(1)中,Scomp表示中药中入血成分与网络中蛋白进行分子对接的得分值;Sreference表示网络中蛋白质复合结构的原始配体与蛋白进行分子对接的最低得分值;EW表示网络中的任何一个节点到下游直接节点的连接的权重;
b)将权重代入式(2)和式(3)中,计算获得网络效能和网络流;
式(2)中,dwij为网络中节点i与节点j之间的最短距离,即节点i到节点j之间所经过的边的权重之和最小;N为网络中所有节点的集合;NE为网络效能,即网络中任何一个节点i到下游节点j的最短路径的倒数和;
式(3)中,dwim为网络中节点i与节点m之间的最短距离,即节点i到节点m之间所经过的边的权重之和最小;N为网络中所有节点的集合;NF为网络流,即网络中任何一个节点i到网络中叶节点m的最短路径的倒数和;LeafNode表示叶节点,是指网络中出度为0的节点;
c)将网络效能和网络流代入式(4)中,计算获得中药中入血成分对所述“失和”生物分子网络的调控能力;
式(4)中,当所述“失和”生物分子网络中所有连接的权重均为初始值1时,计算得到NEmax和NFmax;NEcomp表示中药中单个化合物对所述“失和”生物分子网络的网络效能,NFcomp表示中药中单个化合物对所述“失和”生物分子网络的网络流;RE表示中药中入血成分对所述“失和”生物分子网络的调控能力;
d)将“失和”生物分子网络的调控能力代入式(5)中,计算得到中药对复杂性疾病所致“失和”生物分子网络的调控力度值;
式(5)中,k表示中药中的入血成分,G表示中药中所有入血成分的集合,Ak表示成分k在中药所有入血成分中所占的峰面积百分比,FRE表示中药对复杂性疾病所致“失和”生物分子网络的调控力度值。
其中,入血成分在中药所有入血成分中所占的峰面积百分比根据实验数据获得,具体的实验方法如下:
动物实验:禁食过夜的SD大鼠用于动物实验。灌胃前,收集空白血样约1mL,然后随机分成5组,根据组别分别灌胃给药四逆汤类方,给药量3.2mL/100g体重(具体见表1),给药60min后,乙醚麻醉,通过眼眶静脉收集给药后的血样,4000rpm离心15min,取上清血浆,置-80℃冰箱保存。
液质联用分析用血浆样品制备:取血浆样品500μL,加入乙腈1000μL,涡旋30s,10000rpm离心10min。35℃下离心浓缩,残渣中加入100μL甲醇复溶,10000rpm离心10min,取上清液备用。
液相色谱分离:应用配备二元泵的Agilent 1100型液相色谱仪(Agilent公司,德国)和Zorbax Eclipse XDB-C18色谱柱(5μm,4.6mm×250mm,Agilent公司,美国)开展色谱分离实验。流动相分别为0.05%甲酸-水(A)和0.05%甲酸-乙腈(B)。采用梯度洗脱,梯度设置如下:0-50min,5-30%B;50-75min,30-50%B;75-80min,50-95%B;80-90min,95%B。流速为0.6mL/min,柱温30℃,类方提取物样品进样体积为20μL,血浆样品进样体积为30μL。光电二极管阵列检测器(Photodiode array,PDA)波长设置为190nm到400nm。
质谱分析:应用配备电喷雾离子化接口(ESI)的Finnigan LCQ Deca XPplus离子阱质谱仪(Thermo公司,美国)采集多级质谱数据。仪器实验参数如下:全扫描范围,质荷比(m/z)100-1500;正离子模式电离电压+4.00kV,负离子模式为–3.00kV;毛细管温度350℃;鞘气(N2)60arb;辅助气(N2)20arb。毛细管电压19V;透镜管电压25V。多级质谱数据由Xcalibur 2.1,Qual Browser数据工作站处理分析(Thermo公司,美国)。
表1四逆汤类方动物实验口服单次给药量
a日给药量指每千克大鼠每天服用的剂量,大鼠日给药量是根据体重和体表面积按人体等效剂量进行换算得来。
b口服单次给药量是本次实验中的确切给药量,综合大鼠单次最大口服体积和四逆汤类方的溶解性,口服单次给药量为日给药量的8倍(3.2mL/100g体重)。
(3)比较不同中药之间的调控力度值,评价各中药对所述“失和”生物分子网络的调控能力。
对每个网络所有中药对应的FRE值进行归一化处理,并将归一化后的结果分为A-F六个等级(A:0.8-1;B:0.7-0.8;C:0.6-0.7;D:0.5-0.6;E:0.5-0.4;F:0-0.4)。
采用随机化检验(randomization test)方法检验基于成分与蛋白对接模拟得分值设计的网络评价结果是否可靠。在随机检验过程中,将docking score值进行随机排序,按照计算流程重新计算中药对网络的调控能力,该过程重复10万次。
最终,将计算得到的FRE值与随机获得的10万个FRE值进行单样本t检验。
结果显示p值趋近于0,表明基于成分与蛋白对接模拟得分值设计的网络评价结果是可靠的(不是随机结果)。
应用例1
针对四逆汤类方相关的急性心肌梗塞所致的“失和”生物分子网络,本应用例采用本发明方法来评价五首四逆汤类方对其“失和”生物分子网络的调控力度,并通过斑马鱼实验来验证本发明方法。
一、评价五首四逆汤类方对急性心肌梗塞所致的“失和”生物分子网络的调控力度
(1)获取四逆汤类方对急性心肌梗塞的潜在作用靶点,利用所述潜在作用靶点构建“失和”生物分子网络
A)获取四逆汤类方作用于急性心肌梗塞的芯片数据;
以大鼠为研究对象,构建其急性心肌梗塞模型,对患有急性心肌梗塞的大鼠实施灌胃给药进行疾病的治疗,取大鼠的心肌组织用于基因芯片检测,获得芯片数据。
①实验组别及剂量的设置
根据《伤寒论》中人的日给药量换算得到动物实验中大鼠的日给药量,给药采用的中药共有五种,分别为:四逆汤SIN、四逆加人参汤RSIN、茯苓四逆汤FSIN、通脉四逆汤TSIN和通脉四逆加猪胆汁汤ZSIN,药材组成详见表2;组别设置和给药剂量见表3。
表2四逆汤类方药材组成
临用前每首类方称重,加适量Milli-Q水溶解、混匀并定容至每首类方的给药体积为1mL/100g。
表3分组情况及药物剂量、浓度和给药体积
②大鼠急性心肌梗塞模型的制备
SD大鼠,禁食过夜,自由饮水。实验当天腹腔注射360mg/kg水合氯醛麻醉,每只肌注0.1mL盐酸利多卡因抗心律失常。开胸前连接ALC-V8呼吸机,开胸后迅速于大鼠左侧第三和第四肋间小心掏出心脏,并用6/0线结扎冠状动脉左前降支,其中假手术组大鼠只穿线不结扎。结扎后小心将心脏放回胸腔,排尽胸腔内空气,用4/0线缝合肌肉组织和皮肤。每只大鼠肌肉注射青霉素16万单位。手术后第二天,通过眼眶取血方法对手术后存活的大鼠进行采血,血液室温静置30min后4℃下4000rpm/min离心15min,取上清测定大鼠血清中心肌肌钙蛋白T(cTnT)的浓度。根据cTnT浓度判断大鼠急性心肌梗塞模型的成功与否。
③取材
TTC用心肌组织取材:SD大鼠称重,腹腔注射360mg/kg水合氯醛麻醉后,腹主动脉取血后,用20mL预冷的生理盐水对心脏进行灌流后,剪下心脏于预冷的生理盐水中清洗。剔除结缔组织和右心室等后取结扎线以下的部位置于预冷的2mL冻存管中,再-80℃贮存备用。
基因芯片检测用心肌组织取材:所用手术器械使用之前首先用含有0.1%DEPC水处理,再进行高压灭菌并烘干。用20mL预冷的含0.1%DEPC生理盐水对心肌进行灌流后,剪下心脏置于预冷的含0.1%DEPC生理盐水中清洗,并剔除结缔组织和右心室等。取结扎线以下梗死和正常心肌交界区的心肌组织两份,每份约150mg,用铝箔纸包装后分别置于预冷的2mL冻存管中,迅速于液氮中冻存,之后再-80℃贮存备用。
④心肌梗死范围测定
待心肌组织冻僵硬后,沿冠状面均匀切成4-5片心肌组织并置于2%TTC溶液中,37℃水浴避光加热5min左右,在孵育过程中用镊子小心翻动心肌切片,使得心肌切片的两面都能充分与2%TTC反应。最终心肌非梗死区呈现红色,梗死区则是白色,应用美国国立卫生研究院(NIH)发布的公共图像统计软件ImageJ 1.48对心肌梗死范围(Infarct Size)进行统计。
⑤四逆汤类方对急性心肌梗塞大鼠心肌梗死范围的影响
急性心肌梗塞模型对照组和给药组大鼠给药结束后,取心脏以2%TTC染色法评价大鼠的心肌梗死范围。急性心肌梗塞模型对照组大鼠平均梗死范围为38.7±8.8%。四逆汤、茯苓四逆汤、通脉四逆汤和通脉四逆加猪胆汁汤灌胃给药7d后,能显著减少大鼠冠状动脉左前降支结扎引起的心肌梗死范围,但各中药影响程度各不同,具体参见图2A。四逆加人参汤灌胃给药7d后能减轻大鼠心肌梗死的程度,但未达到统计学显著性差异(RSIN vsModel,p=0.066)。图2B展示了TTC染色标记的大鼠心肌梗死范围。
⑥基因芯片实验及数据预处理
本发明中所采用的芯片为Affymetrix Rat Genome 230 2.0 Array,该芯片平台为大鼠全基因组芯片,共包括31099个探针。共检测35个样本,分为假手术组(Sham)和模型对照组(Model)、四逆汤组(SIN)、四逆加人参汤组(RSIN)、茯苓四逆汤组(FSIN)、通脉四逆汤组(TSIN)和通脉四逆加猪胆汁汤组(ZSIN),每组样本数均为5。
应用mirVanaTM RNA Isolation Kit(Applied Biosystem p/n AM1556)试剂盒提取心肌样品的总RNA。利用NanoDrop ND-2100对样品总RNA进行定量并采用AgilentBioanalyzer 2100检测RNA完整性。
质检合格后,采用Affymetrix表达谱芯片配套试剂盒,参照标准操作流程进行样本标记、芯片的杂交以及洗脱等。经Affymetrix Scanner 3000扫描得到原始图像并利用AGCC软件(Affymetrix GeneChip Command Console version4.0)提取得到原始数据(.CEL格式)。原始数据导入ArrayTrack软件中,应用MAS5算法进行背景校正,再对所有数据进行中位数1000的标准化处理,方便后续芯片之间的比较研究。
⑦基因芯片实验结果验证
对同一批RNA样本,采用实时定量PCR(Real-time Quantitative PCR)实验验证芯片数据结果的可靠性。Txn2、Pkp2、Cxcl12、Cacnb2和GAPDH的引物由上海祥音生物科技有限公司合成,其中GAPDH为内参基因。
实时定量PCR验证实验过程分为3个步骤。
1)RNA逆转录合成cDNA。采用AMBION动物试剂盒抽提大鼠心脏组织中RNA。样品中RNA浓度采用NanoDrop 2000 Spectrophotometer进行测定,再按照逆转录试剂盒的标准操作流程逆转录合成cDNA。
2)样品实时定量PCR测定。引物经适量DEPC水溶解后得到100μM的引物储备液,再根据QuantiFast SYBR Green PCR Kit试剂盒的标准操作流程进行样品的实时定量PCR检测。
3)应用2-ΔΔCt方法计算基因转录水平表达的改变强度(Fold Change),最后比较实时定量PCR和基因芯片的结果。
B)采用cMAP法得到四逆汤类方作用急性心肌梗塞的潜在相关靶蛋白
将上述获取的芯片数据代入cMAP(Connectivity Map)数据库中进行分析,获得潜在相关靶蛋白。
(a)将中药给药组分别与模型对照组相比较,当探针的FC(SNs/Model)绝对值大于1.2,且p<0.05时,定义为显著差异表达的探针;
(b)根据探针改变的幅度大小,挑选每首四逆汤类方芯片数据中上或下调幅度最大的各150个显著差异表达的探针与cMAP数据库中的芯片数据进行相似度匹配,找到与四逆汤类方相关性最强的化合物;
(c)根据化合物-蛋白质相互作用数据(STITCH4.0数据库),得到五首四逆汤类方潜在作用蛋白靶点。
C)采用EoR方法得到四逆汤类方作用急性心肌梗塞的潜在相关基因,具体为:
根据芯片数据计算节点的回调状态;
式(6)中,RL′表示回调状态;LogE为各实验组的基因表达水平取log2对数后的均值;Edrug为中药给药组的芯片数据;Esham为假手术组的芯片数据;Emodel为模型对照组的芯片数据。
再计算节点的回复调控效率(Efficiency of Recovery regulation,EoR)
EoR=100%-|100%-RL′| (7)
式(7)中,EoR表示为回复调控效率(Efficiency of Recovery regulation),是用于表征节点(芯片数据)回调效率的指标;RL′表示回调状态。
根据计算,EoR大于0的节点对复杂性疾病所致“失和”生物分子网络具有复衡贡献,而EoR小于等于0的节点对其贡献值为0。故,定义EoR值大于0的节点为获得的潜在相关基因。
D)整合相关靶蛋白和相关基因,得到四逆汤类方作用急性心肌梗塞的靶蛋白列表,再通过靶蛋白列表构建四逆汤类方相关的急性心肌梗塞所致的“失和”生物分子网络。
a)对潜在相关靶基因和潜在相关靶蛋白之间的对应关系进行合并,去重;
b)在通晓冠心病主要病理环节的基础上,对四逆汤类方治疗急性心肌梗塞相关的潜在作用靶点进行整合分析(见表4)。
表4冠心病发生发展相关的病理过程
c)应用DAVID Bioinformatics Resources 6.7软件中Functional AnnotationTool工具下的Pathways模块对靶蛋白列表进行通路富集分析,研究其涉及的生物学信号通路。其中Pathways模块下的分类选用KEGG数据库,显著通路定义为p<0.05。
d)综合考虑冠心病相关的各病理环节和关键生物分子显著富集通路的p值,选以KEGG通路Calcium signaling pathway(钙离子信号通路)、ErbB signaling pathway(ErbB信号通路)、Vascular smooth muscle contraction(血管平滑肌收缩通路)、Insulinsignaling pathway(胰岛素信号通路)和Steroid hormone biosynthesis(类固醇激素的生物合成)为基础,构建四逆汤类方相关的急性心肌梗塞所致的“失和”生物分子网络(见图3)
(2)计算获得四逆汤类方对复杂性疾病所致“失和”生物分子网络的调控力度值。
a)利用四逆汤类方中入血成分与“失和”生物分子网络中靶蛋白的对接得分值,计算“失和”生物分子网络中各连接的权重;
在PDB蛋白数据库中搜索相关“失和”生物分子网络中四逆汤类方作用的靶蛋白的三维结构,应用MOE软件完成蛋白质和化合物的分子对接工作,得到四逆汤类方中入血成分和生物分子网络中靶蛋白的对接得分值;
设定网络中所有连接的初始权重为1,最大权重为200。每一个网络中蛋白质复合结构中的原始配体与该蛋白进行分子对接模拟,选取最低得分值作为该网络的Sreference;
式(1)中,Scomp表示中药中入血成分与网络中蛋白进行分子对接的得分值;Sreference表示网络中蛋白质复合结构的原始配体与蛋白进行分子对接的最低得分值;EW表示网络中的任何一个节点i到下游直接节点j的连接的权重;
b)将权重代入式(2)和式(3)中,计算获得网络效能和网络流;
式(2)中,dwij为网络中节点i与节点j之间的最短距离,即节点i到节点j之间所经过的边的权重之和最小;N为网络中所有节点的集合;NE为网络效能,即网络中任何一个节点i到下游节点j的最短路径的倒数和;
式(3)中,dwim为网络中节点i与节点m之间的最短距离,即节点i到节点m之间所经过的边的权重之和最小;N为网络中所有节点的集合;NF为网络流,即网络中任何一个节点i到网络中叶节点m的最短路径的倒数和;LeafNode表示叶节点,是指网络中出度为0的节点;
c)将网络效能和网络流代入式(4)中,计算获得四逆汤类方中入血成分对所述“失和”生物分子网络的调控能力;
式(4)中,当所述“失和”生物分子网络中所有连接的权重均为初始值1时,计算得到NEmax和NFmax;NEcomp表示中药中单个化合物对所述“失和”生物分子网络的网络效能,NFcomp表示中药中单个化合物对所述“失和”生物分子网络的网络流;RE表示中药中入血成分对所述“失和”生物分子网络的调控能力;
d)将“失和”生物分子网络的调控能力代入式(5)中,计算得到四逆汤类方对急性心肌梗塞所致“失和”生物分子网络的调控力度值;
式(5)中,k表示中药中的入血成分,G表示四逆汤类方中所有入血成分的集合,Ak表示成分k在四逆汤类方所有入血成分中所占的峰面积百分比,FRE表示四逆汤类方对急性心肌梗塞所致“失和”生物分子网络的调控力度值。
(3)比较不同中药之间的调控力度值,评价各中药对所述“失和”生物分子网络的调控能力。
对每个网络所有四逆汤类方对应的FRE值进行归一化处理,并将归一化后的结果分为A-F六个等级(A:0.8-1;B:0.7-0.8;C:0.6-0.7;D:0.5-0.6;E:0.5-0.4;F:0-0.4)。
采用随机化检验(randomization test)方法检验基于成分与蛋白对接模拟得分值设计的网络评价结果是否可靠。在随机检验过程中,将docking score值进行随机排序,按照计算流程重新计算中药对网络的调控能力,该过程重复10万次。
最终,将计算得到的FRE值与随机获得的10万个FRE值进行单样本t检验。
结果显示p值趋近于0,表明基于成分与蛋白对接模拟得分值设计的网络评价结果是可靠的(不是随机结果)。
实验结果如下:
随机化检验结果显示各个中药对不同网络计算得到的FRE值与相应随机获得的10万个FRE值进行单样本t检验时,p值趋近于0,表明基于成分与蛋白对接模拟得分值设计的网络评价结果是可靠的。
FRE值归一化后的四逆汤类方各个中药结果如图4所示,该图清楚地展示了各个中药对其相关“失和”生物分子网络调控强弱的侧重。五首四逆汤类方对钙离子稳态相关“失和”网络(Net.Calcium)调控力度一致,均能很好的调控缺氧缺血导致的细胞内钙离子紊乱。而四逆加人参汤和茯苓四逆汤对促血管新生相关的血管新生相关“失和”网络(Net.Angiogenesis)具有更强的调控力度。通脉四逆加猪胆汁汤对能量代谢相关“失和”网络(Net.EnergyMetabolism)、血管平滑肌收缩相关“失和”网络(Net.VSMContraction)和类固醇激素合成相关“失和”网络(Net.Steroid)具有更强的调控作用。而通脉四逆汤相比四逆汤而言,对血管新生相关“失和”网络、能量代谢相关“失和”网络、血管平滑肌收缩相关“失和”网络和类固醇激素合成相关“失和”网络的调控力度均有所增加。
如图4所示,四逆加人参汤和茯苓四逆汤对促血管新生相关的血管新生相关“失和”网络(Net.Angiogenesis)具有比其他方剂更强的调控力度。早在上个世纪末,有学者通过基因敲除手段发现ErbB受体在心血管系统有表达,而且在心脏的发育过程中起作用。一方面,HER2被Heregulin激活后,可通过促进内皮细胞生长、增殖、迁移,促进内皮细胞募集平滑肌细胞等作用直接调控血管生成。另一方面,HER2激活后还可通过VEGF依赖性的方式间接调控血管新生。进一步,实验研究表明抗HER2单克隆抗体可通过下调HER2的表达,弱化HER2的二聚化等机制发挥抗血管新生的作用。
如图4所示,通脉四逆加猪胆汁汤对能量代谢相关“失和”网络(Net.EnergyMetabolism)、血管平滑肌收缩相关“失和”网络(Net.VSMContraction)和类固醇激素合成相关“失和”网络(Net.Steroid)具有比其他方剂更强的调控作用。研究表明胰岛素抵抗与心力衰竭之间存在着密切的联系,胰岛素信号调控在缺血性心脏保护以及延缓心衰的发生发展中具有重要作用。心肌缺血引发胰岛素抵抗后,不仅上调了脂肪酸的代谢,同时也抑制了葡萄糖的摄取和代谢,导致心肌耗氧量增加,心肌能力代谢紊乱,从而加剧冠心病的发展,诱发冠心病心力衰竭。此外,当心肌缺血时,活性氧的增加、胞浆内Ca2+浓度的超载和能量代谢紊乱等会引发血管内皮功能失调,NO释放量减少,从而抑制鸟苷酸环化酶的活化,减少平滑肌细胞内环鸟苷酸(cGMP)的产生,导致平滑肌舒张功能障碍,收缩异常,血压升高甚至出现痉挛。而正常血液循环的维持有赖于正常的心肌舒缩功能、正常的血管舒缩功能和足够的循环血量。
针对上述结果,选择以斑马鱼血管缺失模型来验证RSIN和FSIN是否比其他类方具有更强调控血管新生相关“失和”网络的作用。选以斑马鱼心衰模型来验证ZSIN是否比其他类方具有更强调控能量代谢相关“失和”网络、血管平滑肌收缩相关“失和”网络和类固醇激素合成相关“失和”网络的作用。具体内容如下:
1、实验材料
辛伐他汀:日本TCI公司,批号:404-520-2;黄芪甲苷:阿拉丁试剂(上海)有限公司,批号:G1325036;盐酸维拉帕米:上海晶纯实业有限公司,批号:E1216038;地高辛:上海晶纯实业有限公司,批号:D102298。
四逆汤类方粉末:四逆汤、四逆加人参汤、茯苓四逆汤、通脉四逆汤和通脉四逆加猪胆汁汤均自制。深黄色粉末,常温干燥保存,临用前用超纯水配制成工作液,现配现用。
20mM盐酸维拉帕米母液:临用前称取一定量的盐酸维拉帕米粉末,加超纯水配置成20mM母液;0.1mg/mL地高辛母液:临用前称取一定量的地高辛粉末,加DMSO配置成0.1mg/mL母液;70mM辛伐他汀母液:临用前称取一定量的辛伐他汀粉末,加DMSO配置成70mM工作液;30mM黄芪甲苷母液:称取一定量的黄芪甲苷粉末,加DMSO配置成30mM母液,-20℃保存。
解剖显微镜:SZX7,Olympus;心跳血流分析系统:ZebraLab 3.3(PB2084C),Viewpoint life science;精密电子天平:CP214,Ohaus;电动聚焦连续变倍荧光显微镜:AZ100,Nikon。
斑马鱼的品系与来源:实验动物使用许可证号为:SYXK(浙)2012-0171。饲养管理符合国际AAALAC认证的要求。1)野生型AB系斑马鱼,以自然成对交配繁殖方式进行。年龄为受精后2天或3天(2dpf/3dpf,用于四逆汤类方心衰药效评价实验)。2)转基因Fli-1品系斑马鱼,以自然成对交配繁殖方式进行。年龄为受精后1天(1dpf,用于四逆汤类方促血管形成药效评价实验)。
2、四逆汤类方对斑马鱼促血管再生的药效学检验
随机选取转基因Fli-1系1dpf斑马鱼于六孔板中,每孔30尾,用70nM辛伐他汀诱导Fli-1系斑马鱼建立微血管缺失模型,3h后移除诱导剂,分别水溶给予不同浓度的类方药液,四逆汤类方给药量根据其最大无可见毒害剂量作用水平(NOAEL)和药效预实验结果设定。阳性对照组同样条件下给予30μM黄芪甲苷(Ast),每孔3mL药液,同时设置正常对照组和模型对照组。给药24h后移除药物,再继续在斑马鱼生化培养箱中培养24h。培养结束后,每组随机取10尾斑马鱼在显微镜下观察肠下血管面积和分支数,拍照并保存图片。
3、四逆汤类方对斑马鱼心衰模型的药效学检验
随机选取野生型AB系斑马鱼于六孔板中,每孔30尾,200μM维拉帕米诱发AB系斑马鱼心衰模型,分别水溶给予不同浓度的四逆汤类方药物,类方给药量根据其NOAEL水平和药效预实验结果设定。阳性对照组给予0.1μg/mL地高辛,每孔3mL药液,给药4h,同时设置正常对照组和模型对照组。给药结束后,每组随机选取10尾斑马鱼拍照并采集数据,分析统计斑马鱼心脏静脉淤血面积;随机选取10尾并置于心跳血流分析系统下录制斑马鱼血流视频,分析斑马鱼的血流速度。
验证结果如下:
从图5的统计结果来看,RSIN和FSIN对血管缺失斑马鱼肠下血管面积的保护率高达90%以上,TSIN和ZSIN的保护率为60%左右,而SIN的保护率接近50%。因此,该实验证明五首四逆汤类方都具有促血管生成的作用,但RSIN和FSIN的效果最佳,这可能与RSIN和FSIN对血管新生相关“失和”网络的调控力度更大有关。
图6(B)的统计结果显示,ZSIN对心衰斑马鱼心脏静脉淤血的改善效果最佳。图6(C)结果显示,ZSIN和TSIN对心衰斑马鱼血流速度的恢复效果明显优于其他三首方剂。因此,改结果该实验初步证明了ZSIN对治疗心力衰竭具有良好的效果,而且其作用明显优于SIN、RSIN和FSIN,这可能与ZSIN对能量代谢相关“失和”网络、血管平滑肌收缩相关“失和”网络和类固醇激素合成相关“失和”网络同时具有强调控能力有关。
最后,还需要特别注意的是,以上所举例子仅是本发明的具体实施例子。显然,本发明不仅仅限于以上实施例子,还可以有许多变通的情况。本领域的技术人员从本发明公开的内容直接推导出或联想到的所有变通情况,均认为是本发明的保护范围。
Claims (9)
1.一种中药对复杂性疾病所致“失和”网络整合调节作用的定量评价方法,其特征在于,包括以下步骤:
(1)获取中药对复杂性疾病的潜在作用靶点,利用所述潜在作用靶点构建所述“失和”生物分子网络;
(2)计算获得中药对所述“失和”生物分子网络的调控力度值;
a)利用所述中药中入血成分与所述“失和”生物分子网络中靶蛋白的对接得分值,计算所述“失和”生物分子网络中各连接的权重;
式(1)中,Scomp表示中药中入血成分与网络中蛋白进行分子对接的得分值;Sreference表示网络中蛋白质复合结构的原始配体与蛋白进行分子对接的最低得分值;EW表示网络中的任何一个节点到下游直接节点的连接的权重;
b)将所述权重代入式(2)和式(3)中,计算获得网络效能和网络流;
式(2)中,dwij为网络中节点i与节点j之间的最短距离,即节点i到节点j之间所经过的边的权重之和最小;N为网络中所有节点的集合,NE为网络效能,即网络中任何一个节点i到下游节点j的最短路径的倒数和;
式(3)中,dwim为网络中节点i与节点m之间的最短距离,即节点i到节点m之间所经过的边的权重之和最小;N为网络中所有节点的集合,NF为网络流,即网络中任何一个节点i到网络中叶节点m的最短路径的倒数和;LeafNode表示叶节点,是指网络中出度为0的节点;
c)将所述网络效能和网络流代入式(4)中,计算获得中药中入血成分对所述“失和”生物分子网络的调控能力;
式(4)中,当所述“失和”生物分子网络中所有连接的权重均为初始值1时,计算得到NEmax和NFmax;NEcomp表示中药中单个入血成分对所述“失和”生物分子网络的网络效能,NFcomp表示中药中单个入血成分对所述“失和”生物分子网络的网络流;RE表示中药中入血成分对所述“失和”生物分子网络的调控能力;
d)将所述“失和”生物分子网络的调控能力代入式(5)中,计算得到中药对复杂性疾病所致“失和”生物分子网络的调控力度值;
式(5)中,k表示中药中的入血成分,G表示中药中所有入血成分的集合,Ak表示成分k在中药所有入血成分中所占的峰面积百分比,FRE表示中药对复杂性疾病所致“失和”生物分子网络的调控力度值;
(3)比较不同中药之间的调控力度值,评价各中药对所述“失和”生物分子网络的调控能力。
2.如权利要求1所述的评价方法,其特征在于,所述中药为四逆汤类方。
3.如权利要求1所述的评价方法,其特征在于,所述中药为四逆汤、四逆加人参汤、茯苓四逆汤、通脉四逆汤或通脉四逆加猪胆汁汤。
4.如权利要求1所述的评价方法,其特征在于,所述复杂性疾病为心脑血管疾病。
5.如权利要求1所述的评价方法,其特征在于,所述复杂性疾病为急性心肌梗塞。
6.如权利要求1所述的评价方法,其特征在于,步骤(1)中,获取所述潜在作用靶点的方法,包括以下步骤:
A)获得中药作用于复杂性疾病的芯片数据;
B)采用cMAP法得到中药作用复杂性疾病的潜在相关靶蛋白;
C)采用EoR法得到中药作用复杂性疾病的潜在相关基因;
D)整合所述潜在相关靶蛋白和潜在相关基因,得到中药作用复杂性疾病的潜在作用靶点。
7.如权利要求6所述的评价方法,其特征在于,所述cMAP法,包括以下步骤:
(a)计算探针的FC值,当探针的FC绝对值大于1.2,且p<0.05时,定义为显著差异表达的探针;
(b)根据探针改变的幅度大小,挑选中药芯片数据中上调或下调幅度最大的显著差异表达的探针与cMAP数据库中的芯片数据进行相似度匹配,找到与中药相关性最强的化合物;
(c)根据STITCH4.0数据库中化合物-蛋白质相互作用的数据信息,得到潜在作用蛋白靶点。
8.如权利要求6所述的评价方法,其特征在于,所述EoR法,包括以下步骤:
A)根据芯片数据计算节点的回调状态;
式(6)中,RL′表示回调状态;LogE为各组的基因表达水平取log2对数后的均值;Edrug为中药给药组的芯片数据;Esham为假手术组的芯片数据;Emodel为模型对照组的芯片数据;
B)计算节点的回复调控效率
EoR=100%-|100%-RL′| (7)
式(7)中,EoR表示为回复调控效率,是用于表征节点回调效率的指标;RL′表示回调状态;EoR值大于0的节点为所述潜在相关基因。
9.如权利要求1所述的评价方法,其特征在于,所述“失和”生物分子网络的构建,包括以下步骤:
A)对潜在相关靶基因和潜在相关靶蛋白之间的对应关系进行合并,去重;
B)应用DAVID Bioinformatics Resources 6.7软件对靶蛋白列表进行通路富集分析,找到所述复杂性疾病相关的生物学信号通路,定义p<0.05为显著富集通路;
C)以KEGG数据库中所述显著富集通路为基础,构建所述“失和”生物分子网络。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710308991.3A CN107292126B (zh) | 2017-05-04 | 2017-05-04 | 一种中药对复杂性疾病所致“失和”网络整合调节作用的定量评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710308991.3A CN107292126B (zh) | 2017-05-04 | 2017-05-04 | 一种中药对复杂性疾病所致“失和”网络整合调节作用的定量评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107292126A CN107292126A (zh) | 2017-10-24 |
CN107292126B true CN107292126B (zh) | 2019-12-24 |
Family
ID=60094330
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710308991.3A Active CN107292126B (zh) | 2017-05-04 | 2017-05-04 | 一种中药对复杂性疾病所致“失和”网络整合调节作用的定量评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107292126B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109411033B (zh) * | 2018-11-05 | 2021-08-31 | 杭州师范大学 | 一种基于复杂网络的药物疗效筛选方法 |
CN109738552A (zh) * | 2019-03-06 | 2019-05-10 | 湖南中医药大学第一附属医院((中医临床研究所)) | 一种筛选山茱萸降糖活性成分的方法 |
CN111785333B (zh) * | 2020-03-20 | 2022-02-25 | 中国医学科学院阜外医院深圳医院 | 基于生物学网络数据的药物作用靶点筛选方法、装置、电子设备及存储介质 |
CN113053455A (zh) * | 2021-03-05 | 2021-06-29 | 北京中医药大学 | 一种潜在味效关联关键质量属性辨识方法在小儿消积止咳口服液中的应用 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102646168A (zh) * | 2012-04-16 | 2012-08-22 | 南京大学 | 适用中药方剂网络的基于共近邻相似三角形凝聚的层次重叠社区发现方法 |
CN103678954A (zh) * | 2013-12-11 | 2014-03-26 | 深圳先进技术研究院 | 一种由生物芯片数据构建多类别特异表达分子集及类别网的方法及其应用和评价方法 |
CN104666771A (zh) * | 2015-03-27 | 2015-06-03 | 黑龙江中医药大学 | 一种治疗抑郁症的中药胶囊及制备工艺与检测分析方法 |
CN104978474A (zh) * | 2014-04-11 | 2015-10-14 | 中国中医科学院中医临床基础医学研究所 | 一种基于分子网络的药效评价方法及系统 |
CN105641694A (zh) * | 2008-12-19 | 2016-06-08 | H.隆德贝克有限公司 | 用于治疗精神和行为障碍的对Vps10p-结构域受体家族的调节 |
CN105956413A (zh) * | 2016-04-27 | 2016-09-21 | 王�忠 | 识别生物分子网络中关键模块或关键节点的方法 |
-
2017
- 2017-05-04 CN CN201710308991.3A patent/CN107292126B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105641694A (zh) * | 2008-12-19 | 2016-06-08 | H.隆德贝克有限公司 | 用于治疗精神和行为障碍的对Vps10p-结构域受体家族的调节 |
CN102646168A (zh) * | 2012-04-16 | 2012-08-22 | 南京大学 | 适用中药方剂网络的基于共近邻相似三角形凝聚的层次重叠社区发现方法 |
CN103678954A (zh) * | 2013-12-11 | 2014-03-26 | 深圳先进技术研究院 | 一种由生物芯片数据构建多类别特异表达分子集及类别网的方法及其应用和评价方法 |
CN104978474A (zh) * | 2014-04-11 | 2015-10-14 | 中国中医科学院中医临床基础医学研究所 | 一种基于分子网络的药效评价方法及系统 |
CN104666771A (zh) * | 2015-03-27 | 2015-06-03 | 黑龙江中医药大学 | 一种治疗抑郁症的中药胶囊及制备工艺与检测分析方法 |
CN105956413A (zh) * | 2016-04-27 | 2016-09-21 | 王�忠 | 识别生物分子网络中关键模块或关键节点的方法 |
Non-Patent Citations (1)
Title |
---|
四逆汤加减治疗冠心病心绞痛的疗效评价与Meta分析;庄灿;《广西中医药大学学报》;20130915;第16卷(第3期);第132-135页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107292126A (zh) | 2017-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107292126B (zh) | 一种中药对复杂性疾病所致“失和”网络整合调节作用的定量评价方法 | |
Wang et al. | Network pharmacology‑based prediction of the active compounds and mechanism of Buyang Huanwu Decoction for ischemic stroke | |
CN110448636A (zh) | 一种具有抗抑郁功效的百合地黄汤标准煎液制备方法及应用 | |
CN109738552A (zh) | 一种筛选山茱萸降糖活性成分的方法 | |
Cao et al. | Atractylodesin III maintains mitochondrial function and inhibits caspase-3 activity to reverse apoptosis of cardiomyocytes in AMI rats | |
Xin et al. | Mechanism of ‘Invigorating Qi and Promoting Blood Circulation’Drug Pair Ginseng-Danshen on Treatment of Ischemic Heart Disease Based on Network Pharmacology | |
Liu et al. | Identification of PIK3CG as a hub in septic myocardial injury using network pharmacology and weighted gene co‐expression network analysis | |
Zhang et al. | Rheum officinale Baill. Treats zebrafish embryo thrombosis by regulating NOS3 expression in the arginine biosynthesis pathway | |
Yang et al. | Network Pharmacology‐Based Prediction and Verification of the Potential Mechanisms of He’s Yangchao Formula against Diminished Ovarian Reserve | |
Ye et al. | Renshen Baidu powder protects ulcerative colitis via inhibiting the PI3K/Akt/NF-κB signaling pathway | |
Fan et al. | Mechanism of xiezhuo huayu yiqi tongluo formula in the treatment of uric acid nephropathy based on network pharmacology, molecular docking, and in vivo experiments | |
CN114783529A (zh) | 木犀草素-孤独症相关作用基因靶点及其筛选方法 | |
Zhang et al. | Optimized New Shengmai Powder ameliorates myocardial fibrosis in rats with heart failure by inhibition of the MAPK signaling pathway | |
CN111888394A (zh) | 痰热清在制备治疗新型冠状病毒肺炎药物中的应用 | |
Yu et al. | Investigation of the pharmacological effect and mechanism of mountain-cultivated ginseng and garden ginseng in cardiovascular diseases based on network pharmacology and zebrafish experiments | |
Zheng et al. | Bao Yuan decoction alleviates fatigue by restraining inflammation and oxidative stress via the AMPK/CRY2/PER1 signaling pathway | |
Fang et al. | Clinical observation on YiSuiShengXueGranule on treating 156 patients with β-thalassemia major and the molecular mechanism study | |
Bhadresha et al. | Ayurgenomics: A brief note on ayurveda and their cross kingdom genomics | |
CN1919239A (zh) | 一种治疗心脑血管疾病的中药组合物 | |
Ye et al. | Network pharmacology and experimental verification of the potential mechanism of Er-Xian decoction in aplastic anemia | |
Shi et al. | Qi-Po-Sheng-Mai granule ameliorates Ach-CaCl2-induced atrial fibrillation by regulating calcium homeostasis in cardiomyocytes | |
Jiang et al. | The acute toxic effect of Chinese medicine Fuzi is exacerbated in kidney yang deficiency mice due to metabolic difference | |
Jiang et al. | Network pharmacology combined with lipidomics to reveal the regulatory effects and mechanisms of Kangzao granules in the hypothalamus of rats with central precocious puberty | |
Wang et al. | Analysis of network pharmacology and experimental verification to explore the mechanism of Modified Guilu Erxian Glue in the treatment of low immune function | |
CN1919247A (zh) | 一种治疗心脑血管病的中药 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |