1
|
Jiang M, Yu CH, Xu Z, Qin Z. Binding of Carbon Monoxide to Hemoglobin in an Oxygen Environment: Force Field Development for Molecular Dynamics. J Chem Theory Comput 2024; 20:4229-4238. [PMID: 38400860 PMCID: PMC11137813 DOI: 10.1021/acs.jctc.4c00029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2024] [Revised: 02/08/2024] [Accepted: 02/12/2024] [Indexed: 02/26/2024]
Abstract
Carbon monoxide (CO) is a byproduct of the incomplete combustion of carbon-based fuels, such as wood, coal, gasoline, or natural gas. As incomplete combustion in a fire accident or in an engine, massively produced CO leads to a serious life threat because CO competes with oxygen (O2) binding to hemoglobin and makes people suffer from hypoxia. Although there is hyperbaric O2 therapy for patients with CO poisoning, the nanoscale mechanism of CO dissociation in the O2-rich environment is not completely understood. In this study, we construct the classical force field parameters compatible with the CHARMM for simulating the coordination interactions between hemoglobin, CO, and O2, and use the force field to reveal the impact of O2 on the binding strength between hemoglobin and CO. Density functional theory and Car-Parrinello molecular dynamics simulations are used to obtain the bond energy and equilibrium geometry, and we used machine learning enabled via a feedforward neural network model to obtain the classical force field parameters. We used steered molecular dynamics simulations with a force field to characterize the mechanical strength of the hemoglobin-CO bond before rupture under different simulated O2-rich environments. The results show that as O2 approaches the Fe2+ of heme at a distance smaller than ∼2.8 Å, the coordination bond between CO and Fe2+ is reduced to 50% bond strength in terms of the peak force observed in the rupture process. This weakening effect is also shown by the free energy landscape measured by our metadynamics simulation. Our work suggests that the O2-rich environment around the hemoglobin-CO bond effectively weakens the bonding, so that designing of O2 delivery vector to the site is helpful for alleviating CO binding, which may shed light on de novo drug design for CO poisoning.
Collapse
Affiliation(s)
- Mingrui Jiang
- Laboratory
for Multiscale Material Modeling, Syracuse
University, 151L Link Hall, Syracuse, NY 13244, USA
- Department
of Civil and Environmental Engineering, Syracuse University, 151L Link Hall, Syracuse, NY 13244, USA
| | - Chi-Hua Yu
- Department
of Engineering Science, National Cheng Kung
University, No.1, University Road, Tainan City 701, Taiwan
| | - Zhiping Xu
- Applied
Mechanics Laboratory, Department of Engineering Mechanics, Tsinghua University, Beijing 100084, China
| | - Zhao Qin
- Laboratory
for Multiscale Material Modeling, Syracuse
University, 151L Link Hall, Syracuse, NY 13244, USA
- Department
of Civil and Environmental Engineering, Syracuse University, 151L Link Hall, Syracuse, NY 13244, USA
- The
BioInspired Institute, Syracuse University, Syracuse, NY 13244, USA
| |
Collapse
|
2
|
Sohraby F, Nunes-Alves A. Characterization of the Bottlenecks and Pathways for Inhibitor Dissociation from [NiFe] Hydrogenase. J Chem Inf Model 2024; 64:4193-4203. [PMID: 38728115 PMCID: PMC11134402 DOI: 10.1021/acs.jcim.4c00187] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2024] [Revised: 04/24/2024] [Accepted: 04/30/2024] [Indexed: 05/12/2024]
Abstract
[NiFe] hydrogenases can act as efficient catalysts for hydrogen oxidation and biofuel production. However, some [NiFe] hydrogenases are inhibited by gas molecules present in the environment, such as O2 and CO. One strategy to engineer [NiFe] hydrogenases and achieve O2- and CO-tolerant enzymes is by introducing point mutations to block the access of inhibitors to the catalytic site. In this work, we characterized the unbinding pathways of CO in the complex with the wild-type and 10 different mutants of [NiFe] hydrogenase from Desulfovibrio fructosovorans using τ-random accelerated molecular dynamics (τRAMD) to enhance the sampling of unbinding events. The ranking provided by the relative residence times computed with τRAMD is in agreement with experiments. Extensive data analysis of the simulations revealed that from the two bottlenecks proposed in previous studies for the transit of gas molecules (residues 74 and 122 and residues 74 and 476), only one of them (residues 74 and 122) effectively modulates diffusion and residence times for CO. We also computed pathway probabilities for the unbinding of CO, O2, and H2 from the wild-type [NiFe] hydrogenase, and we observed that while the most probable pathways are the same, the secondary pathways are different. We propose that introducing mutations to block the most probable paths, in combination with mutations to open the main secondary path used by H2, can be a feasible strategy to achieve CO and O2 resistance in the [NiFe] hydrogenase from Desulfovibrio fructosovorans.
Collapse
Affiliation(s)
- Farzin Sohraby
- Institute of Chemistry, Technische Universität Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Ariane Nunes-Alves
- Institute of Chemistry, Technische Universität Berlin, Straße des 17. Juni 135, 10623 Berlin, Germany
| |
Collapse
|
3
|
Zare M, Saleheen MS, Singh N, Uline MJ, Faheem M, Heyden A. Liquid-Phase Effects on Adsorption Processes in Heterogeneous Catalysis. JACS AU 2022; 2:2119-2134. [PMID: 36186571 PMCID: PMC9516566 DOI: 10.1021/jacsau.2c00389] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/08/2022] [Revised: 07/21/2022] [Accepted: 07/26/2022] [Indexed: 06/16/2023]
Abstract
Aqueous solvation free energies of adsorption have recently been measured for phenol adsorption on Pt(111). Endergonic solvent effects of ∼1 eV suggest solvents dramatically influence a metal catalyst's activity with significant implications for the catalyst design. However, measurements are indirect and involve adsorption isotherm models, which potentially reduces the reliability of the extracted energy values. Computational, implicit solvation models predict exergonic solvation effects for phenol adsorption, failing to agree with measurements even qualitatively. In this study, an explicit, hybrid quantum mechanical/molecular mechanical approach for computing solvation free energies of adsorption is developed, solvation free energies of phenol adsorption are computed, and experimental data for solvation free energies of phenol adsorption are reanalyzed using multiple adsorption isotherm models. Explicit solvation calculations predict an endergonic solvation free energy for phenol adsorption that agrees well with measurements to within the experimental and force field uncertainties. Computed adsorption free energies of solvation of carbon monoxide, ethylene glycol, benzene, and phenol over the (111) facet of Pt and Cu suggest that liquid water destabilizes all adsorbed species, with the largest impact on the largest adsorbates.
Collapse
Affiliation(s)
- Mehdi Zare
- Department
of Chemical Engineering, University of South
Carolina, 301 Main Street, Columbia, South Carolina 29208, United States
| | - Mohammad S. Saleheen
- Department
of Chemical Engineering, University of South
Carolina, 301 Main Street, Columbia, South Carolina 29208, United States
| | - Nirala Singh
- Department
of Chemical Engineering and Catalysis Science and Technology Institute, University of Michigan, Ann Arbor, Michigan 48109-2136, United States
| | - Mark J. Uline
- Department
of Chemical Engineering, University of South
Carolina, 301 Main Street, Columbia, South Carolina 29208, United States
| | - Muhammad Faheem
- Department
of Chemical Engineering, University of South
Carolina, 301 Main Street, Columbia, South Carolina 29208, United States
- Department
of Chemical Engineering, University of Engineering
& Technology, Lahore 54890, Pakistan
| | - Andreas Heyden
- Department
of Chemical Engineering, University of South
Carolina, 301 Main Street, Columbia, South Carolina 29208, United States
| |
Collapse
|
4
|
Amadei A, Aschi M. Stationary and Time-Dependent Carbon Monoxide Stretching Mode Features in Carboxy Myoglobin: A Theoretical-Computational Reappraisal. J Phys Chem B 2021; 125:13624-13634. [PMID: 34904432 DOI: 10.1021/acs.jpcb.1c05815] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The stationary and time-dependent infrared spectrum (IR) of the CO stretching mode (νCO) in carboxymyoglobin (MbCO), a longstanding problem of biophysical chemistry, has been modeled through a theoretical-computational method specifically designed for simulating quantum observables in complex atomic-molecular systems and based on a combined application of long time scale molecular dynamics simulations and quantum-chemical calculations. This study is basically focused on two aspects: (i) the origin of the stationary IR substates (termed as A0, A1, and A3) and (ii) the modeling and the interpretation of the νCO energy relaxation. The results, strengthened by a more than satisfactory agreement with the experimental data, concisely indicate that (i) the conformational His64-FeCO relevant substates, i.e., characterized by the formation-disruption of the H-bond between the above moieties, are the main responsible of the presence of two distinct and well separated (A0 and A1/A3) spectroscopic regions; (ii) the characteristic bimodal shape of the A1/A3 spectral region, according to our model, is the result of the fluctuation of the electric field pattern as provided by the protein-solvent framework perturbing the bound His64-CO-Heme complex; and (iii) the electric field pattern, in conjunction with the relatively high density of MbCO vibrational states, is also the main determinant of the νCO energy relaxation, characterizing its kinetic efficiency.
Collapse
Affiliation(s)
- Andrea Amadei
- Dipartimento di Scienze e Tecnologie Chimiche, Università di Roma "Tor Vergata", via della Ricerca Scientifica 1, 00 133 Roma, Italia
| | - Massimiliano Aschi
- Dipartimento di Scienze Fisiche e Chimiche, Università de l'Aquila, via Vetoio (Coppito 1), 67 010 l'Aquila, Italia
| |
Collapse
|
5
|
Bacellar C, Kinschel D, Cannelli O, Sorokin B, Katayama T, Mancini GF, Rouxel JR, Obara Y, Nishitani J, Ito H, Ito T, Kurahashi N, Higashimura C, Kudo S, Cirelli C, Knopp G, Nass K, Johnson PJM, Wach A, Szlachetko J, Lima FA, Milne CJ, Yabashi M, Suzuki T, Misawa K, Chergui M. Femtosecond X-ray spectroscopy of haem proteins. Faraday Discuss 2021; 228:312-328. [PMID: 33565544 DOI: 10.1039/d0fd00131g] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We discuss our recently reported femtosecond (fs) X-ray emission spectroscopy results on the ligand dissociation and recombination in nitrosylmyoglobin (MbNO) in the context of previous studies on ferrous haem proteins. We also present a preliminary account of femtosecond X-ray absorption studies on MbNO, pointing to the presence of more than one species formed upon photolysis.
Collapse
Affiliation(s)
- Camila Bacellar
- Laboratoire de Spectroscopie Ultrarapide (LSU), Lausanne Centre for Ultrafast Science (LACUS), École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland.
| | - Dominik Kinschel
- Laboratoire de Spectroscopie Ultrarapide (LSU), Lausanne Centre for Ultrafast Science (LACUS), École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland.
| | - Oliviero Cannelli
- Laboratoire de Spectroscopie Ultrarapide (LSU), Lausanne Centre for Ultrafast Science (LACUS), École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland.
| | - Boris Sorokin
- Laboratoire de Spectroscopie Ultrarapide (LSU), Lausanne Centre for Ultrafast Science (LACUS), École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland.
| | - Tetsuo Katayama
- Japan Synchrotron Radiation Research Institute (JASRI), 1-1-1, Kouto, Sayo-cho Sayo-gun, Hyogo 679-5198, Japan
| | - Giulia F Mancini
- Laboratoire de Spectroscopie Ultrarapide (LSU), Lausanne Centre for Ultrafast Science (LACUS), École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland.
| | - Jeremy R Rouxel
- Laboratoire de Spectroscopie Ultrarapide (LSU), Lausanne Centre for Ultrafast Science (LACUS), École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland.
| | - Yuki Obara
- Tokyo University of Agriculture and Technology (TUAT), 2-24-16 Naka-cho, Koganei, Tokyo 184-8588, Japan
| | - Junichi Nishitani
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa-Oiwakecho, Sakyo-Ku, Kyoto 606-8502, Japan
| | - Hironori Ito
- Tokyo University of Agriculture and Technology (TUAT), 2-24-16 Naka-cho, Koganei, Tokyo 184-8588, Japan
| | - Terumasa Ito
- Tokyo University of Agriculture and Technology (TUAT), 2-24-16 Naka-cho, Koganei, Tokyo 184-8588, Japan
| | - Naoya Kurahashi
- Department of Materials and Life Sciences, Faculty of Science and Technology, Sophia University, Kioicho, 7-1, Chiyoda, 102-8554 Tokyo, Japan
| | - Chika Higashimura
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa-Oiwakecho, Sakyo-Ku, Kyoto 606-8502, Japan
| | - Shotaro Kudo
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa-Oiwakecho, Sakyo-Ku, Kyoto 606-8502, Japan
| | - Claudio Cirelli
- SwissFEL, Paul Scherrer Institut (PSI), 5232 Villigen, Switzerland
| | - Gregor Knopp
- SwissFEL, Paul Scherrer Institut (PSI), 5232 Villigen, Switzerland
| | - Karol Nass
- SwissFEL, Paul Scherrer Institut (PSI), 5232 Villigen, Switzerland
| | | | - Anna Wach
- Institute of Nuclear Physics, Polish Academy of Sciences, 31-342 Krakow, Poland
| | - Jakub Szlachetko
- Institute of Nuclear Physics, Polish Academy of Sciences, 31-342 Krakow, Poland
| | | | | | - Makina Yabashi
- Japan Synchrotron Radiation Research Institute (JASRI), 1-1-1, Kouto, Sayo-cho Sayo-gun, Hyogo 679-5198, Japan
| | - Toshinori Suzuki
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa-Oiwakecho, Sakyo-Ku, Kyoto 606-8502, Japan
| | - Kazuhiko Misawa
- Tokyo University of Agriculture and Technology (TUAT), 2-24-16 Naka-cho, Koganei, Tokyo 184-8588, Japan
| | - Majed Chergui
- Laboratoire de Spectroscopie Ultrarapide (LSU), Lausanne Centre for Ultrafast Science (LACUS), École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland.
| |
Collapse
|
6
|
Wang X, Yan J, Zhang H, Xu Z, Zhang JZH. An electrostatic energy-based charge model for molecular dynamics simulation. J Chem Phys 2021; 154:134107. [PMID: 33832260 DOI: 10.1063/5.0043707] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Abstract
The interactions of the polar chemical bonds such as C=O and N-H with an external electric field were investigated, and a linear relationship between the QM/MM interaction energies and the electric field along the chemical bond is established in the range of weak to intermediate electrical fields. The linear relationship indicates that the electrostatic interactions of a polar group with its surroundings can be described by a simple model of a dipole with constant moment under the action of an electric field. This relationship is employed to develop a general approach to generating an electrostatic energy-based charge (EEC) model for molecules containing single or multiple polar chemical bonds. Benchmark test studies of this model were carried out for (CH3)2-CO and N-methyl acetamide in explicit water, and the result shows that the EEC model gives more accurate electrostatic energies than those given by the widely used charge model based on fitting to the electrostatic potential (ESP) in direct comparison to the energies computed by the QM/MM method. The MD simulations of the electric field at the active site of ketosteroid isomerase based on EEC demonstrated that EEC gave a better representation of the electrostatic interaction in the hydrogen-bonding environment than the Amber14SB force field by comparison with experiment. The current study suggests that EEC should be better suited for molecular dynamics study of molecular systems with polar chemical bonds such as biomolecules than the widely used ESP or RESP (restrained ESP) charge models.
Collapse
Affiliation(s)
- Xianwei Wang
- College of Science, Zhejiang University of Technology, Hangzhou, Zhejiang 310023, China
| | - Jinhua Yan
- College of Science, Zhejiang University of Technology, Hangzhou, Zhejiang 310023, China
| | - Hang Zhang
- College of Science, Zhejiang University of Technology, Hangzhou, Zhejiang 310023, China
| | - Zhousu Xu
- College of Science, Zhejiang University of Technology, Hangzhou, Zhejiang 310023, China
| | - John Z H Zhang
- Shanghai Engineering Research Center for Molecular Therapeutics and New Drug Development, Shanghai Key Laboratory of Green Chemistry and Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| |
Collapse
|
7
|
Patt A, Simon JM, Salazar JM, Picaud S. Adsorption of CO and N 2 molecules at the surface of solid water. A grand canonical Monte Carlo study. J Chem Phys 2020; 153:204502. [PMID: 33261471 DOI: 10.1063/5.0031254] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The adsorption of carbon monoxide and nitrogen molecules at the surface of four forms of solid water is investigated by means of grand canonical Monte Carlo simulations. The trapping ability of crystalline Ih and low-density amorphous ices, along with clathrate hydrates of structures I and II, is compared at temperatures relevant for astrophysics. It is shown that when considering a gas phase that contains mixtures of carbon monoxide and nitrogen, the trapping of carbon monoxide is favored with respect to nitrogen at the surface of all solids, irrespective of the temperature. The results of the calculations also indicate that some amounts of molecules can be incorporated in the bulk of the water structures, and the molecular selectivity of the incorporation process is investigated. Again, it is shown that incorporation of carbon monoxide is favored with respect to nitrogen in most of the situations considered here. In addition, the conclusions of the present simulations emphasize the importance of the strength of the interactions between the guest molecules and the water network. They indicate that the accuracy of the corresponding interaction potentials is a key point, especially for simulating clathrate selectivity. This highlights the necessity of having interaction potential models that are transferable to different water environments.
Collapse
Affiliation(s)
- Antoine Patt
- Institut UTINAM UMR 6213, CNRS/Université de Bourgogne Franche-Comté, Besançon, France
| | - Jean-Marc Simon
- Laboratoire Interdisciplinaire Carnot de Bourgogne (ICB) UMR 6303, CNRS, Université de Bourgogne Franche-Comté, F-21078 Dijon Cedex, France
| | - J Marcos Salazar
- Laboratoire Interdisciplinaire Carnot de Bourgogne (ICB) UMR 6303, CNRS, Université de Bourgogne Franche-Comté, F-21078 Dijon Cedex, France
| | - Sylvain Picaud
- Institut UTINAM UMR 6213, CNRS/Université de Bourgogne Franche-Comté, Besançon, France
| |
Collapse
|
8
|
Koner D, Salehi SM, Mondal P, Meuwly M. Non-conventional force fields for applications in spectroscopy and chemical
reaction dynamics. J Chem Phys 2020; 153:010901. [DOI: 10.1063/5.0009628] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Debasish Koner
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel,
Switzerland
| | - Seyedeh Maryam Salehi
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel,
Switzerland
| | - Padmabati Mondal
- Indian Institute of Science Education and Research (IISER) Tirupati, Karakambadi Road, Mangalam, Tirupati 517507, Andhra
Pradesh, India
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel,
Switzerland and Department of Chemistry, Brown University, Providence, Rhode Island, USA
| |
Collapse
|
9
|
Cho EH, Lin LC. Electrostatic Potential Optimized Molecular Models for Molecular Simulations: CO, CO 2, COS, H 2S, N 2, N 2O, and SO 2. J Chem Theory Comput 2019; 15:6323-6332. [PMID: 31618577 DOI: 10.1021/acs.jctc.9b00653] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular simulations have been widely employed in the discovery of nanoporous materials, such as metal-organic frameworks (MOFs) and zeolite, for energy- and environment-related applications. To achieve simulation predictions with better accuracy, we herein present a collection of molecular models, including carbon monoxide (CO), carbon dioxide (CO2), carbonyl sulfide (COS), hydrogen sulfide (H2S), nitrogen (N2), nitrous oxide (N2O), and sulfur dioxide (SO2). These models, denoted as electrostatic potential optimized molecular models (ESP-MMs), are systematically developed to not only reproduce experimental vapor-liquid equilibrium but also have accurate electrostatic potential representation surrounding the molecules. Our results show that, with accurate electrostatic potential representations, ESP-MMs can offer improved predictions in a variety of adsorption properties for porous materials, including MOFs with open-metal sites and all-silica zeolites. Specifically, by using ESP-MMs, the binding geometry and adsorption energy landscape can be well captured. This enables these models to be employed to unravel the fundamental mechanism of gaseous adsorption in materials of interest as well as to facilitate the parametrization of adsorbent-adsorbate interactions. We also demonstrate that, combined with generic force fields for adsorbents, ESP-MMs can offer reasonable predictions in adsorption isotherms. Although these ESP-MMs use a relatively simple and nonpolarizable potential form for the sake of efficiency and applicability, their accuracy has been extensively validated in this study. Furthermore, the set of Lennard-Jones potentials with static point charges adopted for ESP-MMs can be readily implemented in all available simulation packages. We anticipate that these ESP-MMs can largely facilitate future computational studies of porous materials for gas separation and removal.
Collapse
Affiliation(s)
- Eun Hyun Cho
- William G. Lowrie Department of Chemical and Biomolecular Engineering , The Ohio State University , Columbus , Ohio 43210 , United States
| | - Li-Chiang Lin
- William G. Lowrie Department of Chemical and Biomolecular Engineering , The Ohio State University , Columbus , Ohio 43210 , United States
| |
Collapse
|
10
|
Bringas M, Lombardi LE, Luque FJ, Estrin DA, Capece L. Ligand Binding Rate Constants in Heme Proteins Using Markov State Models and Molecular Dynamics Simulations. Chemphyschem 2019; 20:2451-2460. [PMID: 31365183 DOI: 10.1002/cphc.201900589] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2019] [Revised: 07/30/2019] [Indexed: 11/07/2022]
Abstract
Computer simulation studies of the molecular basis for ligand migration in proteins allow the description of key events such as the transition between docking sites, displacement of existing ligands and solvent molecules, and open/closure of specific "gates", among others. In heme proteins, ligand migration from the solvent to the active site preludes the binding to the heme iron and triggers different functions. In this work, molecular dynamics simulations, a Markov State Model of migration and empirical kinetic equations are combined to study the migration of O2 and NO in two truncated hemoglobins of Mycobacterium tuberculosis (Mt-TrHbN and Mt-TrHbO). For Mt-TrHbN, we show that the difference in the association constant in the oxy and deoxy states relies mainly in the displacement of water molecules anchored in the distal cavity in the deoxy form. The results here provide a valuable approach to study ligand migration in globins.
Collapse
Affiliation(s)
- Mauro Bringas
- Departamento de Química Inorgánica, Analítica y Química Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, C1428EGA, Buenos Aires, Argentina.,Instituto de Química Física de los Materiales, Medio Ambiente y Energía (INQUIMAE-CONICET), C1428EGA, Buenos Aires, Argentina
| | - Leandro E Lombardi
- Instituto de Cálculo, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires - CONICET, C1428EGA, Buenos Aires, Argentina
| | - F Javier Luque
- Department of Nutrition, Food Sciences and Gastronomy, Faculty of Pharmacy and Food Sciences, University of Barcelona, Campus Torribera, 08921, Santa Coloma de Gramenet, Spain.,Institute of Biomedicine (IBUB) and Institute of Theoretical and Computational Chemistry (IQTCUB), University of Barcelona, 08028, Barcelona, Spain
| | - Darío A Estrin
- Departamento de Química Inorgánica, Analítica y Química Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, C1428EGA, Buenos Aires, Argentina.,Instituto de Química Física de los Materiales, Medio Ambiente y Energía (INQUIMAE-CONICET), C1428EGA, Buenos Aires, Argentina
| | - Luciana Capece
- Departamento de Química Inorgánica, Analítica y Química Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, C1428EGA, Buenos Aires, Argentina.,Instituto de Química Física de los Materiales, Medio Ambiente y Energía (INQUIMAE-CONICET), C1428EGA, Buenos Aires, Argentina
| |
Collapse
|
11
|
Azar ANV, Velioglu S, Keskin S. Large-Scale Computational Screening of Metal Organic Framework (MOF) Membranes and MOF-Based Polymer Membranes for H 2/N 2 Separations. ACS SUSTAINABLE CHEMISTRY & ENGINEERING 2019; 7:9525-9536. [PMID: 31157127 PMCID: PMC6537470 DOI: 10.1021/acssuschemeng.9b01020] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/20/2019] [Revised: 04/04/2019] [Indexed: 05/05/2023]
Abstract
Several thousands of metal organic frameworks (MOFs) have been reported to date, but the information on H2/N2 separation performances of MOF membranes is currently very limited in the literature. We report the first large-scale computational screening study that combines state-of-the-art molecular simulations, grand canonical Monte Carlo (GCMC) and molecular dynamics (MD), to predict H2 permeability and H2/N2 selectivity of 3765 different types of MOF membranes. Results showed that MOF membranes offer very high H2 permeabilities, 2.5 × 103 to 1.7 × 106 Barrer, and moderate H2/N2 membrane selectivities up to 7. The top 20 MOF membranes that exceed the polymeric membranes' upper bound for H2/N2 separation were identified based on the results of initial screening performed at infinite dilution condition. Molecular simulations were then carried out considering binary H2/N2 and quaternary H2/N2/CO2/CO mixtures to evaluate the separation performance of MOF membranes under industrial operating conditions. Lower H2 permeabilities and higher N2 permeabilities were obtained at binary mixture conditions compared to the ones obtained at infinite dilution due to the absence of multicomponent mixture effects in the latter. Structure-performance relations of MOFs were also explored to provide molecular-level insights into the development of new MOF membranes that can offer both high H2 permeability and high H2/N2 selectivity. Results showed that the most promising MOF membranes generally have large pore sizes (>6 Å) as well as high surface areas (>3500 m2/g) and high pore volumes (>1 cm3/g). We finally examined H2/N2 separation potentials of the mixed matrix membranes (MMMs) in which the best MOF materials identified from our high-throughput screening were used as fillers in various polymers. Results showed that incorporation of MOFs into polymers almost doubles H2 permeabilities and slightly enhances H2/N2 selectivities of polymer membranes, which can advance the current membrane technology for efficient H2 purification.
Collapse
Affiliation(s)
- Ayda Nemati Vesali Azar
- Department of Chemical and Biological Engineering, Koc University, Rumelifeneri Yolu, Sariyer, 34450 Istanbul, Turkey
| | - Sadiye Velioglu
- Department of Chemical and Biological Engineering, Koc University, Rumelifeneri Yolu, Sariyer, 34450 Istanbul, Turkey
| | - Seda Keskin
- Department of Chemical and Biological Engineering, Koc University, Rumelifeneri Yolu, Sariyer, 34450 Istanbul, Turkey
| |
Collapse
|
12
|
Campbell C, Gomes JRB, Fischer M, Jorge M. New Model for Predicting Adsorption of Polar Molecules in Metal-Organic Frameworks with Unsaturated Metal Sites. J Phys Chem Lett 2018; 9:3544-3553. [PMID: 29886744 DOI: 10.1021/acs.jpclett.8b00967] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
Abstract
Conventional molecular models fail to correctly describe interactions of adsorbates with coordinatively unsaturated sites (CUS) present in a large number of metal-organic frameworks (MOFs). Here, we confirm the failure of these models for a prototypical polar adsorbate, carbon monoxide, and show that simply adjusting their parameters leads to poor agreement with experimental isotherms when outside the fitting conditions. We propose a new approach that combines quantum mechanical density functional theory (DFT) with Monte Carlo simulations to rigorously account for specific interactions at the CUS. By explicitly including electrostatic interactions and employing accurate DFT functionals that describe dispersion interactions, our modeling approach becomes generally applicable to both polar and nonpolar molecules. We demonstrate that this CUS model leads to substantial improvement in carbon monoxide adsorption isotherm predictions, and correctly captures the coordination binding mechanism. This paper represents a major stepping stone in the development of a robust, transferable and generally applicable approach to describe the complex interactions between gas molecules and CUS, with great potential for use in large-scale screening studies.
Collapse
Affiliation(s)
- Christopher Campbell
- Department of Chemical and Process Engineering , University of Strathclyde , 75 Montrose Street , Glasgow G1 1XJ , Scotland , United Kingdom
| | - José R B Gomes
- CICECO- Aveiro Institute of Materials, Department of Chemistry , University of Aveiro , Campus Universitário de Santiago , 3810-193 Aveiro , Portugal
| | - Michael Fischer
- Crystallography Group, Department of Geosciences , University of Bremen , Klagenfurter Straße , 28359 Bremen , Germany
- MAPEX Center for Materials and Processes , University of Bremen , 28359 Bremen , Germany
| | - Miguel Jorge
- Department of Chemical and Process Engineering , University of Strathclyde , 75 Montrose Street , Glasgow G1 1XJ , Scotland , United Kingdom
| |
Collapse
|
13
|
Chen G, Kong X, Lu D, Wu J, Liu Z. Kinetics of CO 2 diffusion in human carbonic anhydrase: a study using molecular dynamics simulations and the Markov-state model. Phys Chem Chem Phys 2018; 19:11690-11697. [PMID: 28435964 DOI: 10.1039/c7cp00887b] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Molecular dynamics (MD) simulations, in combination with the Markov-state model (MSM), were applied to probe CO2 diffusion from an aqueous solution into the active site of human carbonic anhydrase II (hCA-II), an enzyme useful for enhanced CO2 capture and utilization. The diffusion process in the hydrophobic pocket of hCA-II was illustrated in terms of a two-dimensional free-energy landscape. We found that CO2 diffusion in hCA-II is a rate-limiting step in the CO2 diffusion-binding-reaction process. The equilibrium distribution of CO2 shows its preferential accumulation within a hydrophobic domain in the protein core region. An analysis of the committors and reactive fluxes indicates that the main pathway for CO2 diffusion into the active site of hCA-II is through a binding pocket where residue Gln136 contributes to the maximal flux. The simulation results offer a new perspective on the CO2 hydration kinetics and useful insights toward the development of novel biochemical processes for more efficient CO2 sequestration and utilization.
Collapse
Affiliation(s)
- Gong Chen
- Key Laboratory of Industrial Biocatalysis, Ministry of Education, China.
| | | | | | | | | |
Collapse
|
14
|
Solvent Composition Drives the Rebinding Kinetics of Nitric Oxide to Microperoxidase. Sci Rep 2018; 8:5281. [PMID: 29588445 PMCID: PMC5869715 DOI: 10.1038/s41598-018-22944-z] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2017] [Accepted: 02/27/2018] [Indexed: 12/25/2022] Open
Abstract
The rebinding kinetics of NO after photodissociation from microperoxidase (Mp-9) is studied in different solvent environments. In mixed glycerol/water (G/W) mixtures the dissociating ligand rebinds with a yield close to 1 due to the cavities formed by the solvent whereas in pure water the ligand can diffuse into the solvent after photodissociation. In the G/W mixture, only geminate rebinding on the sub-picosecond and 5 ps time scales was found and the rebinding fraction is unity which compares well with available experiments. Contrary to that, simulations in pure water find two time scales – ~10 ps and ~200 ps - indicating that both, geminate rebinding and rebinding after diffusion of NO in the surrounding water contribute. The rebinding fraction is around 0.63 within 1 ns which is in stark contrast with experiment. Including ions (Na and Cl) at 0.15 M concentration in water leads to rebinding kinetics tending to that in the glycerol/water mixture and yields agreement with experiments. The effect of temperature is also probed and found to be non-negligible. The present simulations suggest that NO rebinding in Mp is primarily driven by thermal fluctuations which is consistent with recent resonance Raman spectroscopy experiments and simulations on MbNO.
Collapse
|
15
|
|
16
|
Wang XW, Zhang JZH, He X. Ab initio Quantum Mechanics/Molecular Mechanics Molecular Dynamics Simulation of CO in the Heme Distal Pocket of Myoglobin. CHINESE J CHEM PHYS 2017. [DOI: 10.1063/1674-0068/30/cjcp1709169] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Xian-wei Wang
- College of Science, Zhejiang University of Technology, Hangzhou 310023, China
- College of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- Zhejiang Provincial Collaborative Innovation Center of High-end Laser Manufacturing Equipment, Hangzhou 310014, China
| | - John Z. H. Zhang
- College of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Department of Chemistry, New York University, New York 10003, USA
| | - Xiao He
- College of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
17
|
Antipov SV, Bhattacharyya S, El Hage K, Xu ZH, Meuwly M, Rothlisberger U, Vaníček J. Ultrafast dynamics induced by the interaction of molecules with electromagnetic fields: Several quantum, semiclassical, and classical approaches. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2017; 4:061509. [PMID: 29376107 PMCID: PMC5758379 DOI: 10.1063/1.4996559] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/17/2017] [Accepted: 10/11/2017] [Indexed: 06/07/2023]
Abstract
Several strategies for simulating the ultrafast dynamics of molecules induced by interactions with electromagnetic fields are presented. After a brief overview of the theory of molecule-field interaction, we present several representative examples of quantum, semiclassical, and classical approaches to describe the ultrafast molecular dynamics, including the multiconfiguration time-dependent Hartree method, Bohmian dynamics, local control theory, semiclassical thawed Gaussian approximation, phase averaging, dephasing representation, molecular mechanics with proton transfer, and multipolar force fields. In addition to the general overview, some focus is given to the description of nuclear quantum effects and to the direct dynamics, in which the ab initio energies and forces acting on the nuclei are evaluated on the fly. Several practical applications, performed within the framework of the Swiss National Center of Competence in Research "Molecular Ultrafast Science and Technology," are presented: These include Bohmian dynamics description of the collision of H with H2, local control theory applied to the photoinduced ultrafast intramolecular proton transfer, semiclassical evaluation of vibrationally resolved electronic absorption, emission, photoelectron, and time-resolved stimulated emission spectra, infrared spectroscopy of H-bonding systems, and multipolar force fields applications in the condensed phase.
Collapse
Affiliation(s)
- Sergey V Antipov
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Swarnendu Bhattacharyya
- Laboratory of Computational Chemistry and Biochemistry, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Krystel El Hage
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Zhen-Hao Xu
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
18
|
Waga M, Takeda S, Sakata R. A method for extracting carboxy-myoglobin from beef. Anim Sci J 2017; 89:467-473. [PMID: 29052941 DOI: 10.1111/asj.12942] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2015] [Accepted: 08/17/2017] [Indexed: 11/30/2022]
Abstract
The present study was conducted to devise a method for the effective extraction of carboxy-myoglobin (COMb) from beef without carbon monoxide dissociation. The ratio of COMb to myoglobin was computed at absorptions of wavelengths 541 and 551 nm, which characterize COMb and the isosbestic point between COMb and deoxy-myoglobin, respectively. The COMb extraction rate was found to vary with temperature, pH and oxygen conditions. The decrease observed in this rate was inversely proportional to the rise in extraction temperature. The COMb extraction rate was also affected by pH, and the stability of COMb in the extract solution was the highest at pH 8.0-9.0. Moreover, the presence of oxygen was found to disturb COMb extraction. According to these results, nearly all COMb could be extracted from carbon-monoxide-treated beef under stirring conditions in pH 8.5 deoxidized buffer, at 1°C, and under N2 flow with the improved extraction method in this study (98.1 ± 2.7%). The decrement of COMb in the extract was accelerated by light, and the COMb was stable for 20 min in the dark, at 1°C. The extraction conditions for COMb described above should allow the accurate evaluation of COMb in meat tissue.
Collapse
Affiliation(s)
- Masahiro Waga
- Graduate School of Veterinary Science, Azabu University, Sagamihara, Japan.,Central Research Institute, Itoham Foods Inc., Moriya, Japan
| | - Shiro Takeda
- Department of Animal Science and Biotechnology, School of Veterinary Medicine, Azabu University, Sagamihara, Japan
| | - Ryoichi Sakata
- Graduate School of Veterinary Science, Azabu University, Sagamihara, Japan.,Department of Animal Science and Biotechnology, School of Veterinary Medicine, Azabu University, Sagamihara, Japan
| |
Collapse
|
19
|
Awoonor-Williams E, Rowley CN. The hydration structure of carbon monoxide byab initiomethods. J Chem Phys 2017; 146:034503. [DOI: 10.1063/1.4974164] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
|
20
|
Zhao W, Francisco JS, Zeng XC. CO Separation from H 2 via Hydrate Formation in Single-Walled Carbon Nanotubes. J Phys Chem Lett 2016; 7:4911-4915. [PMID: 27934039 DOI: 10.1021/acs.jpclett.6b02443] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Hydrogen is an alternative fuel without generating greenhouse gas or other harmful emissions. Industrial hydrogen production, however, always contains a small fraction of carbon monoxide (CO) (∼0.5-2%) that must be removed for use in fuel cells. Here, we present molecular dynamics simulation evidence on facile separation of CO from H2 at ambient pressure via the formation of quasi-one-dimensional (Q1D) clathrate hydrates within single-walled carbon nanotubes (SW-CNTs). At ambient pressure, Q1D CO (or H2) clathrates in SW-CNTs are formed spontaneously when the SW-CNTs are immersed in CO (or H2) aqueous solution. More interestingly, for the CO/H2 aqueous solution, highly preferential adsorption of CO over H2 occurs within the octagonal or nonagonal ice nanotubes inside of SW-CNTs. These results suggest that the formation of Q1D hydrates within SW-CNTs can be a viable and safe method for the separation of CO from H2, which can be exploited for hydrogen purification in fuel cells.
Collapse
Affiliation(s)
- Wenhui Zhao
- Department of Physics, Ningbo University , Ningbo, Zhejiang 315211, China
- Collaborative Innovation Center of Chemistry for Energy Materials, Department of Chemical Physics, University of Science and Technology of China , Hefei, Anhui 230026, China
| | - Joseph S Francisco
- Department of Chemistry, University of Nebraska-Lincoln , Lincoln, Nebraska 68588, United States
| | - Xiao Cheng Zeng
- Collaborative Innovation Center of Chemistry for Energy Materials, Department of Chemical Physics, University of Science and Technology of China , Hefei, Anhui 230026, China
- Department of Chemistry, University of Nebraska-Lincoln , Lincoln, Nebraska 68588, United States
| |
Collapse
|
21
|
Andrew CR, Petrova ON, Lamarre I, Lambry JC, Rappaport F, Negrerie M. The Dynamics Behind the Affinity: Controlling Heme-Gas Affinity via Geminate Recombination and Heme Propionate Conformation in the NO Carrier Cytochrome c'. ACS Chem Biol 2016; 11:3191-3201. [PMID: 27709886 DOI: 10.1021/acschembio.6b00599] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Nitric oxide (NO) sensors are heme proteins which may also bind CO and O2. Control of heme-gas affinity and their discrimination are achieved by the structural properties and reactivity of the heme and its distal and proximal environments, leading to several energy barriers. In the bacterial NO sensor cytochrome c' from Alcaligenes xylosoxidans (AXCP), the single Leu16Ala distal mutation boosts the affinity for gas ligands by a remarkable 106-108-fold, transforming AXCP from one of the lowest affinity gas binding proteins to one of the highest. Here, we report the dynamics of diatomics after photodissociation from wild type and L16A-AXCP over 12 orders of magnitude in time. For the L16A variant, the picosecond geminate rebinding of both CO and NO appears with an unprecedented 100% yield, and no exit of these ligands from protein to solvent could be observed. Molecular dynamic simulations saliently demonstrate that dissociated CO stays within 4 Å from Fe2+, in contrast to wild-type AXCP. The L16A mutation confers a heme propionate conformation and docking site which traps the diatomics, maximizing the probability of recombination and directly explaining the ultrahigh affinities for CO, NO, and O2. Overall, our results point to a novel mechanism for modulating heme-gas affinities in proteins.
Collapse
Affiliation(s)
- Colin R. Andrew
- Department
of Chemistry and Biochemistry, Eastern Oregon University, La Grande, Oregon 97850, United States
| | - Olga N. Petrova
- Laboratoire
d’Optique et Biosciences, INSERM, Ecole Polytechnique, 91128 Palaiseau, France
| | - Isabelle Lamarre
- Laboratoire
d’Optique et Biosciences, INSERM, Ecole Polytechnique, 91128 Palaiseau, France
| | - Jean-Christophe Lambry
- Laboratoire
d’Optique et Biosciences, INSERM, Ecole Polytechnique, 91128 Palaiseau, France
| | - Fabrice Rappaport
- Laboratoire
de Physiologie Membranaire et Moléculaire du Chloroplaste, CNRS, Université Pierre et Marie Curie, 75005 Paris, France
| | - Michel Negrerie
- Laboratoire
d’Optique et Biosciences, INSERM, Ecole Polytechnique, 91128 Palaiseau, France
| |
Collapse
|
22
|
Xiong Y, Fan YZ, Damasceno Borges D, Chen CX, Wei ZW, Wang HP, Pan M, Jiang JJ, Maurin G, Su CY. Ligand and Metal Effects on the Stability and Adsorption Properties of an Isoreticular Series of MOFs Based on T-Shaped Ligands and Paddle-Wheel Secondary Building Units. Chemistry 2016; 22:16147-16156. [DOI: 10.1002/chem.201603299] [Citation(s) in RCA: 42] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2016] [Indexed: 12/21/2022]
Affiliation(s)
- Ying Xiong
- MOE Laboratory of Bioinorganic and Synthetic Chemistry; Lehn Institute of Functional Materials; Sun Yat-Sen University; Guangzhou 510275 P.R. China
| | - Yan-Zhong Fan
- MOE Laboratory of Bioinorganic and Synthetic Chemistry; Lehn Institute of Functional Materials; Sun Yat-Sen University; Guangzhou 510275 P.R. China
| | - Daiane Damasceno Borges
- Institut Charles Gerhardt Montpellier; UMR 5253 CNRS, UM, ENSCM; Université Montpellier; Pl. E. Bataillon 34095 Montpellier Cedex 05 France
- Institute of Physics “Gleb Wataghin“; University of Campinas; Campinas - SP 13083-970 Brazil
| | - Cheng-Xia Chen
- MOE Laboratory of Bioinorganic and Synthetic Chemistry; Lehn Institute of Functional Materials; Sun Yat-Sen University; Guangzhou 510275 P.R. China
| | - Zhang-Wen Wei
- MOE Laboratory of Bioinorganic and Synthetic Chemistry; Lehn Institute of Functional Materials; Sun Yat-Sen University; Guangzhou 510275 P.R. China
| | - Hai-Ping Wang
- Institut Charles Gerhardt Montpellier; UMR 5253 CNRS, UM, ENSCM; Université Montpellier; Pl. E. Bataillon 34095 Montpellier Cedex 05 France
| | - Mei Pan
- MOE Laboratory of Bioinorganic and Synthetic Chemistry; Lehn Institute of Functional Materials; Sun Yat-Sen University; Guangzhou 510275 P.R. China
| | - Ji-Jun Jiang
- MOE Laboratory of Bioinorganic and Synthetic Chemistry; Lehn Institute of Functional Materials; Sun Yat-Sen University; Guangzhou 510275 P.R. China
| | - Guillaume Maurin
- Institut Charles Gerhardt Montpellier; UMR 5253 CNRS, UM, ENSCM; Université Montpellier; Pl. E. Bataillon 34095 Montpellier Cedex 05 France
| | - Cheng-Yong Su
- MOE Laboratory of Bioinorganic and Synthetic Chemistry; Lehn Institute of Functional Materials; Sun Yat-Sen University; Guangzhou 510275 P.R. China
- State Key Laboratory of Applied Organic Chemistry; Lanzhou University; Lanzhou 730000 P.R. China
| |
Collapse
|
23
|
Ultrafast anisotropic protein quake propagation after CO photodissociation in myoglobin. Proc Natl Acad Sci U S A 2016; 113:10565-70. [PMID: 27601659 DOI: 10.1073/pnas.1603539113] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
"Protein quake" denotes the dissipation of excess energy across a protein, in response to a local perturbation such as the breaking of a chemical bond or the absorption of a photon. Femtosecond time-resolved small- and wide-angle X-ray scattering (TR-SWAXS) is capable of tracking such ultrafast protein dynamics. However, because the structural interpretation of the experiments is complicated, a molecular picture of protein quakes has remained elusive. In addition, new questions arose from recent TR-SWAXS data that were interpreted as underdamped oscillations of an entire protein, thus challenging the long-standing concept of overdamped global protein dynamics. Based on molecular-dynamics simulations, we present a detailed molecular movie of the protein quake after carbon monoxide (CO) photodissociation in myoglobin. The simulations suggest that the protein quake is characterized by a single pressure peak that propagates anisotropically within 500 fs across the protein and further into the solvent. By computing TR-SWAXS patterns from the simulations, we could interpret features in the reciprocal-space SWAXS signals as specific real-space dynamics, such as CO displacement and pressure wave propagation. Remarkably, we found that the small-angle data primarily detect modulations of the solvent density but not oscillations of the bare protein, thereby reconciling recent TR-SWAXS experiments with the notion of overdamped global protein dynamics.
Collapse
|
24
|
Schubert A, Falvo C, Meier C. Mixed quantum-classical simulations of the vibrational relaxation of photolyzed carbon monoxide in a hemoprotein. J Chem Phys 2016; 145:054108. [DOI: 10.1063/1.4959859] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alexander Schubert
- Laboratoire Collisions Agrégats et Réactivité, IRSAMC, UMR CNRS 5589, Université Paul Sabatier, 31062 Toulouse, France
| | - Cyril Falvo
- Institut des Sciences Moléculaires d’Orsay (ISMO), CNRS, Univ. Paris-Sud, Université Paris-Saclay, 91405 Orsay, France
| | - Christoph Meier
- Laboratoire Collisions Agrégats et Réactivité, IRSAMC, UMR CNRS 5589, Université Paul Sabatier, 31062 Toulouse, France
| |
Collapse
|
25
|
El Hage K, Bereau T, Jakobsen S, Meuwly M. Impact of Quadrupolar Electrostatics on Atoms Adjacent to the Sigma-Hole in Condensed-Phase Simulations. J Chem Theory Comput 2016; 12:3008-19. [PMID: 27158892 DOI: 10.1021/acs.jctc.6b00202] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Halogenation is one of the cases for which advanced molecular simulation methods are mandatory for quantitative and predictive studies. The present work provides a systematic investigation of the importance of higher-order multipoles on specific sites of halobenzenes, other than the halogen, for static and dynamic properties in condensed-phase simulations. For that purpose, solute-solvent interactions using point charge (PC), multipole (MTP), and hybrid point charge/multipole (HYB) electrostatic models are analyzed in regions of halogen bonding and extended to regions of π orbitals of phenyl carbons. Using molecular dynamics simulations and quantum chemical methods, it is found that the sigma-hole does not only affect the halogen and the carbon bound to it but its effect extends to the carbons adjacent to the CX bond. This effect increases with the magnitude of the positive potential of the sigma-hole. With the MTP and HYB3 models, all hydration free energies of the PhX compounds are reproduced within 0.1 kcal/mol. Analysis of pair distribution functions and hydration free energies of halogenated benzenes provides a microscopic explanation why "point charge"-based representations with off-site charges fail in reproducing thermodynamic properties of the sigma-hole. Application of the hybrid models to study protein-ligand binding demonstrates both their accuracy and computational efficiency.
Collapse
Affiliation(s)
- Krystel El Hage
- Department of Chemistry, University of Basel , Klingelbergstrasse 80, 4056 Basel, Switzerland
| | - Tristan Bereau
- Max Planck Institute for Polymer Research , Ackermannweg 10, 55128 Mainz, Germany
| | - Sofie Jakobsen
- Department of Chemistry, University of Basel , Klingelbergstrasse 80, 4056 Basel, Switzerland.,Department of Chemistry, Aarhus University , Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Markus Meuwly
- Department of Chemistry, University of Basel , Klingelbergstrasse 80, 4056 Basel, Switzerland
| |
Collapse
|
26
|
Nicolaides A, Soulimane T, Varotsis C. Nanosecond ligand migration and functional protein relaxation in ba3 oxidoreductase: Structures of the B0, B1 and B2 intermediate states. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2016; 1857:1534-1540. [PMID: 27207588 DOI: 10.1016/j.bbabio.2016.05.002] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/11/2016] [Revised: 04/13/2016] [Accepted: 05/15/2016] [Indexed: 11/28/2022]
Abstract
Nanosecond time-resolved step-scan FTIR spectroscopy (nTRS (2) -FTIR) has been applied to literally probe the active site of the carbon monoxide (CO)-bound thermophilic ba3 heme-copper oxidoreductase as it executes its function. The nTRS (2) - snapshots of the photolysed heme a3 Fe-CO/CuB species captured a "transition state" whose side chains prevent the photolysed CO to enter the docking cavity. There are three sets of ba3 photoproduct bands of docked CO with different orientation exhibiting different kinetics. The trajectories of the "docked" CO at 2122, 2129 and 2137cm(-1) is referred to in the literature as B2, B1 and B0 intermediate states, respectively. The present data provided direct evidence for the role of water in controlling ligand orientation in an intracavity protein environment.
Collapse
Affiliation(s)
- Antonis Nicolaides
- Cyprus University of Technology, Department of Environmental Science and Technology, P.O. Box 50329, 3603 Lemesos, Cyprus
| | - Tewfik Soulimane
- Chemical and Environmental Science Department, University of Limerick, Limerick, Ireland; Materials & Surface Science Institute, University of Limerick, Limerick, Ireland
| | - Constantinos Varotsis
- Cyprus University of Technology, Department of Environmental Science and Technology, P.O. Box 50329, 3603 Lemesos, Cyprus.
| |
Collapse
|
27
|
Danielsson J, Meuwly M. Atomistic Simulation of Adiabatic Reactive Processes Based on Multi-State Potential Energy Surfaces. J Chem Theory Comput 2015; 4:1083-93. [PMID: 26636362 DOI: 10.1021/ct800066q] [Citation(s) in RCA: 64] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The adiabatic reactive molecular dynamics (ARMD) method provides a framework to study chemical reactions using molecular dynamics simulations with minimal computational overhead. Here, ARMD is generalized to an arbitrary reactive process between two states in which reactants and products can be treated by an atomistic force field. The implementation is described, and the method is applied to two systems: the kinetics of NO rebinding to myoglobin (Mb) as a validation system and the conformational transition in neuroglobin (Ngb) which explores the full functionality of ARMD. For MbNO, the nonexponential kinetics observed both in experiment and earlier ARMD studies is reproduced. Furthermore, the sensitivity of the results with respect to the asymptotic separation between the two potential energy surfaces (NO bound and unbound) is studied.
Collapse
Affiliation(s)
- Jonas Danielsson
- Chemistry Department, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Markus Meuwly
- Chemistry Department, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| |
Collapse
|
28
|
De Sancho D, Kubas A, Wang PH, Blumberger J, Best RB. Identification of Mutational Hot Spots for Substrate Diffusion: Application to Myoglobin. J Chem Theory Comput 2015; 11:1919-27. [PMID: 26574395 PMCID: PMC6132223 DOI: 10.1021/ct5011455] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The pathways by which small molecules (substrates or inhibitors) access active sites are a key aspect of the function of enzymes and other proteins. A key problem in designing or altering such proteins is to identify sites for mutation that will have the desired effect on the substrate transport properties. While specific access channels have been invoked in the past, molecular simulations suggest that multiple routes are possible, complicating the analysis. This complexity, however, can be captured by a Markov State Model (MSM) of the ligand diffusion process. We have developed a sensitivity analysis of the resulting rate matrix, which identifies the locations where mutations should have the largest effect on the diffusive on rate. We apply this method to myoglobin, which is the best characterized example both from experiment and simulation. We validate the approach by translating the sensitivity parameter obtained from this method into the CO binding rates in myoglobin upon mutation, resulting in a semi-quantitative correlation with experiments. The model is further validated against an explicit simulation for one of the experimental mutants.
Collapse
Affiliation(s)
- David De Sancho
- Department of Chemistry, University of Cambridge , Lensfield Road, Cambridge CB2 1EW, United Kingdom
- CIC nanoGUNE , Tolosa Hiribidea 76, 20018 Donostia-San Sebastian, Spain
- IKERBASQUE, Basque Foundation for Science , Maria Diaz de Haro 3, 48013 Bilbao, Spain
| | - Adam Kubas
- Department of Physics and Astronomy, University College London , Gower Street, London WC1E 6BT, United Kingdom
| | - Po-Hung Wang
- Department of Physics and Astronomy, University College London , Gower Street, London WC1E 6BT, United Kingdom
- Theoretical Molecular Science Laboratory , 2-1 Hirosawa, Wako-shi, Saitama 351-0198, Japan
| | - Jochen Blumberger
- Department of Physics and Astronomy, University College London , Gower Street, London WC1E 6BT, United Kingdom
| | - Robert B Best
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health , Bethesda, Maryland 20892-0520, United States
| |
Collapse
|
29
|
Ghesquière P, Mineva T, Talbi D, Theulé P, Noble JA, Chiavassa T. Diffusion of molecules in the bulk of a low density amorphous ice from molecular dynamics simulations. Phys Chem Chem Phys 2015; 17:11455-68. [DOI: 10.1039/c5cp00558b] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Arguments for a solvent driven mechanism for the diffusion of CO, CO2, NH3, and H2CO in a LDA water ice.
Collapse
Affiliation(s)
- P. Ghesquière
- Laboratoire Univers et Particules de Montpellier UMR 5299
- CNRS et Université de Montpellier
- Montpellier Cedex 05
- France
| | - T. Mineva
- Institut Charles Gerhardt
- UMR 5253 CNRS/ENSCM/UM2/UM1
- 34296 Montpellier Cedex 05
- France
| | - D. Talbi
- Laboratoire Univers et Particules de Montpellier UMR 5299
- CNRS et Université de Montpellier
- Montpellier Cedex 05
- France
| | - P. Theulé
- Aix-Marseille Univ
- CNRS
- Marseille
- France
| | | | | |
Collapse
|
30
|
Takayanagi M, Kurisaki I, Nagaoka M. Non-site-specific allosteric effect of oxygen on human hemoglobin under high oxygen partial pressure. Sci Rep 2014; 4:4601. [PMID: 24710521 PMCID: PMC3978498 DOI: 10.1038/srep04601] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2013] [Accepted: 03/12/2014] [Indexed: 12/11/2022] Open
Abstract
Protein allostery is essential for vital activities. Allosteric regulation of human hemoglobin (HbA) with two quaternary states T and R has been a paradigm of allosteric structural regulation of proteins. It is widely accepted that oxygen molecules (O2) act as a “site-specific” homotropic effector, or the successive O2 binding to the heme brings about the quaternary regulation. However, here we show that the site-specific allosteric effect is not necessarily only a unique mechanism of O2 allostery. Our simulation results revealed that the solution environment of high O2 partial pressure enhances the quaternary change from T to R without binding to the heme, suggesting an additional “non-site-specific” allosteric effect of O2. The latter effect should play a complementary role in the quaternary change by affecting the intersubunit contacts. This analysis must become a milestone in comprehensive understanding of the allosteric regulation of HbA from the molecular point of view.
Collapse
Affiliation(s)
- Masayoshi Takayanagi
- 1] Venture Business Laboratory, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan [2] Graduate School of Information Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan [3] Core Research for Evolutional Science and Technology, Japan Science and Technology Agency, Honmachi, Kawaguchi 332-0012, Japan
| | - Ikuo Kurisaki
- Graduate School of Information Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan
| | - Masataka Nagaoka
- 1] Graduate School of Information Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan [2] Core Research for Evolutional Science and Technology, Japan Science and Technology Agency, Honmachi, Kawaguchi 332-0012, Japan
| |
Collapse
|
31
|
Lee MW, Carr JK, Göllner M, Hamm P, Meuwly M. 2D IR spectra of cyanide in water investigated by molecular dynamics simulations. J Chem Phys 2014; 139:054506. [PMID: 23927269 DOI: 10.1063/1.4815969] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
Using classical molecular dynamics simulations, the 2D infrared (IR) spectroscopy of CN(-) solvated in D2O is investigated. Depending on the force field parametrizations, most of which are based on multipolar interactions for the CN(-) molecule, the frequency-frequency correlation function and observables computed from it differ. Most notably, models based on multipoles for CN(-) and TIP3P for water yield quantitatively correct results when compared with experiments. Furthermore, the recent finding that T1 times are sensitive to the van der Waals ranges on the CN(-) is confirmed in the present study. For the linear IR spectrum, the best model reproduces the full widths at half maximum almost quantitatively (13.0 cm(-1) vs. 14.9 cm(-1)) if the rotational contribution to the linewidth is included. Without the rotational contribution, the lines are too narrow by about a factor of two, which agrees with Raman and IR experiments. The computed and experimental tilt angles (or nodal slopes) α as a function of the 2D IR waiting time compare favorably with the measured ones and the frequency fluctuation correlation function is invariably found to contain three time scales: a sub-ps, 1 ps, and one on the 10-ps time scale. These time scales are discussed in terms of the structural dynamics of the surrounding solvent and it is found that the longest time scale (≈10 ps) most likely corresponds to solvent exchange between the first and second solvation shell, in agreement with interpretations from nuclear magnetic resonance measurements.
Collapse
Affiliation(s)
- Myung Won Lee
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | | | | | | | | |
Collapse
|
32
|
Zhang P, Ahn SW, Straub JE. “Strange Kinetics” in the Temperature Dependence of Methionine Ligand Rebinding Dynamics in Cytochrome c. J Phys Chem B 2013; 117:7190-202. [DOI: 10.1021/jp400481m] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Ping Zhang
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| | - Steven Wooseok Ahn
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| | - John E. Straub
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| |
Collapse
|
33
|
Cazade PA, Meuwly M. Oxygen migration pathways in NO-bound truncated hemoglobin. Chemphyschem 2012; 13:4276-86. [PMID: 23161831 DOI: 10.1002/cphc.201200608] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2012] [Revised: 09/28/2012] [Indexed: 11/10/2022]
Abstract
Atomistic simulations of dioxygen (O(2)) dynamics and migration in nitric oxide-bound truncated Hemoglobin N (trHbN) of Mycobacterium tuberculosis are reported. From more than 100 ns of simulations the connectivity network involving the metastable states for localization of the O(2) ligand is built and analyzed. It is found that channel I is the primary entrance point for O(2) whereas channel II is predominantly an exit path although access to the protein active site is also possible. For O(2) a new site compared to nitric oxide, from which reaction with the heme group can occur, was found. As this site is close to the heme iron, it could play an important role in the dioxygenation mechanism as O(2) can remain there for hundreds of picoseconds after which it can eventually leave the protein, while NO is localized in Xe2. The present study supports recent experimental work which proposed that O(2) docks in alternative pockets than Xe close to the reactive site. Similar to other proteins, a phenylalanine residue (Phe62) plays the role of a gate along the access route in channel I. The most highly connected site is the Xe3 pocket which is a "hub" and free energy barriers between the different metastable states are ≈1.5 kcal mol(-1) which allows facile O(2) migration within the protein.
Collapse
Affiliation(s)
- Pierre-André Cazade
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056, Basel, Switzerland
| | | |
Collapse
|
34
|
Pérez-Yáñez S, Beobide G, Castillo O, Fischer M, Hoffmann F, Fröba M, Cepeda J, Luque A. Gas Adsorption Properties and Selectivity in CuII/Adeninato/Carboxylato Metal-Biomolecule Frameworks. Eur J Inorg Chem 2012. [DOI: 10.1002/ejic.201201028] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
35
|
Tsuduki T, Tomita A, Koshihara SY, Adachi SI, Yamato T. Ligand migration in myoglobin: a combined study of computer simulation and x-ray crystallography. J Chem Phys 2012; 136:165101. [PMID: 22559505 DOI: 10.1063/1.4704586] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
A ligand-migration mechanism of myoglobin was studied by a multidisciplinary approach that used x-ray crystallography and molecular dynamics simulation. The former revealed the structural changes of the protein along with the ligand migration, and the latter provided the statistical ensemble of protein conformations around the thermal average. We developed a novel computational method, homogeneous ensemble displacement, and generated the conformational ensemble of ligand-detached species from that of ligand-bound species. The thermally averaged ligand-protein interaction was illustrated in terms of the potential of mean force. Although the structural changes were small, the presence of the ligand molecule in the protein matrix significantly affected the 3D scalar field of the potential of mean force, in accordance with the self-opening model proposed in the previous x-ray study.
Collapse
Affiliation(s)
- Takayuki Tsuduki
- Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan
| | | | | | | | | |
Collapse
|
36
|
Plattner N, Meuwly M. Quantifying the importance of protein conformation on ligand migration in myoglobin. Biophys J 2012; 102:333-41. [PMID: 22339870 DOI: 10.1016/j.bpj.2011.10.058] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2011] [Revised: 10/24/2011] [Accepted: 10/26/2011] [Indexed: 10/14/2022] Open
Abstract
Myoglobin (Mb) is a model system for ligand binding and migration. The energy barriers (ΔG) for ligand migration in Mb have been studied in the past by experiment and theory and significant differences between different approaches were found. From experiment, it is known that Mb can assume a large number of conformational substates. In this work, these substates are investigated as a possible source of the differences in migration barriers. We show that the initial structure significantly affects the calculated ΔG for a particular transition and that fluctuations in barrier heights δΔG are of similar magnitude as the free energy barriers themselves. The sensitivity of ΔG to the initial structure is compared to other sources of errors. Different protein structures can affect the calculated ΔG by up to 4 kcal/mol, whereas differences between simple point charge models and more elaborate multipolar charge models for the CO-ligand are smaller by a factor of two. Analysis of the structural changes underlying the large effect of the conformational substate reveals the importance of coupling between protein and ligand motion for migration.
Collapse
Affiliation(s)
- Nuria Plattner
- Chemistry Department, Brown University, Providence, Rhode Island, USA
| | | |
Collapse
|
37
|
Kim S, Park J, Lee T, Lim M. Direct Observation of Ligand Rebinding Pathways in Hemoglobin Using Femtosecond Mid-IR Spectroscopy. J Phys Chem B 2012; 116:6346-55. [PMID: 22587393 DOI: 10.1021/jp3026495] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Seongheun Kim
- Department of Chemistry and Chemistry
Institute for
Functional Materials, Pusan National University, Busan 609-735, Korea
| | - Jaeheung Park
- Department of Chemistry and Chemistry
Institute for
Functional Materials, Pusan National University, Busan 609-735, Korea
| | - Taegon Lee
- Department of Chemistry and Chemistry
Institute for
Functional Materials, Pusan National University, Busan 609-735, Korea
| | - Manho Lim
- Department of Chemistry and Chemistry
Institute for
Functional Materials, Pusan National University, Busan 609-735, Korea
| |
Collapse
|
38
|
Zhang P, Małolepsza E, Straub JE. Dynamics of Methionine Ligand Rebinding in Cytochrome c. J Phys Chem B 2012; 116:6980-90. [DOI: 10.1021/jp300783j] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ping Zhang
- Department
of Chemistry, Boston University, Boston,
Massachusetts 02215, United States
| | - Edyta Małolepsza
- Department
of Chemistry, Boston University, Boston,
Massachusetts 02215, United States
| | - John E. Straub
- Department
of Chemistry, Boston University, Boston,
Massachusetts 02215, United States
| |
Collapse
|
39
|
Fischer M, Hoffmann F, Fröba M. Metal–organic frameworks and related materials for hydrogen purification: Interplay of pore size and pore wall polarity. RSC Adv 2012. [DOI: 10.1039/c2ra01239a] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
|
40
|
Huang H, Zhang W, Liu D, Liu B, Chen G, Zhong C. Effect of temperature on gas adsorption and separation in ZIF-8: A combined experimental and molecular simulation study. Chem Eng Sci 2011. [DOI: 10.1016/j.ces.2011.09.009] [Citation(s) in RCA: 131] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
41
|
Csányi É, Ható Z, Kristóf T. Molecular simulation of water removal from simple gases with zeolite NaA. J Mol Model 2011; 18:2349-56. [DOI: 10.1007/s00894-011-1253-7] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2011] [Accepted: 09/21/2011] [Indexed: 11/24/2022]
|
42
|
|
43
|
Anselmi M, Di Nola A, Amadei A. Kinetics of carbon monoxide migration and binding in solvated neuroglobin as revealed by molecular dynamics simulations and quantum mechanical calculations. J Phys Chem B 2011; 115:2436-46. [PMID: 21332165 DOI: 10.1021/jp110833v] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Neuroglobin (Ngb) is a globular protein that reversibly binds small ligands at the six coordination position of the heme. With respect to other globins similar to myoglobin, Ngb displays some peculiarities as the topological reorganization of the internal cavities coupled to the sliding of the heme, or the binding of the endogenous distal histidine to the heme in the absence of an exogenous ligand. In this Article, by using multiple (independent) molecular dynamics trajectories (about 500 ns in total), the migration pathways of photolized carbon monoxide (CO) within solvated Ngb were analyzed, and a quantitative description of CO migration and corresponding kinetics was obtained. MD results, combined with quantum mechanical calculations on the CO-heme binding-unbinding reaction step in Ngb, allowed construction of a quantitative model representing the relevant steps of CO migration and rebinding.
Collapse
|
44
|
Keyes T, Napoleon RL. Extending Classical Molecular Theory with Polarization. J Phys Chem B 2010; 115:522-31. [DOI: 10.1021/jp105595q] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Tom Keyes
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| | - Raeanne L. Napoleon
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| |
Collapse
|
45
|
Anselmi M, Di Nola A, Amadei A. The effects of the L29F mutation on the ligand migration kinetics in crystallized myoglobin as revealed by molecular dynamics simulations. Proteins 2010; 79:867-79. [DOI: 10.1002/prot.22924] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2010] [Accepted: 10/13/2010] [Indexed: 11/09/2022]
|
46
|
Sharma S, Singh H, Harvey JN, Balint-Kurti GG. Design of an infrared laser pulse to control the multiphoton dissociation of the Fe-CO bond in CO-heme compounds. J Chem Phys 2010; 133:174103. [PMID: 21054002 DOI: 10.1063/1.3494543] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Optimal control theory is used to design a laser pulse for the multiphoton dissociation of the Fe-CO bond in the CO-heme compounds. The study uses a hexacoordinated iron-porphyrin-imidazole-CO complex in its ground electronic state as a model for CO liganded to the heme group. The potential energy and dipole moment surfaces for the interaction of the CO ligand with the heme group are calculated using density functional theory. Optimal control theory, combined with a time-dependent quantum dynamical treatment of the laser-molecule interaction, is then used to design a laser pulse capable of efficiently dissociating the CO-heme complex model. The genetic algorithm method is used within the mathematical framework of optimal control theory to perform the optimization process. This method provides good control over the parameters of the laser pulse, allowing optimized pulses with simple time and frequency structures to be designed. The dependence of photodissociation yield on the choice of initial vibrational state and of initial laser field parameters is also investigated. The current work uses a reduced dimensionality model in which only the Fe-C and C-O stretching coordinates are explicitly taken into account in the time-dependent quantum dynamical calculations. The limitations arising from this are discussed in Sec. IV.
Collapse
Affiliation(s)
- Sitansh Sharma
- Center for Computational Natural Sciences and Bioinformatics, International institute of Information Technology, Hyderabad 500032, India.
| | | | | | | |
Collapse
|
47
|
Cwiklik L, Devlin JP. Hindering of rotational motion of guest molecules in the Type I clathrate hydrate. Chem Phys Lett 2010. [DOI: 10.1016/j.cplett.2010.06.025] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
48
|
Wittich B, Deiters UK. Calculating Thermodynamic Properties of an Ionic Liquid with Monte Carlo Simulations with an Orthorhombic and a Cubic Simulation Box. J Phys Chem B 2010; 114:8954-60. [DOI: 10.1021/jp101676m] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Björn Wittich
- Department of Physical Chemistry, University of Cologne, Cologne, Germany
| | - Ulrich K. Deiters
- Department of Physical Chemistry, University of Cologne, Cologne, Germany
| |
Collapse
|
49
|
D’Abramo M, Di Nola A, Amadei A. Kinetics of Carbon Monoxide Migration and Binding in Solvated Myoglobin as Revealed by Molecular Dynamics Simulations and Quantum Mechanical Calculations. J Phys Chem B 2009; 113:16346-53. [DOI: 10.1021/jp903165p] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Marco D’Abramo
- Institut de Recerca Biomèdica, Parc Cientific de Barcelona Josep Samitier 1-5, Barcelona 08028 and Barcelona Supercomputing Center Jordi Girona 29, Barcelona 08034, Spain, Departament de Bioquimica, Facultat de Biologia, Universitat de Barcelona, Avgda Diagonal 647 Barcelona 08028, Spain, Dipartimento di Chimica, Università di Roma “La Sapienza”, P.le A. Moro 5 00185 Rome, Italy, and Departimento di Scienze e Tecnologie Chimiche, Università di Roma “Tor Vergata”, via della Ricerca Scientifica 00133 Rome,
| | - Alfredo Di Nola
- Institut de Recerca Biomèdica, Parc Cientific de Barcelona Josep Samitier 1-5, Barcelona 08028 and Barcelona Supercomputing Center Jordi Girona 29, Barcelona 08034, Spain, Departament de Bioquimica, Facultat de Biologia, Universitat de Barcelona, Avgda Diagonal 647 Barcelona 08028, Spain, Dipartimento di Chimica, Università di Roma “La Sapienza”, P.le A. Moro 5 00185 Rome, Italy, and Departimento di Scienze e Tecnologie Chimiche, Università di Roma “Tor Vergata”, via della Ricerca Scientifica 00133 Rome,
| | - Andrea Amadei
- Institut de Recerca Biomèdica, Parc Cientific de Barcelona Josep Samitier 1-5, Barcelona 08028 and Barcelona Supercomputing Center Jordi Girona 29, Barcelona 08034, Spain, Departament de Bioquimica, Facultat de Biologia, Universitat de Barcelona, Avgda Diagonal 647 Barcelona 08028, Spain, Dipartimento di Chimica, Università di Roma “La Sapienza”, P.le A. Moro 5 00185 Rome, Italy, and Departimento di Scienze e Tecnologie Chimiche, Università di Roma “Tor Vergata”, via della Ricerca Scientifica 00133 Rome,
| |
Collapse
|
50
|
Lutz S, Nienhaus K, Nienhaus GU, Meuwly M. Ligand Migration between Internal Docking Sites in Photodissociated Carbonmonoxy Neuroglobin. J Phys Chem B 2009; 113:15334-43. [DOI: 10.1021/jp905673p] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Stephan Lutz
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland; Institute of Biophysics, University of Ulm, Albert-Einstein-Allee 11, 89081 Ulm, Germany; Institute of Applied Physics and Center for Functional Nanostructures, Karlsruhe Institute of Technology, Karlsruhe, Germany; and Department of Physics, University of Illinois at Urbana−Champaign, 1110 West Green Street, Urbana, Illinois 61801
| | - Karin Nienhaus
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland; Institute of Biophysics, University of Ulm, Albert-Einstein-Allee 11, 89081 Ulm, Germany; Institute of Applied Physics and Center for Functional Nanostructures, Karlsruhe Institute of Technology, Karlsruhe, Germany; and Department of Physics, University of Illinois at Urbana−Champaign, 1110 West Green Street, Urbana, Illinois 61801
| | - G. Ulrich Nienhaus
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland; Institute of Biophysics, University of Ulm, Albert-Einstein-Allee 11, 89081 Ulm, Germany; Institute of Applied Physics and Center for Functional Nanostructures, Karlsruhe Institute of Technology, Karlsruhe, Germany; and Department of Physics, University of Illinois at Urbana−Champaign, 1110 West Green Street, Urbana, Illinois 61801
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland; Institute of Biophysics, University of Ulm, Albert-Einstein-Allee 11, 89081 Ulm, Germany; Institute of Applied Physics and Center for Functional Nanostructures, Karlsruhe Institute of Technology, Karlsruhe, Germany; and Department of Physics, University of Illinois at Urbana−Champaign, 1110 West Green Street, Urbana, Illinois 61801
| |
Collapse
|