摘要
分子识别的研究具有重要的生物学和药理学意义。例如DNA的转录,激素的作用机制,免疫系统中的抗原—抗体识别,化学中的酶催化反应,以及很多药物发挥药效的机制都与之密切相关。而绝对结合自由能计算对于研究分子识别机制及其动态作用是十分重要的。自由能的正负决定了化学反应的方向,其大小决定了反应的趋势强弱。国际上有很多的实验小组投入到了自由能的计算模型的研究。Jorgensen提出的双湮灭方法被认为是一种很经典的计算绝对结合自由能的方法,并且已经得到了广泛的应用。在论文的第一部分,我们用双湮灭方法计算了谷氨酰胺结合蛋白(GlnBP)单体和底物Gln结合的绝对结合自由能。模拟得到的结果是-14.59 kcal/mol,与实验值-8.025 kcal/mol接近,并分析了双湮灭方法的优点和不足,并对该方法的改进提出了一些建议。
在过去的几十年里,分子动力学(MD)模拟已经成为研究蛋白质物理化学特性的一个重要工具。MD模拟通过给出生物大分子在原子水平上的相互作用,提供生物大分子的涨落和构象变化的详细信息。MD模拟对于研究蛋白质等生物大分子的结构与功能的关系具有重要的意义。论文的第二部分工作是选择了GlnBP单体和谷氨酰胺结合蛋白(GlnBP-Gln)复合物体系为研究对象。使用GROMOS96程序及其力场分别进行了600 ps的MD模拟。通过两个模拟所得到的轨迹,我们从体系的分子结构、自由能、铰链区柔性、口袋区相互作用及底物结合的特异性等几个方面分析了该体系的构象开合机制,得到了与实验数据是相吻合的结果。
溶剂化效应的研究是国际上异常活跃的研究热点。而静电相互作用是溶剂
北京工业大学工学硕士学位论文
一
化效应中最主要的贡献,尤其当溶质分子带电荷或高度极化时更是如此。一般
来讲,广义的PB方程可以描述蛋白质溶液体系的静电效应。本课题组以前已
把用有限差分方法求解PB方程的方法与随机动力学相结合,得到FDSD模拟
程序,FDSD模拟方法在动力学模拟过程中把大量的水分子代之以连续介质。该
程序在较小的蛋白体系中是较为成功的。在论文的最后一部分,我们选择了R6
态胰岛素六聚体这个较大的蛋白为研究对象,把FDSD模拟结果与MD模拟和
通常的真空SD模拟结果作了比较。结果表明:FDSD模拟方法对于研究分子
的结构和动力学性质均比真空SD模拟有很明显的改善,与MD模拟结果接近。
Molecular recognition is of central importance in biology and pharmacology. It underlies the action of hormones, the control of DNA transcription, the recognition of Antigens-Antibody in the immune system, the catalysis of chemical reactions by enzymes and the actions of many drugs. Calculations of absolute binding free energy play a very important role in the study of the mechanism and dynamics of the molecular recognition. The orientation of chemical reactions depends on the sign of binding free energy and the tendency depending on its corresponding value. Many groups have invested ingenuity and effort in the development of such models, double annihilation method developed by Jorgensen was proposed as a classical way of computing the absolute binding free energy of a complex in solution, which has been applied in a number of studies. In the first part of this thesis, The binding free energy between GlnBP monomer and its ligand Gin was computed by using the classical double annihilation method. It is found that the result (-14.59 kcal/mol) agrees well with experimental data (-8.025 kcal/mol). The advantage and disadvantage of this methods are analyzed.
During the past several decades, molecular dynamics (MD) simulation of proteins has become a widely used tool to deepen our understanding of the molecules. Computer simulation techniques can be used to understand the properties of a molecular system in terms of interactions at the atomic level. The results of MD simulations are very useful to understand the relation between the macromolecular
conformational changes and its biological function. In the second part of this thesis, we have chosen the protein monomer GlnBP and complex GlnBP-Gln as our study systems. The MD simulations were carried out with the GROMOS96 package and its force field for 600 ps. The MD simulation data were analyzed in terms of the overall structure, the binding free energy, the flexibility hi the region around hinges, the interactions in the binding site, pocket ligand binding specificity, and so on. The results obtained from the biochemistry experiment can be correctly identified by our MD simulations.
The study on solvation effect is the focus of the international scientific research fields. The electrostatic interaction is the most important part hi the solvation effect, especially when the solute is charged or highly polarized. Commonly, the general Poisson-Boltzmannn (PB) equation can well describe the electrostatic property of a protein in solution. We use a finite difference method to solve PB equation and incorporate the solution into SD simulation, which is called as FDSD simulation procedure and developed by our group before. The solvent is considered as a continuum in FDSD simulation procedure, which gained success in simulation of smaller proteins. In the last part of this thesis, the R-state human insulin hexamer has been selected as a target. A comparison has been made among the results obtained with FDSD simulation, MD simulation and SD simulation. The results show that the FDSD simulation has obvious improvement over the SD simulation for structure and dynamics properties, and has a good agreement with that of the MD simulation.
引文
1 Goodfellow J M, Moss D S, Computer Modeling of Biomolecular Process. New York: Ellis Horwood, 1992.
2 Warshel A, Computer Modeling of Chemical Reactions in Enzymes and Solutions. New York: Jonh Wiley & Sons, 1991.
3 Bash P A, Singh U C, Brown F K, Langridge R, and Kollman P A, Calculation of the Relative Change in Binding Free Energy of a Protein-inhibitor Complex. Science, 1987, 235:574-576.
4 Ran B G, Tilton R F, and Singh U C, Free Energy Perturbation Studies on inhibitor Binding to HIV-1 Proteinase. J. Am. Chem. Soc., 1992, 114: 4447-4452.
5 Warshel A, Sussman F and Hwang J K, Evaluation of Catalytic Free Energies in Genetically Modified Protein. J. Mol. Biol., 1988, 201: 139-159.
6 B Z Lu, W Z Chen, C X Wang and X J Xu, Protein Molecular Dynamics with Electrostatic Force Entirely Determined by a Single Poisson-Boltzmann Calculation. protein, 2002, 48,497-504.
7 Van Gunsteren W F, Protein Eng. 1988, 2: 5.
8 Jorgensen W L, Acc. Chem. Res. 1989, 22: 184.
9 Kollman, P Chem. Rev. 1993, 93: 2395.
10 Williams, D H & Cox, J P L & Doig, A J et al. or. Am. Chem. Soc. 1991, 113: 7020.
11 Discover User Guide, part1. Versions 2.9.7&95.0/3.00. Biosym/MSI. 1995,6-2.
12 Jaroslav Vacek, Peter A Kollman, J. Phys. Chem. A 1999, 103: 10015-10020.
13 Gilson M K, Given J A, Bush B L, McCammon J A, Biophys J. 1997, 72: 1047-1069.
14 Jan Hermans and Lu Wang, Inclusion of Translational and Rotational Freedom in Theoretical Estimates of Free Energies of Binding .Application to a complex of Benzene and Mutant T4
Lysozyme. J. Am. Chem. Soc. 1997, 119: 2707-2714.
15 王存新,宋伟,陈慰祖,蛋白质与极性配体复合物的绝对结合自由能计算,生物化学与生物物理进展,2003,30:143-146.
16 Berendsen H J C, Postma J P M, van Gunsteren W F, et al. Intermolecular Forces. Dordrecht:Reidel, 1981. 331.
17 Kollman P A, Massova I, Reyes C M, Kuhn B, Huo S, et al. Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. Acc Chem Res. 2001, 33: 889-897.
18 Smith K C, Honig B, Evaluation of the conformational free energies of loops in proteins, Proteins Struct Funct Genet. 1994, 18:119-132.
19 Srinivasan J, Cheatham T E, Cieplak P, Kollman P A and Case D A, Continuum solvent studies of the stability of DNA, RNA, and phosphoramidate—DNA helices. JAm Chem Soc.1998, 120: 9401-9409.
20 Yang A S, Hitz B, and Honig B, Free energy determinants of secondary structure formation.Ⅲ. Beta-turns and their role in protein folding. J Mol Biol. 1996, 259: 873-882.
21 Yang A S, Honig B. Free energy determinants of secondary structure formation. Ⅱ.Antiparallel beta-sheets. J Mol Biol. 1995, 252: 366-376.
22 Honig B, Calculating total electrostatic energies with the nonlinear Poisson-Boltzmann equation. J Phys Chem. 1990, 94: 7684-7692.
23 Sharp K A, Honig B, Electrostatic interactions in macromolecules—theory and applications. Annu Rev Biophys Biophys. Chem. 1990, 19: 301-332.
24 Lee B and Richards F M. The interpretation of protein structure: Estimation of static accessibility. J Mol Biol. 1971, 55: 379-400.
25 Aqvist J; Medina C and Samuelsson J E, Protein Eng. 1994, 7: 385.
26 S Weiner, P A Kollman, D T Nguyen and D A Case, ar. Comp. Chem., 1986, 7: 230.
27 W L Jorgenson and J Tirado-Rives, J. Am. Chem. Soc., 1988, 110: 1657.
28 B R Brooks, R E Bruccoleri, B D Olafson, S Swaminathan and M Karplus, J. Comp. Chem.,1983, 4: 187.
29 W F van Gunsteren, S R Billeter A, A Eising, P H H(?)nenberger, P Kr(?)ger, A E Mark, W R P Scott and I P Tironi. Biomolecular simulation: the GROMOS96 manual and user guide.Z(?)rich, Switzerland: vdf Hochschuverlag, 1996.
30 P A Madden and M Wilson, Chem. Soc. Rev., 1996, 25: 339.
31 M A Van der Hoef and P A Madden, Mol. Phys. 1998, 94: 417.
32 A J Stone, The theory of intermolecular forces, Oxford, UK: Clarendon, 1996.
33 Pangali C, et al., J. Chem. Phys., 1979, 71: 2982-2990.
34 C X Wang, W Z Chen, V Tran. and R. Douillard, Analysis of interfacial water structure and dynamics in α-maltose solution by molecular dynamics simulation. Chemical Physics Letters, 1996, 251:268-274.
35 Turq P. et al., J. Chem. Phys., 1977,66:3039-3044.
36 R Kubo, Rep. Progr. Phys., 1966, 29: 255.
37 卢本卓,分子模拟中的静电计算方法及蛋白质折叠的研究,博士论文,2002.
38 Ames G F L, Bacterial periplasmic transport system: structure, mechanism, and evolution.Annu. Rev. Biochem. 1986, 55: 397-425.
39 Weiner J H, & Heppel L A, A binding protein for glutamine and its relation to active transport in Escherichia coli. J. Biol. Chem. 1971, 246: 6933-6491.
40 Weiner J H, Furlong C E & Heppel L A, J. Biol. Chem. 1971, 246: 715-717.
41 Tjandra N, Simplaceanu V, Cottam P F & Ho C., Multidimensional ~1H and ~(15)N NMR investigation of glutamine-binding protein of Escherichia coli. J. Biomol. NMR, 1992, 2:149-160.
42 Chwan-Deng Hsiao, Yuh-Ju Sun, John Rose and Bi-Cheng Wang, J. Mol. Biol. 1996, 262:
225-242.
43 Yuh-Ju Sun, John Rose, Bi-Cheng Wang and Chwan-Deng Hsiao, J. Mol. Biol. 1998, 278:219-229.
44 Oh B -H, Pandit J, Kang C H, Nikaido K, Gokcen S, Ames, G F-L & Kim S-H, Three-dimensional structure of the periplasmic lysine/arginine/ornithine-binding protein with and without a ligand. J. Biol. Chem. 1993,268: 11348-11355.
45 Sharff A J, Rodseth L E, Spurlino J C & Quiocho F A, Crystallographic evidence of a large domains of the maltodextrin binding protein involved in active transport and chemotaxis.Biochemistry, 1992,31: 10657-10633.
46 Shen Q, Simplaceanu V, Cottam P F, & Ho C, Proton nuclear magnetic resonance studies on glutamine-binding protein from Escherichia coli:formation of intermolecular and intramolecular hydrogen bonds upon ligandbinding. J. Mol. Biol. 1989, 210:849-857.
47 Gilson M K, Given J A, Bush B L, McCammon J A, Biophys J. 1997, 72: 1047-1069.
48 Hermans J, Yun R H, Leech J, Cavanaugh D. Department of Biochemistry and Biophysics,University of North Carolina. SIGMA Documentation. 1999.
49 Ryckaert J P, Ciccotti G, Berendsen H J C, Numerical integration of the cartesian equation of motion of a system with constraints: molecular dynamics of n-alkanes. J Comput Phys. 1977, 23: 327-341.
50 Ames., G. F.-L. Bacterial periplasmic transport system: structure, mechanism, and evolution. Annu. Rev. Biochem. 1986,55: 397-425.
51 Karplus M, Petsko G A Molecular dynamics Simulations in biology. Nature, 1990, 347: 631-639.
52 Smith L J, Dobson C M, et al. Side-chain conformational disorder in a moltenglobule:molecular dynamics simulations of the A-state of human α-lactaklbumin, J. Mol. Biol, 1999,286: 1567-1580.
53 Van Gunsteren W F, Billeter S R, Eising A. A, and h(?)nenberger P. H., Kr(?)ger P., Mark A E, Scott W R P & Tironi I G. Biomolecular Simulation: The GROMOS96 manualAnd user guide. Vdf Hochschulverlag AG an der ETH Z(?)rich, Z(?)rich, Switzerlan, 1996
54 Berendsen H J C, Postma J P M, van gunsteren V F, et al. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81(8): 3684-3690.
55 Nanzer A P, Poulsen F M, Van Gunsteren W F, et al. A Reassessment of the structure of chymotrypsin inhibitor 2 (CI-2) using time-averaged NMR restraints. Biochemistry, 1994, 33:14503-14511.
56 Nanzer A P, Parametrisation of time averaged distance restrains in MD simulations. J. Biol. NMR, 1995, 6:313-320.
57 B Lee, F M Richards, The interpretation of protein structures: estimation of static accessibility. J. Mol. Biology. 1971, 55:379-400.
58 D Van der Spoel, A R Van Buuren, E Apol, et al. Biomolecular4 simulations: the GROMACS user manual. Nijenborgh 4, 9747 AG Groningenm, Netherlands, 1991
59 龚庆国,吴季辉,刘勤,施蕴渝,崔大敷,徐来根等,乙肝病毒表面抗原Pre-S1的N端片断溶液构象研究.中国科学 1999,29
60 Hamaguhi K, The Protein Molecule, Conformation, Stability and Folding. Tokyo: Japan Scientific Societies Press, 1992, 71-84.
61 M J Adams, T L Blundell et al. Structure of Rhombohedral 2 Zinc Insulin Crystals. Nature,1969, 224: 491-495.
62 刘新奇,金雷,张英,王大成,不同聚合态胰岛素的晶体学研究:重组突变体 B9(Ser—Asp)人胰岛素的系列结晶及其初步晶体学研究.科学通报,44(4):403—406.
63 北京胰岛素结构研究组.胰岛素晶体结构的研究:1.8埃分辨率的胰岛素分子。中国科学 1974,6:591-606.
64 Blundell T L, Cutfield J F, et al. Atomic positions in rhombohedral 2-zinc insulin crystals.
Nature, 1971, 231: 591-606.
65 B Honig, K Sharp and A S Yang, J. Phys. Chem., 1993, 97: 1101.
66 B Honig and A Nicholls, Science, 1995, 268:1144.
67 S C Harvey, Proteins, 1989, 5: 78.
68 M E Davis and J A McCammon, Chem. Rev., 1990, 90: 509.
69 K A Sharp and B Honig, Annu. Rev. Biophys. Biochem., 1990, 19: 301.
70 Warshel A, & Aqvist J, Electrostatic energy and macromolecular function. Annu. Rev. Biophys. Chem., 1991, 20: 267-298.
71 Warshel A, & Russell S T, Q. Rev. Biophys., 1984, 17: 283-297.
72 Russell S T & Warshel A, J. Mol. Biol., 1985, 185: 389-406.
73 Warshel A & Levitt M, Science, 1995, 268: 1144-1146.
74 Xu Y W, Wang C X. & Shi Y Y, J. Comp. Chem., 1992, 13:1109.
75 Warwicker J & Watson H C, J. Mol. Biol., 1982, 157: 671.
76 McCammon J A, et al., Chem. Rev., 1990, 90: 509.
77 Xiang Z X, Shi Y Y & Xu Y W,J. Comp. Chem., 1995, 16: 200.
78 Cheng W, Wang C X, Chen W Z, Xu Y W & Shi Y Y, Eur. Biophys. J., 1998, 27: 105.
79 王维宴,卢本卓,陈慰祖,王存新,胰岛素六聚体在水溶液中的稳定性研究,化学学报,2002,60:2129—2134。
80 王维宴,卢本卓,陈慰祖,王存新,胰岛素六聚体在水溶液中的构象柔性研究,生物物理学报,2002,19:35-40.
81 W Kabsch, C Sandder. Dictionary of protein secondary structure pattern recongnition of hydrogen bonded and geometrical features. Biopolymers, 1983, 22:2577-2637.
82 C.X.. Wang, A. R. Bizzarri, Y. W. Xu and S. Cannistraro, Chem. Phys., 1994, 183:155.
83 C N Schutz and A Warshel. Proteins, 2001, 44:400.
84 王维宴.胰岛素六聚体的分子动力学模拟研究.硕士论文.2002.
85 王志珍,梁栋材.生物化学杂志,1985,1:9-22.
86 张洪涛,陈劲秋,胡美浩.生物化学杂志,1995,11:354-358.
87 Qing Xin Hua, S E Shoelson et al. Receptor binding redefined by a structural switch in a mutant human insulin. Nature, 1991, 354:238-241.
88 梁栋材,宋浪舟,万柱礼,常文瑞.中国科学(B辑),1991,7:715-723.