Susceptible gene of stasis-stagnation constitution from genome-wide association study related to cardiovascular disturbance and possible regulated traditional Chinese medicine

BMC Complementary and Alternative Medicine

Fuu-Jen Tsai⁵  Yu-Chian Chen⁴  Rong-Hsing Chen³  Tzu-Yuan Wang³  Chwen-Tzuei Chang³  Ching-Chu Chen³  Hung-Jin Huang²  Kuo-Chin Huang¹ 

[1] School of Chinese Medicine, College of Chinese Medicine, China Medical University, Taichung 40402, Taiwan;Department of Chinese Pharmaceutical Sciences and Chinese Medicine Resources, College of Pharmacy, China Medical University, Taichung 40402, Taiwan;Division of Endocrinology and Metabolism, Department of Medicine, China Medical University Hospital, Taichung 40447, Taiwan;Department of Biotechnology, Asia University, Taichung 41354, Taiwan;Department of Medical Genetics, Medical Research and Pediatrics, China Medical University Hospital, No. 2, Yuh Der Road, Taichung, Taiwan

关键词: Molecular dynamics (MD) simulation;    Type 2 diabetes;    Traditional Chinese medicine;    Body constitution;    Genome-wide association study;    Type 2 diabetes;   

Others  :  1220034
DOI  :  10.1186/s12906-015-0761-x

 received in 2014-07-08, accepted in 2015-07-02,  发布年份 2015

【 摘 要 】

Background

This study identified susceptible loci related to the Yu-Zhi (YZ) constitution, which indicates stasis-stagnation, found in a genome-wide association study (GWAS) in patients with type 2 diabetes and possible regulated traditional Chinese medicine (TCM) using docking and molecular dynamics (MD) simulation.

Methods

Non-aboriginal Taiwanese with type 2 diabetes were recruited. Components of the YZ constitution were assessed by a self-reported questionnaire. Genome-wide SNP genotypes were obtained using the Illumina HumanHap550 platform. The world’s largest TCM database (http://tcm.cmu.edu.tw/) was employed to investigate potential compounds for PON2 interactions.

Results

The study involved 1,021 unrelated individuals with type 2 diabetes. Genotyping data were obtained from 947 of the 1,021 participants. The GWAS identified 22 susceptible single nucleotide polymorphisms on 13 regions of 11 chromosomes for the YZ constitution. Genotypic distribution showed that PON2 on chromosome 7 was most significantly associated with the risk of the YZ constitution. Docking and MD simulation indicated 13-hydroxy-(9E_11E)-octadecadienoic acid was the most stable TCM ligand.

Conclusions

Risk loci occurred in PON2, which has antioxidant properties that might protect against atherosclerosis and hyperglycemia, showing it is a susceptible gene for the YZ constitution and possible regulation by 13-hydroxy-(9E_11E)-octadecadienoic acid.

【 授权许可】

   
2015 Huang et al.

【 图 表 】

Fig. 1.

Fig. 2.

Fig. 3.

Fig. 4.

Fig. 5.

Fig. 6.

Fig. 7.

Fig. 8.

Fig. 9.

Fig. 10.

Fig. 11.

Fig. 12.

Fig. 13.

【 参考文献 】

• [1]Theories of Physical Constitutions of Traditional Chinese Medicine. 1st ed. People’s Medical Publishing House, Beijing; 2005.
• [2]Chen S, Lv F, Gao J, Lin J, Liu Z, Fu Y et al.. HLA class II polymorphisms associated with the physiologic characteristics defined by Traditional Chinese Medicine: linking modern genetics with an ancient medicine. J Altern Complement Med. 2007; 13(2):231-239.
• [3]Wang Q, Yao S. Molecular basis for cold-intolerant yang-deficient constitution of traditional Chinese medicine. Am J Chin Med. 2008; 36(5):827-834.
• [4]Wu Y, Cun Y, Dong J, Shao J, Luo S, Nie S et al.. Polymorphisms in PPARD, PPARG and APM1 associated with four types of traditional Chinese medicine constitutions. J Genet Genomics. 2010; 37(6):371-379.
• [5]Grundy SM, Benjamin IJ, Burke GL, Chait A, Eckel RH, Howard BV et al.. Diabetes and cardiovascular disease: a statement for healthcare professionals from the American Heart Association. Circulation. 1999; 100(10):1134-1146.
• [6]Sarwar N, Gao P, Seshasai SR, Gobin R, Kaptoge S, Di Angelantonio E et al.. Diabetes mellitus, fasting blood glucose concentration, and risk of vascular disease: a collaborative meta-analysis of 102 prospective studies. Lancet. 2010; 375(9733):2215-2222.
• [7]Gao ZY, Xu H, Shi DZ, Wen C, Liu BY. Analysis on outcome of 5284 patients with coronary artery disease: the role of integrative medicine. J Ethnopharmacol. 2012; 141(2):578-583.
• [8]Lei Y, Wang ZH, Zhao H, Liu JG. Study of the relationship between carotid intima-media thickness and traditional Chinese medicine syndrome of dyslipidemia. Chin J Integr Med. 2009; 15(2):112-116.
• [9]Min MA. Expression of platelet CD62p gene and leucocyte HSP70 gene in patients with blood stasis syndrome MA. Zhongguo Zhong Xi Yi Jie He Za Zhi. 2005; 25(4):307-310.
• [10]Ma XJ, Yin HJ, Chen KJ. [Investigation of gene expression profiles in patients with blood stasis syndrome]. Zhong Xi Yi Jie He Xue Bao. 2008; 6(4):355-360.
• [11]Chang Y-M, Velmurugan BK, Kuo W-W, Chen Y-S, Ho T-J, Tsai C-T et al.. Inhibitory effect of alpinate Oxyphyllae fructus extracts on Ang II-induced cardiac pathological remodeling-related pathways in H9c2 cardiomyoblast cells. BioMedicine. 2013; 3(4):148-152.
• [12]Jao C-L, Huang S-L, Hsu K-C. Angiotensin I-converting enzyme inhibitory peptides: Inhibition mode, bioavailability, and antihypertensive effects. BioMedicine. 2012; 2(4):130-136.
• [13]Lin M-C, Tsai S-Y, Wang F-Y, Liu F-H, Syu J-N, Tang F-Y. Leptin induces cell invasion and the upregulation of matrilysin in human colon cancer cells. BioMedicine. 2013; 3(4):174-180.
• [14]Huang H-J, Yu HW, Chen C-Y, Hsu C-H, Chen H-Y, Lee K-J et al.. Current developments of computer-aided drug design. J Taiwan Inst Chem Eng. 2010; 41(6):623-635.
• [15]Chen CY. A novel integrated framework and improved methodology of computer-aided drug design. Curr Top Med Chem. 2013; 13(9):965-988.
• [16]Chen KC, Sun MF, Yang SC, Chang SS, Chen HY, Tsai FJ et al.. Investigation into potent inflammation inhibitors from traditional Chinese medicine. Chem Biol Drug Des. 2011; 78(4):679-688.
• [17]Chang SS, Huang HJ, Chen CY. Two birds with one stone? Possible dual-targeting H1N1 inhibitors from traditional Chinese medicine. PLoS Comput Biol. 2011; 7(12):e1002315.
• [18]Chang SS, Huang HJ, Chen CY. High performance screening, structural and molecular dynamics analysis to identify H1 inhibitors from TCM Database@Taiwan. Mol Biosyst. 2011; 7(12):3366-3374.
• [19]Tou WI, Chang SS, Lee CC, Chen CY. Drug design for neuropathic pain regulation from traditional Chinese medicine. Sci Rep. 2013; 3:844.
• [20]Chen KC, Chang SS, Huang HJ, Lin TL, Wu YJ, Chen CY. Three-in-one agonists for PPAR-alpha, PPAR-gamma, and PPAR-delta from traditional Chinese medicine. J Biomol Struct Dyn. 2012; 30(6):662-683.
• [21]Chen KY, Chang SS, Chen CY. In silico identification of potent pancreatic triacylglycerol lipase inhibitors from traditional Chinese medicine. PLoS One. 2012; 7(9):e43932.
• [22]Chang TT, Chen KC, Chang KW, Chen HY, Tsai FJ, Sun MF et al.. In silico pharmacology suggests ginger extracts may reduce stroke risks. Mol Biosyst. 2011; 7(9):2702-2710.
• [23]Chen KC, Chang KW, Chen HY, Chen CY. Traditional Chinese medicine, a solution for reducing dual stroke risk factors at once? Mol Biosyst. 2011; 7(9):2711-2719.
• [24]Chen K-C, Yu-Chian CC. Stroke prevention by traditional Chinese medicine? A genetic algorithm, support vector machine and molecular dynamics approach. Soft Matter. 2011; 7(8):4001-4008.
• [25]Tou WI, Chen CY. In silico investigation of potential SRC kinase ligands from traditional Chinese medicine. PLoS One. 2012; 7(3):e33728.
• [26]Tsou YA, Chen KC, Lin HC, Chang SS, Chen CY. Uroporphyrinogen decarboxylase as a potential target for specific components of traditional Chinese medicine: a virtual screening and molecular dynamics study. PLoS One. 2012; 7(11): Article ID e50087
• [27]Yang SC, Chang SS, Chen CY. Identifying HER2 inhibitors from natural products database. PLoS One. 2011; 6(12): Article ID e28793
• [28]Yang SC, Chang SS, Chen HY, Chen CY. Identification of potent EGFR inhibitors from TCM Database@Taiwan. PLoS Comput Biol. 2011; 7(10): Article ID e1002189
• [29]Chen CY. TCM Database@Taiwan: the world's largest traditional Chinese medicine database for drug screening in silico. PLoS One. 2011; 6(1): Article ID e15939
• [30]Diagnosis and Classification of Diabetes Mellitus. Diabetes Care. 1997; 20(7):1183-1197.
• [31]Lin JD, Lin JS, Chen LL, Chang CH, Huang YC, Su YC. BCQs: A Body Constitution Questionnaire to assess Stasis in traditional Chinese medicine. Eur J Integr Med. 2012; 4(4):E379-E391.
• [32]Huang K-C, Chen C-C, Su Y-C, Lin J-S, Chang C-T, Wang T-Y et al.. The Relationship between Stasis-Stagnation Constitution and Peripheral Arterial Disease in Patients with Type 2 Diabetes. Evid Based Complement Alternat Med. 2014; 2014:6.
• [33]Tsai FJ, Yang CF, Chen CC, Chuang LM, Lu CH, Chang CT et al.. A genome-wide association study identifies susceptibility variants for type 2 diabetes in Han Chinese. PLoS Genet. 2010; 6(2): Article ID e1000847
• [34]Zhang Y. I-TASSER server for protein 3D structure prediction. BMC Bioinformatics. 2008; 9:40. BioMed Central Full Text
• [35]Roy A, Kucukural A, Zhang Y. I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc. 2010; 5(4):725-738.
• [36]Roy A, Yang J, Zhang Y. COFACTOR: an accurate comparative algorithm for structure-based protein function annotation. Nucleic Acids Res. 2012; 40:471-477.
• [37]Lovell SC, Davis IW, Arendall WB, de Bakker PIW, Word JM, Prisant MG et al.. Structure validation by Cα geometry: ϕ, ψ and Cβ deviation. Proteins: Struct Funct Bioinf. 2003; 50(3):437-450.
• [38]Xue B, Dunbrack RL, Williams RW, Dunker AK. Uversky VN. PONDR-FIT: a meta-predictor of intrinsically disordered amino acids. Biochim Biophys Acta. 2010; 1804(4):996-1010.
• [39]Altenhofer S, Witte I, Teiber JF, Wilgenbus P, Pautz A, Li H et al.. One enzyme, two functions: PON2 prevents mitochondrial superoxide formation and apoptosis independent from its lactonase activity. J Biol Chem. 2010; 285(32):24398-24403.
• [40]Brooks BR, Brooks CL, Mackerell AD, Nilsson L, Petrella RJ, Roux B et al.. CHARMM: the biomolecular simulation program. J Comput Chem. 2009; 30(10):1545-1614.
• [41]Pronk S, Pall S, Schulz R, Larsson P, Bjelkmar P, Apostolov R et al.. GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit. Bioinformatics. 2013; 29(7):845-854.
• [42]Priyakumar UD, MacKerell AD. Comparison of the CHARMM27, AMBER4.1 and BMS nucleic acid force fields via free energy calculations of base flipping. Abstr Pap Am Chem S. 2005; 230:U1391-U1392.
• [43]Darden T, York D, Pedersen L. Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems. J Chem Phys. 1993; 98(12):10089-10092.
• [44]Essmann U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LG. A smooth particle mesh Ewald method. J Chem Phys. 1995; 103(19):8577-8593.
• [45]Zoete V, Cuendet MA, Grosdidier A, Michielin O. SwissParam: a fast force field generation tool for small organic molecules. J Comput Chem. 2011; 32(11):2359-2368.
• [46]Tou WI, Chen CY. May disordered protein cause serious drug side effect? Drug Discov Today. 2013. doi:S1359-6446(13)00387-5 10.1016/j.drudis.2013.10.020.
• [47]Chen CY, Tou WI. How to design a drug for the disordered proteins? Drug Discov Today. 2013; 18(19–20):910-915.
• [48]Chovancova E, Pavelka A, Benes P, Strnad O, Brezovsky J, Kozlikova B et al.. CAVER 3.0: a tool for the analysis of transport pathways in dynamic protein structures. PLoS Comput Biol. 2012; 10:e1002708.
• [49]Kim BY, Jin HJ, Kim JY. Genome-wide association analysis of Sasang constitution in the Korean population. J Altern Complement Med. 2012; 18(3):262-269.
• [50]Primo-Parmo SL, Sorenson RC, Teiber J, La Du BN. The human serum paraoxonase/arylesterase gene (PON1) is one member of a multigene family. Genomics. 1996; 33(3):498-507.
• [51]She ZG, Zheng W, Wei YS, Chen HZ, Wang AB, Li HL et al.. Human paraoxonase gene cluster transgenic overexpression represses atherogenesis and promotes atherosclerotic plaque stability in ApoE-null mice. Circ Res. 2009; 104(10):1160-1168.
• [52]She ZG, Chen HZ, Yan Y, Li H, Liu DP. The human paraoxonase gene cluster as a target in the treatment of atherosclerosis. Antioxid Redox Signal. 2012; 16(6):597-632.
• [53]Ng CJ, Wadleigh DJ, Gangopadhyay A, Hama S, Grijalva VR, Navab M et al.. Paraoxonase-2 is a ubiquitously expressed protein with antioxidant properties and is capable of preventing cell-mediated oxidative modification of low density lipoprotein. J Biol Chem. 2001; 276(48):44444-44449.
• [54]Ng CJ, Hama SY, Bourquard N, Navab M, Reddy ST. Adenovirus mediated expression of human paraoxonase 2 protects against the development of atherosclerosis in apolipoprotein E-deficient mice. Mol Genet Metab. 2006; 89(4):368-373.
• [55]Ng CJ, Bourquard N, Grijalva V, Hama S, Shih DM, Navab M et al.. Paraoxonase-2 deficiency aggravates atherosclerosis in mice despite lower apolipoprotein-B-containing lipoproteins: anti-atherogenic role for paraoxonase-2. J Biol Chem. 2006; 281(40):29491-29500.
• [56]Mochizuki H, Scherer SW, Xi T, Nickle DC, Majer M, Huizenga JJ et al.. Human PON2 gene at 7q21.3: cloning, multiple mRNA forms, and missense polymorphisms in the coding sequence. Gene. 1998; 213(1–2):149-157.
• [57]Devarajan A, Bourquard N, Hama S, Navab M, Grijalva VR, Morvardi S et al.. Paraoxonase 2 deficiency alters mitochondrial function and exacerbates the development of atherosclerosis. Antioxid Redox Signal. 2011; 14(3):341-351.
• [58]Madamanchi NR, Runge MS. Mitochondrial dysfunction in atherosclerosis. Circ Res. 2007; 100(4):460-473.
• [59]Sanghera DK, Aston CE, Saha N, Kamboh MI. DNA polymorphisms in two paraoxonase genes (PON1 and PON2) are associated with the risk of coronary heart disease. Am J Hum Genet. 1998; 62(1):36-44.
• [60]Bourquard N, Ng CJ, Reddy ST. Impaired hepatic insulin signalling in PON2-deficient mice: a novel role for the PON2/apoE axis on the macrophage inflammatory response. Biochem J. 2011; 436(1):91-100.
• [61]Hegele RA, Connelly PW, Scherer SW, Hanley AJ, Harris SB, Tsui LC et al.. Paraoxonase-2 gene (PON2) G148 variant associated with elevated fasting plasma glucose in noninsulin-dependent diabetes mellitus. J Clin Endocrinol Metab. 1997; 82(10):3373-3377.
• [62]Yu F, Gao J, Zeng Y, Liu CX. Effects of adlay seed oil on blood lipids and antioxidant capacity in hyperlipidemic rats. J Sci Food Agric. 2011; 91(10):1843-1848.
• [63]Laity JH, Lee BM, Wright PE. Zinc finger proteins: new insights into structural and functional diversity. Curr Opin Struct Biol. 2001; 11(1):39-46.
• [64]Kinlay S, Libby P, Ganz P. Endothelial function and coronary artery disease. Curr Opin Lipidol. 2001; 12(4):383-389.
• [65]Sorrentino SA, Bahlmann FH, Besler C, Muller M, Schulz S, Kirchhoff N et al.. Oxidant stress impairs in vivo reendothelialization capacity of endothelial progenitor cells from patients with type 2 diabetes mellitus: restoration by the peroxisome proliferator-activated receptor-gamma agonist rosiglitazone. Circulation. 2007; 116(2):163-173.
• [66]Pavlakis E, Chiotaki R, Chalepakis G. The role of Fras1/Frem proteins in the structure and function of basement membrane. Int J Biochem Cell Biol. 2011; 43(4):487-495.
• [67]Timmer JR, Mak TW, Manova K, Anderson KV, Niswander L. Tissue morphogenesis and vascular stability require the Frem2 protein, product of the mouse myelencephalic blebs gene. Proc Natl Acad Sci U S A. 2005; 102(33):11746-11750.
• [68]Jadeja S, Smyth I, Pitera JE, Taylor MS, van Haelst M. Bentley E et al. Identification of a new gene mutated in Fraser syndrome and mouse myelencephalic blebs. Nat Genet. 2005; 37(5):520-525.
• [69]Xiao R, Xu XZ. Mechanosensitive channels: in touch with Piezo. Curr Biol. 2010; 20(21):R936-R938.
• [70]Coste B, Mathur J, Schmidt M, Earley TJ, Ranade S, Petrus MJ et al.. Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels. Science. 2010; 330(6000):55-60.
• [71]Li Z, Chim JC, Yang M, Ye J, Wong BC, Qiao L. Role of PCDH10 and its hypermethylation in human gastric cancer. Biochim Biophys Acta. 2012; 1823(2):298-305.
• [72]Williams EO, Sickles HM, Dooley AL, Palumbos S, Bisogni AJ, Lin DM. Delta Protocadherin 10 is Regulated by Activity in the Mouse Main Olfactory System. Front Neural Circuits. 2011; 5:9.