1
|
Hao H, Ruiz Pestana L, Qian J, Liu M, Xu Q, Head‐Gordon T. Chemical transformations and transport phenomena at interfaces. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1639] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Hongxia Hao
- Kenneth S. Pitzer Theory Center and Department of Chemistry University of California Berkeley California USA
- Chemical Sciences Division Lawrence Berkeley National Laboratory Berkeley California USA
| | - Luis Ruiz Pestana
- Department of Civil and Architectural Engineering University of Miami Coral Gables Florida USA
| | - Jin Qian
- Chemical Sciences Division Lawrence Berkeley National Laboratory Berkeley California USA
| | - Meili Liu
- Department of Civil and Architectural Engineering University of Miami Coral Gables Florida USA
| | - Qiang Xu
- Chemical Sciences Division Lawrence Berkeley National Laboratory Berkeley California USA
| | - Teresa Head‐Gordon
- Kenneth S. Pitzer Theory Center and Department of Chemistry University of California Berkeley California USA
- Chemical Sciences Division Lawrence Berkeley National Laboratory Berkeley California USA
- Department of Bioengineering and Chemical and Biomolecular Engineering University of California Berkeley California USA
| |
Collapse
|
2
|
Water structure in glycerol: Spectroscopic and computer simulation investigation of hydrogen bonding and water clustering. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.118916] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
3
|
Konermann L, Kim S. Grotthuss Molecular Dynamics Simulations for Modeling Proton Hopping in Electrosprayed Water Droplets. J Chem Theory Comput 2022; 18:3781-3794. [PMID: 35544700 DOI: 10.1021/acs.jctc.2c00001] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Excess protons in water exhibit unique transport properties because they can rapidly hop along H-bonded water wires. Considerable progress has been made in unraveling this Grotthuss diffusion mechanism using quantum mechanical-based computational techniques. Unfortunately, high computational cost tends to restrict those techniques to small systems and short times. Molecular dynamics (MD) simulations can be applied to much larger systems and longer time windows. However, standard MD methods do not permit the dissociation/formation of covalent bonds, such that Grotthuss diffusion cannot be captured. Here, we bridge this gap by combining atomistic MD simulations (using Gromacs and TIP4P/2005 water) with proton hopping. Excess protons are modeled as hydronium ions that undergo H3O+ + H2O → H2O + H3O+ transitions. In accordance with ab initio MD data, these Grotthuss hopping events are executed in "bursts" with quasi-instantaneous hopping across one or more waters. The bursts are separated by regular MD periods during which H3O+ ions undergo Brownian diffusion. The resulting proton diffusion coefficient agrees with the literature value. We apply this Grotthuss MD technique to highly charged water droplets that are in a size regime encountered during electrospray ionization (5 nm radius, ∼17,000 H2O). The droplets undergo rapid solvent evaporation and occasional H3O+ ejection, keeping them at ca. 81% of the Rayleigh limit. The simulated behavior is consistent with phase Doppler anemometry data. The Grotthuss MD technique developed here should be useful for modeling the behavior of various proton-containing systems that are too large for high-level computational approaches. In particular, we envision future applications related to electrospray processes, where earlier simulations used metal cations while in reality excess protons dominate.
Collapse
Affiliation(s)
- Lars Konermann
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| | - Scott Kim
- Department of Chemistry, The University of Western Ontario, London, Ontario N6A 5B7, Canada
| |
Collapse
|
4
|
Shangguan J, Hensley AJ, Morgenstern L, Li Z, McEwen JS, Ma W, Cathy China YH. Brønsted Acidity of H-adatoms at Protic Solvent-Transition Metal Interfaces and its Kinetic Consequences in Electrophilic Addition Reactions. J Catal 2021. [DOI: 10.1016/j.jcat.2021.12.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
5
|
Abstract
Machine learning (ML) techniques applied to chemical reactions have a long history. The present contribution discusses applications ranging from small molecule reaction dynamics to computational platforms for reaction planning. ML-based techniques can be particularly relevant for problems involving both computation and experiments. For one, Bayesian inference is a powerful approach to develop models consistent with knowledge from experiments. Second, ML-based methods can also be used to handle problems that are formally intractable using conventional approaches, such as exhaustive characterization of state-to-state information in reactive collisions. Finally, the explicit simulation of reactive networks as they occur in combustion has become possible using machine-learned neural network potentials. This review provides an overview of the questions that can and have been addressed using machine learning techniques, and an outlook discusses challenges in this diverse and stimulating field. It is concluded that ML applied to chemistry problems as practiced and conceived today has the potential to transform the way with which the field approaches problems involving chemical reactions, in both research and academic teaching.
Collapse
Affiliation(s)
- Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland.,Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
6
|
Liu P, Yan Z, Mei D. Insights into protonation for cyclohexanol/water mixtures at the zeolitic Brønsted acid site. Phys Chem Chem Phys 2021; 23:10395-10401. [PMID: 33889887 DOI: 10.1039/d0cp06523d] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Proton transfer from Brønsted acid sites (BASs) to alcohol molecules ignites the acid-catalyzed alcohol dehydration reactions. For aqueous phase dehydration reactions in zeolites, the coexisting water molecules around BASs in the zeolite pores significantly affect the alcohol dehydration activity. In the present work, proton transfer processes among the BASs of H-BEA zeolites, the adsorbed cyclohexanol and surrounding water clusters with different sizes up to 8 water molecules were investigated using ab initio molecular dynamics (AIMD) simulations combined with the multiple-walker well-tempered metadynamics algorithm. The plausible proton locations and proton transfer processes were characterized using two/three-dimensional free energy landscapes. The strong proton affinity makes the protonated cyclohexanol stable species until a water trimer is formed. The proton either is shared between protonated cyclohexanol and the water trimer or remains with the water trimer (H7O3+). With a further increase in water concentrations, the proton prefers to remain with the water clusters.
Collapse
Affiliation(s)
- Peng Liu
- School of Chemistry and Chemical Engineering, State Key Laboratory of Separation Membranes and Membrane Processes, Tiangong University, Tianjin 300387, P. R. China.
| | - Zhenxin Yan
- School of Chemistry and Chemical Engineering, State Key Laboratory of Separation Membranes and Membrane Processes, Tiangong University, Tianjin 300387, P. R. China.
| | - Donghai Mei
- School of Chemistry and Chemical Engineering, State Key Laboratory of Separation Membranes and Membrane Processes, Tiangong University, Tianjin 300387, P. R. China. and School of Environmental Science and Engineering, Tiangong University, Tianjin 300387, P. R. China
| |
Collapse
|
7
|
Protonation Dynamics in the K-Channel of Cytochrome c Oxidase Estimated from Molecular Dynamics Simulations. Processes (Basel) 2021. [DOI: 10.3390/pr9020265] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
Proton transfer reactions are one of the most fundamental processes in biochemistry. We present a simplistic approach for estimating proton transfer probabilities in a membrane protein, cytochrome c oxidase. We combine short molecular dynamics simulations at discrete protonation states with a Monte Carlo approach to exchange between those states. Requesting for a proton transfer the existence of a hydrogen-bonded connection between the two source and target residues of the exchange, restricts the acceptance of transfers to only those in which a proton-relay is possible. Together with an analysis of the hydrogen-bonded connectivity in one of the proton-conducting channels of cytochrome c oxidase, this approach gives insight into the protonation dynamics of the hydrogen-bonded networks. The connectivity and directionality of the networks are coupled to the conformation of an important protein residue in the channel, K362, rendering proton transfer in the entire channel feasible in only one of the two major conformations. Proton transport in the channel can thus be regulated by K362 not only through its possible role as a proton carrier itself, but also by allowing or preventing proton transport via water residues.
Collapse
|
8
|
Synergistic Computational Modeling Approaches as Team Players in the Game of Solubility Predictions. J Pharm Sci 2020; 110:22-34. [PMID: 33217423 DOI: 10.1016/j.xphs.2020.10.068] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2020] [Revised: 10/23/2020] [Accepted: 10/28/2020] [Indexed: 11/23/2022]
Abstract
Several approaches to predict and model drug solubility have been used in the drug discovery and development processes during the last decades. Each of these approaches have their own benefits and place, and are typically used as standalone approaches rather than in concert. The synergistic effects of these are often overlooked, partly due to the need of computational experts to perform the modeling and simulations as well as analyzing the data obtained. Here we provide our views on how these different approaches can be used to retrieve more information on drug solubility, ranging from multivariate data analysis over thermodynamic cycle modeling to molecular dynamics simulations. We are discussing aqueous solubility as well as solubility in more complex mixed solvents and media with colloidal structures present. We conclude that the field of computational pharmaceutics is in its early days but with a bright future ahead. However, education of computational formulators with broad knowledge of modeling and simulation approaches is imperative if computational pharmaceutics is to reach its full potential.
Collapse
|
9
|
Leven I, Hao H, Das AK, Head-Gordon T. A Reactive Force Field with Coarse-Grained Electrons for Liquid Water. J Phys Chem Lett 2020; 11:9240-9247. [PMID: 33073998 DOI: 10.1021/acs.jpclett.0c02516] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Nonreactive force fields are defined by perturbations of electron density that are relatively small, whereas chemical reactivity involves wholesale electronic rearrangements that make and break bonds. Thus, reactive force fields are incredibly difficult to develop compared to nonreactive force fields, yet at the same time, they fill a critical need when ab initio molecular dynamics methods are not affordable. We introduce a new reactive force field model for water that combines modified nonbonded terms of the ReaxFF model and its embedding in the electrostatic interactions described by our recently introduced coarse-grained electron model (C-GeM). The ReaxFF/C-GeM force field is characterized for many energetic and dissociative water properties for water clusters, structure, and dynamical properties under ambient conditions in the condensed phase, as well as the temperature dependence of density and water diffusion, with very good agreement with experiment. The ReaxFF/C-GeM force field should be more transferable and more broadly applicable to a range of reactive systems involving both proton and electron transfer in the condensed phase.
Collapse
Affiliation(s)
- Itai Leven
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Hongxia Hao
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | | | - Teresa Head-Gordon
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
10
|
Carpenter WB, Yu Q, Hack JH, Dereka B, Bowman JM, Tokmakoff A. Decoding the 2D IR spectrum of the aqueous proton with high-level VSCF/VCI calculations. J Chem Phys 2020; 153:124506. [PMID: 33003749 DOI: 10.1063/5.0020279] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
The aqueous proton is a common and long-studied species in chemistry, yet there is currently intense interest devoted to understanding its hydration structure and transport dynamics. Typically described in terms of two limiting structures observed in gas-phase clusters, the Zundel H5O2 + and Eigen H9O4 + ions, the aqueous structure is less clear due to the heterogeneity of hydrogen bonding environments and room-temperature structural fluctuations in water. The linear infrared (IR) spectrum, which reports on structural configurations, is challenging to interpret because it appears as a continuum of absorption, and the underlying vibrational modes are strongly anharmonically coupled to each other. Recent two-dimensional IR (2D IR) experiments presented strong evidence for asymmetric Zundel-like motifs in solution, but true structure-spectrum correlations are missing and complicated by the anharmonicity of the system. In this study, we employ high-level vibrational self-consistent field/virtual state configuration interaction calculations to demonstrate that the 2D IR spectrum reports on a broad distribution of geometric configurations of the aqueous proton. We find that the diagonal 2D IR spectrum around 1200 cm-1 is dominated by the proton stretch vibrations of Zundel-like and intermediate geometries, broadened by the heterogeneity of aqueous configurations. There is a wide distribution of multidimensional potential shapes for the proton stretching vibration with varying degrees of potential asymmetry and confinement. Finally, we find specific cross peak patterns due to aqueous Zundel-like species. These studies provide clarity on highly debated spectral assignments and stringent spectroscopic benchmarks for future simulations.
Collapse
Affiliation(s)
- William B Carpenter
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Qi Yu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - John H Hack
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Bogdan Dereka
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Andrei Tokmakoff
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
11
|
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
|
12
|
Zeng Y, Li A, Yan T. Hydrogen Bond Dynamics in the Solvation Shell on Proton Transfer in Aqueous Solution. J Phys Chem B 2020; 124:1817-1823. [DOI: 10.1021/acs.jpcb.0c00990] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Affiliation(s)
- Yonghui Zeng
- Institute of New Energy Material Chemistry, School of Materials Science and Engineering, National Institute for Advanced Materials, Nankai University, Tianjin 300350, China
| | - Ailin Li
- College of Science, Civil Aviation University of China, Tianjin 300300, China
| | - Tianying Yan
- Institute of New Energy Material Chemistry, School of Materials Science and Engineering, National Institute for Advanced Materials, Nankai University, Tianjin 300350, China
| |
Collapse
|