Improvement of Magnesium Cation-Chloride and Acetate Anion Interaction Lennard-Jones Parameters by Experimental Osmotic Pressure: Difficulty Arises from Slow Water Exchange Rate

以作者查詢圖書館館藏

、以作者查詢臺灣博碩士

、以作者查詢全國書目

、勘誤回報

、線上人數：65

、訪客IP：18.190.207.122

姓名

高立昇(Kao-Li Sheng) 查詢紙本館藏

畢業系所

化學學系

論文名稱

(Improvement of Magnesium Cation-Chloride and Acetate Anion Interaction Lennard-Jones Parameters by Experimental Osmotic Pressure: Difficulty Arises from Slow Water Exchange Rate)

相關論文

★ 嗜甲烷菌內甲烷單氧化酵素中催化反應中心三核銅模擬分子之合成與光譜分析	★ 烷烴氧化菌及氧化酵素之純化與功能性探討
★ 以電腦模擬研究香蕉型液晶元的分子交互作用力	★ 利用時間相關的電子密度泛函理論研究反式-二苯乙烯胺的光化學行為
★ 以生物資訊法研究穩定Asparagine在左手螺旋形下的交互作用力	★ 葛蘭氏陰性菌脂質A之結構研究
★ 五苯荑衍生之苯乙炔寡聚物之合成與光物理性質研究	★ 紫質三元件系統的金屬化作用對遠端氫鍵調控的影響
★ 非鍵結作用力的理論研究： (1)質子化與氧化三元件系統遠端調控氫鍵的比較 (2)π- π與CH- π作用力的取代基效應	★ 利用時間相關的密度泛涵理論研究HBI分子及其衍生物在第一激發態的位能曲線
★ Replica-Exchange分子動態模擬法研究類澱粉胜肽25-35 嵌入膜與折疊的行為	★ 抗菌胜肽資料庫分析及利用分子動態模擬法探討抗菌胜肽Indolicidin於生物膜上的行為
★ 網頁圖形界面在分子模擬上的應用	★ 類澱粉胜肽Abeta(25-35) 序列影響該類胜肽在水-膜環境下的組態: 強調多樣性的神經毒性
★ 以分子動態模擬法研究陽離子-負電磷脂質雙層的配位網絡結構：延伸應用於膜融合機制	★ 染料敏化太陽能電池吸光性質的計算研究

檔案

[Endnote RIS 格式]

[Bibtex 格式]

[相關文章]

[文章引用]

[完整記錄]

[館藏目錄]

[檢視]

[下載]

本電子論文使用權限為同意立即開放。
已達開放權限電子全文僅授權使用者為學術研究之目的，進行個人非營利性質之檢索、閱讀、列印。
請遵守中華民國著作權法之相關規定，切勿任意重製、散佈、改作、轉貼、播送，以免觸法。

摘要(中)

在計算模擬實驗中，離子可利用於多樣的生物系統之中，因此離子在水溶液中時，各項數據準確度是非常重要的。值得一提的是鈣離子與鎂離子在細胞膜生理功能上扮演極具影響的角色。早前的古典力場分子模擬研究指出高濃度電解液中會觀察到過多離子簇生成，其原因來自於使用Lorentz-Berthelot combining rule的方法估算陽離子與陰離子的蘭納-瓊斯作用力，此估算方法使用在系統預設值來進行分子動態模擬產生出的陽、陰離子相互作用力與實際狀況是有誤差的。
此次的研究，我們採用Luo與 Roux所提出的方法來改善鎂離子與氯離子、醋酸根離子上的氧原子之間的蘭納-瓊斯作用力參數，有兩組力場參數需要進行優化，第一組為CHARMM36預設的鎂離子與氯離子、醋酸根離子上的氧原子之間的蘭納-瓊斯作用力參數，第二組則為Lim的鎂離子參數。模擬結果發現在氯化鎂溶液系統中兩種參數可以重現在高濃度下的滲透壓實驗值且鎂離子與氯離子蘭納-瓊斯作用參數改進後還能重現大範圍濃度滲透壓的結果。同樣地，利用CHARMM36力場參數的鎂離子與醋酸根離子蘭納-瓊斯作用力參數改進也可沿用上述方法，並利用修正後的參數來模擬與鎂離子鍵結的蛋白質以測試其正確性。然而我們發現到鎂離子與氯離子、醋酸根離子上的氧原子之間作用力參數優化過長的現象，可能的原因為其緩慢的水交換速率。實驗測得鎂離子的水交換率高達〜1s，這表明我們的模擬時間太短而無法達到熱平衡。

摘要(英)

An accurate potential function of ion solvation in aqueous solution is crucial for computer simulations of various biological systems. Previous studies have shown that the excess ion cluster formation in concentrated electrolyte solutions observed in classical molecular dynamics simulations (MD) is arisen from the improper cation-anion Lennard-Jones (LJ) interaction parameters, commonly approximated by Lorentz-Berthelot combining rule, which are generally default used in MD simulations.
In this study, we follow the same methodology proposed by Luo and Roux to determine the Mg2+Cl and Mg2+O atom of OAcLJ interaction parameter. Two sets of Mg2+ LJ parameters are used for optimization: default CHARMM36 and Lim’s parameters. The two sets of Mg2+Cl LJ interaction parameters are able to be optimized to reproduce experimentally-measured osmotic pressure at high concentration. Using the optimized Mg2+Cl LJ interaction parameter, the osmotic pressures at a wide range of MgCl2 concentration are well reproduced. Similarly, the Mg2+oxygen atom of OAc interaction LJ parameter based on CHARMM36 FF can be optimized. Nevertheless, the Rmin values are too large to be applied widely. We ascribe the bottleneck of Mg2+ solution simulations as well as its parameter optimization to its slow water exchange rate. The experimentally-measured water exchange rate of Mg2+ is up to ~1 s indicated that our simulations are obviously too short to reach thermal equilibrium.

關鍵字(中)

★ 滲透壓
★ 計算化學

關鍵字(英)

論文目次

摘要 I
Abstract II
誌謝…………………………………………………………………………………………...III

Contents III
List of Figures V
List of Tables VIII
Chapter 1-Introduction 1
Chapter 2- Computational Methods 5
2.1 Methods for Optimizing the Mg2+Cland Mg2+Oxygen Atom of OAc LJ Parameters 5
2.2 Simulation Methods 7
2.3 Calculations of Osmotic Pressure 10
Chapter 3-Results and Discussion 13
3.1 Performance of Mg2+-Cl- LJ parameter Before Optimization 13
3.2 Optimization of Mg2+-Cl- LJ Interaction Parameter 18
3.3 Optimization of Mg2+Oxygen Atom of cLJ Interaction Parameter 27
3.4 Performance of Optimized Mg2+- Oxygen Atom of OAc- LJ Parameter Using Default CHARMM36 Force Field 30
3.5 Optimization of Mg2+Oxygen Atom of OAc LJ Parameter in Terms of Lim’s Force Field 39
Chapter 4-Conclusions and Summary 44
References 45

參考文獻

1. (a) Megyes, T.; Grosz, T.; Radnai, T.; Bako, I.; Palinkas, G., Solvation of calcium ion in polar solvents: An X-ray diffraction and ab initio study. Journal of Physical Chemistry A 2004, 108 (35), 7261-7271; (b) Pham, V. T.; Fulton, J. L., Ion-pairing in aqueous CaCl2 and RbBr solutions: simultaneous structural refinement of XAFS and XRD data. J Chem Phys 2013, 138 (4), 044201.
2. (a) Babu, C. S.; Lim, C., Empirical force fields for biologically active divalent metal cations in water. J Phys Chem A 2006, 110 (2), 691-9; (b) Mundy, M. D. B. a. C. J., Local Aqueous Solvation Structure Around Ca2+ During Ca2+···Cl− Pair
Formation. The Journal of Physical Chemistry B 2016; (c) Saxena, A.; Garcia, A. E., Multisite ion model in concentrated solutions of divalent cations (MgCl2 and CaCl2): osmotic pressure calculations. J Phys Chem B 2015, 119 (1), 219-27; (d) Tsai, H. H.; Chang, C. M.; Lee, J. B., Multi-step formation of a hemifusion diaphragm for vesicle fusion revealed by all-atom molecular dynamics simulations. Biochim. Biophys. Acta 2014, 1838 (6), 1529-35; (e) Tsai, H. H.; Juang, W. F.; Chang, C. M.; Hou, T. Y.; Lee, J. B., Molecular mechanism of Ca(2+)-catalyzed fusion of phospholipid micelles. Biochim. Biophys. Acta 2013, 1828 (11), 2729-38; (f) Tsai, H. H.; Lai, W. X.; Lin, H. D.; Lee, J. B.; Juang, W. F.; Tseng, W. H., Molecular dynamics simulation of cation-phospholipid clustering in phospholipid bilayers: possible role in stalk formation during membrane fusion. Biochim. Biophys. Acta 2012, 1818 (11), 2742-55.
3. Draper, D. E., RNA Folding: Thermodynamic and Molecular Descriptions of the Roles of Ions. Biophysical journal 2008, 95 (12), 5489-5495.
4. Guzin, K.; Goynumer, G.; Gokdagli, F.; Turkgeldi, E.; Gunduz, G.; Kayabasoglu, F., The effect of magnesium sulfate treatment on blood biochemistry and bleeding time in patients with severe preeclampsia. J Matern Fetal Neonatal Med 2010, 23 (5), 399-402.
5. Brooks, B. R.; Bruccoleri, R. E.; Olafson, B. D.; States, D. J.; Swaminathan, S.; Karplus, M., Charmm - a Program for Macromolecular Energy, Minimization, and Dynamics Calculations. J. Comput. Chem. 1983, 4 (2), 187-217.
6. Weiner, P. K.; Kollman, P. A., Amber - Assisted Model-Building with Energy Refinement - a General Program for Modeling Molecules and Their Interactions. J. Comput. Chem. 1981, 2 (3), 287-303.
7. Scott, W. R. P.; Hünenberger, P. H.; Tironi, I. G.; Mark, A. E.; Billeter, S. R.; Fennen, J.; Torda, A. E.; Huber, T.; Krüger, P.; van Gunsteren, W. F., The GROMOS Biomolecular Simulation Program Package. The Journal of Physical Chemistry A 1999, 103 (19), 3596-3607.
8. Marcus, Y., A simple empirical model describing the thermodynamics of hydration of ions of widely varying charges, sizes, and shapes. Biophysical chemistry 1994, 51, 11-127.
9. Chitra, R.; Smith, P. E., Molecular Association in Solution: A Kirkwood−Buff Analysis of Sodium Chloride, Ammonium Sulfate, Guanidinium Chloride, Urea, and 2,2,2-Trifluoroethanol in Water. The Journal of Physical Chemistry B 2002, 106 (6), 1491-1500.
10. (a) Luo, Y.; Roux, B., Simulation of Osmotic Pressure in Concentrated Aqueous Salt Solutions. The Journal of Physical Chemistry Letters 2010, 1 (1), 183-189; (b) Chen, A. A.; Pappu, R. V., Parameters of monovalent ions in the AMBER-99 forcefield: assessment of inaccuracies and proposed improvements. J Phys Chem B 2007, 111 (41), 11884-7; (c) Auffinger, P.; Cheatham, T. E.; Vaiana, A. C., Spontaneous Formation of KCl Aggregates in Biomolecular Simulations: A Force Field Issue? J Chem Theory Comput 2007, 3 (5), 1851-9; (d) Yoo, J.; Wilson, J.; Aksimentiev, A., Improved model of hydrated calcium ion for molecular dynamics simulations using classical biomolecular force fields. Biopolymers 2016, 105 (10), 752-63.
11. Kale, L.; Skeel, R.; Bhandarkar, M.; Brunner, R.; Gursoy, A.; Krawetz, N.; Phillips, J.; Shinozaki, A.; Varadarajan, K.; Schulten, K., NAMD2: Greater scalability for parallel molecular dynamics. J. Comp. Phys. 1999, 151 (1), 283-312.
12. Dudev, T.; Lim, C., Principles Governing Mg, Ca, and Zn Binding and Selectivity in Proteins. Chemical Reviews 2003, 103 (3), 773-788.
13. Sakharov, D. V.; Lim, C., Force fields including charge transfer and local polarization effects: Application to proteins containing multi/heavy metal ions. J Comput Chem 2009, 30 (2), 191-202.
14. Robinson, R. A.; Stokes, R. H., Electrolyte solutions. Courier Corporation: 2002.
15. Marchand, S.; Roux, B., Molecular dynamics study of calbindin D9k in the apo and singly and doubly calcium-loaded states. Proteins: Structure, Function, and Bioinformatics 1998, 33 (2), 265-284.
16. Humphrey, W.; Dalke, A.; Schulten, K., VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14 (1), 33-38.
17. Klauda, J. B.; Venable, R. M.; Freites, J. A.; O′Connor, J. W.; Tobias, D. J.; Mondragon-Ramirez, C.; Vorobyov, I.; MacKerell, A. D.; Pastor, R. W., Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid Types. J. Phys. Chem. B 2010, 114 (23), 7830-7843.
18. Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L., Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79 (2), 926-935.
19. Feller, S. E.; Zhang, Y. H.; Pastor, R. W.; Brooks, B. R., Constant-Pressure Molecular-Dynamics Simulation - the Langevin Piston Method. J. Chem. Phys. 1995, 103 (11), 4613-4621.
20. Steinbach, P. J.; Brooks, B. R., New Spherical-Cutoff Methods for Long-Range Forces in Macromolecular Simulation. J. Comput. Chem. 1994, 15 (7), 667-683.
21. Ryckaert, J.-P.; Ciccotti, G.; Berendsen, H. J. C., Numerical integration of the cartesian equations of motion of a system with constraints: Molecular dynamics of n-alkanes. J. Comp. Phys. 1977, 23 (3), 327-341.
22. (a) Hewish, N. A.; Neilson, G. W.; Enderby, J. E., Environment of Ca2+ ions in aqueous solvent. Nature 1982, 297 (5862), 138-139; (b) Probst, M. M.; Radnai, T.; Heinzinger, J. K.; Bopp, P.; Rode, B. M., Molecular Dynamics and X-ray Investigation of an Aqueous CaCI, Solution. J. Phys. Chem. 1985, 89, 753-759; (c) Spångberg, D.; Hermansson, K.; Lindqvist-Reis, P.; Jalilehvand, F.; Sandström, M.; Persson, I., Model Extended X-ray Absorption Fine Structure (EXAFS) Spectra from Molecular Dynamics Data for Ca2+ and Al3+ Aqueous Solutions. The Journal of Physical Chemistry B 2000, 104 (45), 10467-10472; (d) Fulton, J. L.; Heald, S. M.; Badyal, Y. S.; Simonson, J., Understanding the effects of concentration on the solvation structure of Ca2+ in aqueous solution. I: The perspective on local structure from EXAFS and XANES. The Journal of Physical Chemistry A 2003, 107 (23), 4688-4696.
23. Brooks, B. R.; Brooks, C. L.; Mackerell, A. D.; Nilsson, L.; Petrella, R. J.; Roux, B.; Won, Y.; Archontis, G.; Bartels, C.; Boresch, S.; Caflisch, A.; Caves, L.; Cui, Q.; Dinner, A. R.; Feig, M.; Fischer, S.; Gao, J.; Hodoscek, M.; Im, W.; Kuczera, K.; Lazaridis, T.; Ma, J.; Ovchinnikov, V.; Paci, E.; Pastor, R. W.; Post, C. B.; Pu, J. Z.; Schaefer, M.; Tidor, B.; Venable, R. M.; Woodcock, H. L.; Wu, X.; Yang, W.; York, D. M.; Karplus, M., CHARMM: The biomolecular simulation program. Journal of Computational Chemistry 2009, 30 (10), 1545-1614.
24. (a) Beglov, D.; Roux, B., Finite representation of an infinite bulk system: Solvent boundary potential for computer simulations. The Journal of Chemical Physics 1994, 100 (12), 9050-9063; (b) Roux, B., Valence selectivity of the gramicidin channel: a molecular dynamics free energy perturbation study. Biophysical journal 1996, 71 (6), 3177-3185.
25. Karen M. Callahan,† Nadia N., Solvation of Magnesium Dication: Molecular Dynamics Simulation and Vibrational Spectroscopic Study of Magnesium Chloride in Aqueous Solutions J. Phys. Chem. A 2010, 114, 5141–5148
26. Todor Dudev;Carmay Lim , Principles Governing Mg, Ca, and Zn Binding and Selectivity in Proteins. Chem. Rev. 2003, 103, 773787
27. Olof Allnér Lennart Nilsson Alessandra Villa Magnesium Ion–Water Coordination and Exchange in Biomolecular Simulations J. Chem. Theory Comput. 2012, 8, 4, 1493-1502
28. Serrano, L.,Coll, M. Magnesium binding to the bacterial chemotaxis protein CheY results in large conformational changes involving its functional surface. J.Mol.Biol. 1994 238:489-495 29. Quax, T.E.F.,Altegoer, F. Structure and function of the archaeal response regulator CheY. Proc. Natl. Acad. Sci. U.S.A.2018 115: E1259-E1268 30. Andersson, M.,Malmendal Structural basis for the negative allostery between Ca(2+)- and Mg(2+)-binding in the intracellular Ca(2+)-receptor calbindin D9k. Protein Sci.1997 6: 1139-1147

指導教授

蔡惠旭

審核日期

2020-1-14

推文