化学工程联合国家重点实验室,清华大学化学工程系,北京 100084
彭雪(1994—),女,博士研究生,px17@mails.tsinghua.edu.cn
卢滇楠(1978—),男,博士,副教授,ludiannan@tsinghua.edu.cn
收稿:2019-07-30,
修回:2019-09-03,
网络首发:2020-08-14,
纸质出版:2020-02-05
移动端阅览
彭雪, 芦琛璘, 卢滇楠. 氧气和一氧化碳在人血红蛋白迁移过程研究[J]. 化工学报, 2020, 71(2): 724-735
PENG Xue, LU Chenlin, LU Diannan. Investigation on migration process of oxygen and carbon monoxide in human hemoglobin[J]. CIESC Journal, 2020, 71(2): 724-735
彭雪, 芦琛璘, 卢滇楠. 氧气和一氧化碳在人血红蛋白迁移过程研究[J]. 化工学报, 2020, 71(2): 724-735 DOI: 10.11949/0438-1157.20190862.
PENG Xue, LU Chenlin, LU Diannan. Investigation on migration process of oxygen and carbon monoxide in human hemoglobin[J]. CIESC Journal, 2020, 71(2): 724-735 DOI: 10.11949/0438-1157.20190862.
为了揭示CO和O
2
竞争性结合人血红蛋白血红素位点的机制及其与人血红蛋白结构转换之间的关系,本文采用全原子分子动力学模拟(MD)结合马尔科夫状态模型(MSMs)研究氧气(O
2
)和一氧化碳(CO)分子从水溶液迁移进入人血红蛋白四聚体α链和β链的全过程。分子动力学模拟揭示了O
2
和CO结合α链和β链的稳态结合位点和瞬态结合位点、迁移通道以及α链的结构变化。结果显示,分子模拟不仅仅能够再现全部实验中所观察到的离散Xe结合位点和分子扩散通道,而且揭示了实验中无法观测的瞬态结合位点和多重气体迁移途径。上述结果表明人血红蛋白因其结构柔性所形成的瞬态通道对于气体分子迁移过程的重要性。除此之外,利用MSM和过渡路径理论(TPT)构建了人血红蛋白α链结构变化与气体分子迁移之间的关系,阐释了血红蛋白中影响气体迁移的关键结构及其微观机制。
To reveal the relationship between CO and O
2
competitive binding to the heme site of human hemoglobin (HbA) and its relationship with HbA structural transition
we used all-atom molecular dynamics simulations (MD) combined with Markov state models (MSMs) to study the migration of O
2
and CO molecules from the aqueous solution into the heme site of α and β chains of HbA. Molecular dy
namics simulations revealed the steady binding sites and instantaneous binding sites to O
2
and CO of α and β chains of HbA
migration pathways and structural changes of HbA. It is shown that all discrete Xe binding sites and O
2
migration pathways observed in experiments can be reproduced by molecular dynamics simulations. Molecular dynamics simulations revealed structural changes in the homeostatic binding sites and transient binding sites
migration channels
and alpha chains of O
2
and CO binding α
and β chains. Above results demonstrate the importance of instantaneous migration pathways formed due to the flexibility of HbA for gas migration. In addition
MSMs and transition path theory(TPT) were used to investigate the conformational transitions among different microstates in α chain
which provides the mechanism between gases migration behavior and HbA conformational transitions.
Feitelson J , McLendon G . Migration of small molecules through the structure of hemoglobin: evidence for gating in a protein electrontransfer reaction [J ] . Bioc ., 1991 , 30 : 5051 - 5055 .
Leitner D M . Energy flow in proteins [J ] . Annu. Rev. Phys. Chem ., 2008 , 59 : 233 - 259 .
Savino C , Miele A E , Draghi F , et al . Pattern of cavities in globins: the case of human hemoglobin [J ] . Biopolymers , 2009 , 91 : 1097 - 1107 .
Park S Y , Yokoyama T , Shibayama N , et al . 1.25 Å resolution crystal structures of human haemoglobin inthe oxy, deoxy and carbonmonoxy forms [J ] . J. Mol. Biol ., 2006 , 360 : 690 - 701 .
Schotte F , Lim M , Jackson T , et al . Watching a protein as it functions with 150-ps time-resolved X-ray crystallography [J ] . Science , 2003 , 300 : 1944 - 1947 .
Schotte F , Soman J , Olson J S , et al . Picosecond time-resolved X-ray crystallography: probing protein function in real time [J ] . Struct. Biol ., 2004 , 147 : 235 - 246 .
Schmidt M , Nienhaus K , Pahl R , et al . Ligand migration pathway and protein dynamics in myoglobin: a time-resolved crystallographic study on L29W MbCO [J ] . Proc. Natl. Acad ., 2005 , 102 : 11704 - 11709 .
Perutz M F , Mathews F S . An X-ray study of azide methaemoglobin [J ] . Biol ., 1966 , 21 : 199 - 202 .
Lucas M F , Guallar V . An atomistic view on human hemoglobin carbon monoxide migration processes [J ] . Biophys ., 2012 , 102 : 887 - 896 .
Lepeshkevich S V , Biziuk S A , Lemeza A M . The kinetics of molecular oxygen migration in the isolated alpha chains of human hemoglobin as revealed by molecular dynamics simulations and laser kinetic spectroscopy [J ] . Biochimica et Biophysica Acta , 2011 , 1814 : 1279 - 1288 .
Lepeshkevich S V , Gilevich S N , Parkhats M V , et al . Molecular oxygen migration through the xenon docking sites of human hemoglobin in the R-state [J ] . Biochimica et Biophysica Acta (BBA) - Proteins and Proteomics , 2016 , 1864 ( 9 ): 1110 - 1121 .
Hummer G , Schotte F , Anfinrud P A , et al . Unveiling functional protein motions with picosecond X-ray crystallography and molecular dynamics simulations [J ] . Proceedings of the National Academy of Sciences of the United States of America , 2004 , 101 ( 43 ): 15330 - 15334 .
Ruscio J Z , Kumar D , Shukla M . Atomic level computational identification of ligand migration pathways between solvent and binding site in myoglobin [J ] . Proc. Natl ., 2008 , 105 : 9204 - 9209 .
Elber R . Ligand diffusion in globins simulations versu experiment [J ] . Biol ., 2010 , 20 : 162 - 167 .
Ostermann A , Waschipky R , Parak F G , et al . Ligand binding and conformational motions in myoglobin [J ] . Nature , 2000 , 404 ( 6774 ): 205 - 208 .
Cordone L , Cottone G , Giuffrida S , et al . Thermal evolution of the CO stretching band in carboxy-myoglobin in the light of neutron scattering and molecular dynamics simulations [J ] . Chemical Physics , 2008 , 345 ( 2 ): 275 - 282 .
Yu T Q , Lapelosa M , Vandeneijnden E , et al . Full kinetics of CO entry, internal diffusion, and exit in myoglobin from transition-path theory simulations [J ] . Journal of the American Chemical Society , 2015 , 137 ( 8 ): 3041 - 3050 .
Alessandro B , Massimiliano B , Michele P . Wiley interdisciplinary reviews: computational molecular [J ] . Journal of Computational Chemistry , 2011 , 8 : 826 - 843 .
Leitner D M . Energy flow in proteins [J ] . Annu. Rev. Phys. Chem ., 2008 , 59 : 233 - 259 .
Elber R , Karplus M . Enhanced sampling in molecular dynamics: use of the time-dependent hartree approximation for a simulation of carbon monoxide diffusion through myoglobin [J ] . Journal of the American Chemical Society , 1990 , 112 : 9161 - 9175 .
Meller J , Elber R . Computer simulations of carbon monoxide photodissociation in myoglobin: structural interpretation of the Bstates [J ] . Biophys ., 1998 , 74 : 789 - 802 .
Bossa C , Anselmi M , Roccatano D , et al . Extended molecular dynamics simulation of thecarbon monoxide migration in sperm whale myoglobin [J ] . Biophys ., 2004 , 86 : 3855 - 3862 .
Cohen J , Arkhipov A , Braun R , et al . Imaging the migration pathways for O 2 , CO, NO, and Xe inside myoglobin [J ] . Biophysical Journal , 2006 , 91 ( 5 ): 1844 - 1857 .
Sancho D D , Kubas A , Wang P H , et al . Identification of mutational hot spots for substrate diffusion: application to myoglobin [J ] . Journal of Chemical Theory & Computation , 2015 , 11 ( 4 ): 1919 - 1927 .
Shadrina M S , English A M , Peslherbe G H . Effective simulations of gas diffusion through kinetically accessible tunnels in multisubunit proteins: O 2 pathways and escape routes in T-state deoxyhemoglobin [J ] . J. Am. Chem. Soc ., 2012 , 134 : 11177 - 11184 .
Takayanagi M , Kurisaki I , Nagaoka M . Oxygen entry through multiple pathways in T-state human hemoglobin [J ] . J. Phys. Chem. B , 2013 , 117 : 6082 - 6091 .
Wang P H , Best R B , Blumberger J . Multiscale simulation reveals multiple pathways for H 2 and O 2 transport in a [NiFe ] -hydrogenase [J ] . Journal of the American Chemical Society , 2013 , 10 : 3548 - 3556 .
Bustamante J P , Szretter M E , Sued M , et al . A quantitative model for oxygen uptake and release in a family of hemeproteins [J ] . Bioinformatics , 2016 , 32 ( 12 ): 1805 - 1813 .
Liang J , Edelsbrunner H , Fu P , et al . Analytical shape computation of macromolecules (Ⅱ): Inacessiblecavities in proteins [J ] . Proteins: Struct . Func. Gen ., 1998 , 33 : 18 - 29 .
Shaanan B . Structure of human oxyhaemoglobin at 2.1 Å resolution [J ] . Mol. Biol ., 1983 , 171 : 31 - 59 .
Wang P H , De S D , Best R B . Computation of rate constants for diffusion of small ligands to and from buried protein active sites [J ] . Methods in Enzymology , 2016 , 578 : 299 - 326 .
Birukou I , Maillett D H , Birukova A . Modulating distal cavities in the α and β subunits of human HbA reveals the primary ligand migration pathway [J ] . Biochemistry , 2011 , 50 ( 34 ): 7361 - 7374 .
Takayanagi M , Nagaoka M . Incipient structural and vibrational relaxation process of photolyzed carbonmonoxy myoglobin: statistical analysis by perturbation ensemble molecular dynamics method [J ] . Theor. Chem. Acc ., 2011 , 130 : 1115 - 1129 .
Estarellas M C , Seira C C , Luque F J , et al . Understanding the kinetics of ligand binding to globins with molecular dynamics simulations: the necessity of multiple state models [J ] . Drug Discovery Today Technologies , 2015 , 17 ( 70 ): 22 - 27 .
Brooks B R , Bruccoleri R E , Olafson B D , et al . CHARMM: a program for macromolecular energy, minimization, and dynamics calculations [J ] . Comp. Chem ., 1983 , 4 : 187 - 217 .
Shadrina M S , Peslherbe G H , English A M . O 2 and water migration pathways between the solvent and heme pockets of hemoglobin with open and closed conformations of the distal hisE7 [J ] . Biochemistry , 2015 , 54 ( 34 ): 5279 - 5289 .
Leitner D M , Straub J E . Proteins: Energy, Heat and Signal Flow (Computation in Chemistry) [M ] . Boca Raton , FL: CRC Press , 2009 .
Cameron A , Giovanni B . Enhanced sampling in molecular dynamics using metadynamics, replica exchange, and temperature-acceleration [J ] . Entropy , 2013 , 16 : 163 - 199 .
Berenbrink M . Evolution of vertebrate haemoglobins: histidine side chains, specific buffer value and Bohr effect [J ] . Respiratory Physiology & Neurobiology , 2006 , 154 ( 1 ): 165 - 184 .
Huber K P , Herzberg G . Molecular Spectra and Molecular Structure [M ] . US : Springer , 1979 .
0
浏览量
0
下载量
1
CSCD
关联资源
相关文章
相关作者
相关机构
京公网安备11010102001995号