1
|
France-Lanord A, Vroylandt H, Salanne M, Rotenberg B, Saitta AM, Pietrucci F. Data-Driven Path Collective Variables. J Chem Theory Comput 2024; 20:3069-3084. [PMID: 38619076 DOI: 10.1021/acs.jctc.4c00123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/16/2024]
Abstract
Identifying optimal collective variables to model transformations using atomic-scale simulations is a long-standing challenge. We propose a new method for the generation, optimization, and comparison of collective variables that can be thought of as a data-driven generalization of the path collective variable concept. It consists of a kernel ridge regression of the committor probability, which encodes a transformation's progress. The resulting collective variable is one-dimensional, interpretable, and differentiable, making it appropriate for enhanced sampling simulations requiring biasing. We demonstrate the validity of the method on two different applications: a precipitation model and the association of Li+ and F- in water. For the former, we show that global descriptors such as the permutation invariant vector allow reaching an accuracy far from the one achieved via simpler, more intuitive variables. For the latter, we show that information correlated with the transformation mechanism is contained in the first solvation shell only and that inertial effects prevent the derivation of optimal collective variables from the atomic positions only.
Collapse
Affiliation(s)
- Arthur France-Lanord
- Institut des Sciences du Calcul et des Données, ISCD, Sorbonne Université, F-75005 Paris, France
- Muséum National d'Histoire Naturelle, UMR CNRS 7590, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, IMPMC, Sorbonne Université, F-75005 Paris, France
| | - Hadrien Vroylandt
- Institut des Sciences du Calcul et des Données, ISCD, Sorbonne Université, F-75005 Paris, France
| | - Mathieu Salanne
- Physicochimie des Électrolytes et Nanosystèmes Interfaciaux, Sorbonne Université, CNRS, 4 Place Jussieu, F-75005 Paris, France
- Institut Universitaire de France (IUF), 75231 Paris, France
| | - Benjamin Rotenberg
- Physicochimie des Électrolytes et Nanosystèmes Interfaciaux, Sorbonne Université, CNRS, 4 Place Jussieu, F-75005 Paris, France
| | - A Marco Saitta
- Muséum National d'Histoire Naturelle, UMR CNRS 7590, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, IMPMC, Sorbonne Université, F-75005 Paris, France
| | - Fabio Pietrucci
- Muséum National d'Histoire Naturelle, UMR CNRS 7590, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, IMPMC, Sorbonne Université, F-75005 Paris, France
| |
Collapse
|
2
|
Naleem N, Abreu CRA, Warmuz K, Tong M, Kirmizialtin S, Tuckerman ME. An exploration of machine learning models for the determination of reaction coordinates associated with conformational transitions. J Chem Phys 2023; 159:034102. [PMID: 37458344 DOI: 10.1063/5.0147597] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2023] [Accepted: 06/23/2023] [Indexed: 07/20/2023] Open
Abstract
Determining collective variables (CVs) for conformational transitions is crucial to understanding their dynamics and targeting them in enhanced sampling simulations. Often, CVs are proposed based on intuition or prior knowledge of a system. However, the problem of systematically determining a proper reaction coordinate (RC) for a specific process in terms of a set of putative CVs can be achieved using committor analysis (CA). Identifying essential degrees of freedom that govern such transitions using CA remains elusive because of the high dimensionality of the conformational space. Various schemes exist to leverage the power of machine learning (ML) to extract an RC from CA. Here, we extend these studies and compare the ability of 17 different ML schemes to identify accurate RCs associated with conformational transitions. We tested these methods on an alanine dipeptide in vacuum and on a sarcosine dipeptoid in an implicit solvent. Our comparison revealed that the light gradient boosting machine method outperforms other methods. In order to extract key features from the models, we employed Shapley Additive exPlanations analysis and compared its interpretation with the "feature importance" approach. For the alanine dipeptide, our methodology identifies ϕ and θ dihedrals as essential degrees of freedom in the C7ax to C7eq transition. For the sarcosine dipeptoid system, the dihedrals ψ and ω are the most important for the cisαD to transαD transition. We further argue that analysis of the full dynamical pathway, and not just endpoint states, is essential for identifying key degrees of freedom governing transitions.
Collapse
Affiliation(s)
- Nawavi Naleem
- Chemistry Program, Science Division, New York University, Abu Dhabi, UAE
| | - Charlles R A Abreu
- Chemical Engineering Department, Escola de Química, Universidade Federal do Rio de Janeiro, 21941-909 Rio de Janeiro, RJ, Brazil
| | - Krzysztof Warmuz
- Computer Science Program, Science Division, New York University, Abu Dhabi, UAE
| | - Muchen Tong
- Department of Chemistry, New York University (NYU), New York, New York 10003, USA
| | - Serdal Kirmizialtin
- Chemistry Program, Science Division, New York University, Abu Dhabi, UAE
- Department of Chemistry, New York University (NYU), New York, New York 10003, USA
- Center for Smart Engineering Materials, New York University, Abu Dhabi, UAE
| | - Mark E Tuckerman
- Department of Chemistry, New York University (NYU), New York, New York 10003, USA
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Rd. North, Shanghai 200062, China
- Simons Center for Computational Physical Chemistry at New York University, New York, New York 10003, USA
| |
Collapse
|
3
|
Burgin T, Ellis S, Mayes HB. ATESA: An Automated Aimless Shooting Workflow. J Chem Theory Comput 2023; 19:235-244. [PMID: 36520006 DOI: 10.1021/acs.jctc.2c00543] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Transition path sampling methods are powerful tools for studying the dynamics of rare events in molecular simulations. However, these methods are generally restricted to experts with the knowledge and resources to properly set up and analyze the often hundreds of thousands of simulations that constitute a complete study. Aimless Transition Ensemble Sampling and Analysis (ATESA) is a new open-source software program written in Python that automates a full transition path sampling workflow based on the aimless shooting algorithm, streamlining the process and reducing the barrier to use for researchers new to this approach. This introduction to ATESA includes a demonstration of a complete transition path sampling process flow for an example reaction, including finding an initial transition state, sampling with aimless shooting, building a reaction coordinate with inertial likelihood maximization, verifying that coordinate with committor analysis, and measuring the reaction energy profile with umbrella sampling. We also describe our implementation of a termination criterion for aimless shooting based on the Godambe information calculated during model building with likelihood maximization as well as a novel approach to constraining simulations to the desired rare event pathway during umbrella sampling.
Collapse
Affiliation(s)
- Tucker Burgin
- Department of Chemical Engineering, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Samuel Ellis
- The Molecular Sciences Software Institute, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24060, United States.,Department of Chemistry, Virginia Tech University, Blacksburg, Virginia 24061, United States
| | - Heather B Mayes
- Department of Chemical Engineering, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
4
|
Roy S, Bocharova V, Stack AG, Bryantsev VS. Nucleation Rate Theory for Coordination Number: Elucidating Water-Mediated Formation of a Zigzag Na 2SO 4 Morphology. ACS APPLIED MATERIALS & INTERFACES 2022; 14:53213-53227. [PMID: 36395432 DOI: 10.1021/acsami.2c17475] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
Predicting and controlling nanostructure formation during nucleation can pave the way to synthesizing novel energy materials via crystallization. However, such control over nucleation and crystallization remains challenging due to an inadequate understanding of critical factors that govern evolving atomistic structures and dynamics. Herein, we utilize coordination number as a reaction coordinate and rate theory to investigate how sodium sulfate, commonly known as a phase-change energy material, nucleates in a supersaturated aqueous solution. In conjunction with ab initio and force field-based molecular dynamics simulation, the rate theoretical analysis reveals that sodium sulfate from an initially dissolved metastable state transits to a heterogeneous mixture of prenucleated clusters and finally to a large cylindrical zigzag morphology. Measurements of Raman spectra and their ab initio modeling confirm that this nucleated morphology contains a few waters for every sulfate. Rate processes such as solvent exchange and desolvation exhibit high sensitivity to the evolving prenucleation/nucleation structures, providing a means to distinguish between critical nucleation precursors. Desolvation and forming the first-shell interionic coordination structure via monomer-by-monomer addition around sulfates are found to explain the formation of large nuclei. Thus, a detailed understanding of the step-by-step structure formation across scales has been achieved. This can be leveraged to predict nucleation-related structures and dynamics and potentially control the synthesis of novel phase-change materials for energy applications.
Collapse
Affiliation(s)
- Santanu Roy
- Chemical Sciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee37830, United States
| | - Vera Bocharova
- Chemical Sciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee37830, United States
| | - Andrew G Stack
- Chemical Sciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee37830, United States
| | - Vyacheslav S Bryantsev
- Chemical Sciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee37830, United States
| |
Collapse
|
5
|
Lei YK, Zhang Z, Han X, Yang YI, Zhang J, Gao YQ. Locating Transition Zone in Phase Space. J Chem Theory Comput 2022; 18:6124-6133. [PMID: 36135927 DOI: 10.1021/acs.jctc.2c00385] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Understanding the reaction mechanism is required for better control of chemical reactions and is usually achieved by locating transition states (TSs) along a proper one-dimensional coordinate called reaction coordinate (RC). The identification of RC can be very difficult for high-dimensional realistic systems. A number of methods have been proposed to tackle this problem. A machine learning method is developed here to incorporate the influence of velocity on the reaction process. The method is also free of the unbalanced label problem resulting from the rather low fraction of configurations near the TS and can be easily extended to large systems. It locates the transition zone in the phase space and defines the dividing surface with a high transmission coefficient. Moreover, considering that the reaction environment can not only change the reaction path but also activate the reactive mode through energy transfer, we devise two measures to quantify the influence of these two factors on the reaction process and find that solvents can assist the reaction by directly doing work along the reactive mode. Not surprisingly, there is a positive correlation between the efficiency of energy transfer into the reactive mode and the reaction rate.
Collapse
Affiliation(s)
- Yao-Kun Lei
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, 518055 Shenzhen, China.,Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, 100871 Beijing, China
| | - Zhen Zhang
- School of Physics and Technology, Tangshan Normal University, 063000 Tangshan, China
| | - Xu Han
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, 518055 Shenzhen, China.,Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, 100871 Beijing, China
| | - Yi Isaac Yang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, 518055 Shenzhen, China
| | - Jun Zhang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, 518055 Shenzhen, China
| | - Yi Qin Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, 518055 Shenzhen, China.,Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, 100871 Beijing, China.,Biomedical Pioneering Innovation Center, Peking University, 100871 Beijing, China
| |
Collapse
|
6
|
Kikutsuji T, Mori Y, Okazaki KI, Mori T, Kim K, Matubayasi N. Explaining reaction coordinates of alanine dipeptide isomerization obtained from deep neural networks using Explainable Artificial Intelligence (XAI). J Chem Phys 2022; 156:154108. [DOI: 10.1063/5.0087310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A method for obtaining appropriate reaction coordinates is required to identify transition states distinguishing product and reactant in complex molecular systems. Recently, abundant research has been devoted to obtaining reaction coordinates using artificial neural networks from deep learning literature, where many collective variables are typically utilized in the input layer. However, it is difficult to explain the details of which collective variables contribute to the predicted reaction coordinates owing to the complexity of the nonlinear functions in deep neural networks. To overcome this limitation, we used Explainable Artificial Intelligence (XAI) methods of the Local Interpretable Model-agnostic Explanation (LIME) and the game theory-based framework known as Shapley Additive exPlanations (SHAP). We demonstrated that XAI enables us to obtain the degree of contribution of each collective variable to reaction coordinates that is determined by nonlinear regressions with deep learning for the committor of the alanine dipeptide isomerization in vacuum. In particular, both LIME and SHAP provide important features to the predicted reaction coordinates, which are characterized by appropriate dihedral angles consistent with those previously reported from the committor test analysis. The present study offers an AI-aided framework to explain the appropriate reaction coordinates, which acquires considerable significance when the number of degrees of freedom increases.
Collapse
Affiliation(s)
| | | | - Kei-ichi Okazaki
- Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, Japan
| | - Toshifumi Mori
- Kyushu University Institute for Materials Chemistry and Engineering, Japan
| | - Kang Kim
- Graduate School of Engineering Science, Osaka University - Toyonaka Campus, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Japan
| |
Collapse
|
7
|
Paul TK, Taraphder S. Nonlinear Reaction Coordinate of an Enzyme Catalyzed Proton Transfer Reaction. J Phys Chem B 2022; 126:1413-1425. [PMID: 35138854 DOI: 10.1021/acs.jpcb.1c08760] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
We present an in-depth study on the theoretical calculation of an optimum reaction coordinate as a linear or nonlinear combination of important collective variables (CVs) sampled from an ensemble of reactive transition paths for an intramolecular proton transfer reaction catalyzed by the enzyme human carbonic anhydrase (HCA) II. The linear models are optimized by likelihood maximization for a given number of CVs. The nonlinear models are based on an artificial neural network with the same number of CVs and optimized by minimizing the root-mean-square error in comparison to a training set of committor estimators generated for the given transition. The nonlinear reaction coordinate thus obtained yields the free energy of activation and rate constant as 9.46 kcal mol-1 and 1.25 × 106 s-1, respectively. These estimates are found to be in quantitative agreement with the known experimental results. We have also used an extended autoencoder model to show that a similar analysis can be carried out using a single CV only. The resultant free energies and kinetics of the reaction slightly overestimate the experimental data. The implications of these results are discussed using a detailed microkinetic scheme of the proton transfer reaction catalyzed by HCA II.
Collapse
Affiliation(s)
- Tanmoy Kumar Paul
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Srabani Taraphder
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| |
Collapse
|
8
|
Bolhuis PG, Swenson DWH. Transition Path Sampling as Markov Chain Monte Carlo of Trajectories: Recent Algorithms, Software, Applications, and Future Outlook. ADVANCED THEORY AND SIMULATIONS 2021. [DOI: 10.1002/adts.202000237] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Peter G. Bolhuis
- Amsterdam Center for Multiscale Modeling van 't Hoff Institute for Molecular Sciences University of Amsterdam PO Box 94157 1090 GD Amsterdam The Netherlands
| | - David W. H. Swenson
- Centre Blaise Pascal Ecole Normale Superieure 46, allée d'Italie 69364 Lyon Cedex 07 France
| |
Collapse
|
9
|
Quiton SJ, Mitra U, Mallikarjun Sharada S. A matrix completion algorithm to recover modes orthogonal to the minimum energy path in chemical reactions. J Chem Phys 2020; 153:054122. [PMID: 32770897 DOI: 10.1063/5.0018326] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Structured statistical methods are promising for recovering or completing information from noisy and incomplete data with high fidelity. In particular, matrix completion exploits underlying structural properties such as rank or sparsity. Our objective is to employ matrix completion to reduce computational effort associated with the calculation of multiple quantum chemical Hessians, which are necessary for identification of temperature-dependent free energy maxima under canonical variational transition state theory (VTST). We demonstrate proof-of-principle of an algebraic variety-based matrix completion method for recovering missing elements in a matrix of transverse Hessian eigenvalues constituting the minimum energy path (MEP) of a reaction. The algorithm, named harmonic variety-based matrix completion (HVMC), utilizes the fact that the points lying on the MEP of a reaction step constitute an algebraic variety in the reaction path Hamiltonian representation. We demonstrate that, with as low as 30% random sampling of matrix elements for the largest system in our test set (46 atoms), the complete matrix of eigenvalues can be recovered. We further establish algorithm performance for VTST rate calculations by quantifying errors in zero-point energies and vibrational free energies. Motivated by this success, we outline next steps toward developing a practical HVMC algorithm, which utilizes a gradient-based sampling protocol for low-cost VTST rate computations.
Collapse
Affiliation(s)
- Stephen Jon Quiton
- Mork Family Department of Chemical Engineering and Materials Science, University of Southern California, Los Angeles, California 90089, USA
| | - Urbashi Mitra
- Ming Hsieh Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, California 90089, USA
| | - Shaama Mallikarjun Sharada
- Mork Family Department of Chemical Engineering and Materials Science, University of Southern California, Los Angeles, California 90089, USA
| |
Collapse
|
10
|
Mori Y, Okazaki KI, Mori T, Kim K, Matubayasi N. Learning reaction coordinates via cross-entropy minimization: Application to alanine dipeptide. J Chem Phys 2020; 153:054115. [DOI: 10.1063/5.0009066] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Affiliation(s)
- Yusuke Mori
- Division of Chemical Engineering, Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | | | - Toshifumi Mori
- Institute for Molecular Science, Okazaki, Aichi 444-8585, Japan
- The Graduate University for Advanced Studies, Okazaki, Aichi 444-8585, Japan
| | - Kang Kim
- Division of Chemical Engineering, Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
- Institute for Molecular Science, Okazaki, Aichi 444-8585, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| |
Collapse
|
11
|
Leitold C, Mundy CJ, Baer MD, Schenter GK, Peters B. Solvent reaction coordinate for an SN2 reaction. J Chem Phys 2020; 153:024103. [DOI: 10.1063/5.0002766] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Christian Leitold
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Christopher J. Mundy
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - Marcel D. Baer
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - Gregory K. Schenter
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - Baron Peters
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
- Department of Chemistry and Biochemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
12
|
Roy S, Schenter GK, Napoli JA, Baer MD, Markland TE, Mundy CJ. Resolving Heterogeneous Dynamics of Excess Protons in Aqueous Solution with Rate Theory. J Phys Chem B 2020; 124:5665-5675. [DOI: 10.1021/acs.jpcb.0c02649] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Santanu Roy
- Chemical Sciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee 37830, United States
| | - Gregory K. Schenter
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Blvd., Richland, Washington 99352, United States
| | - Joseph A. Napoli
- Department of Chemistry, Stanford University, 333 Campus Drive, Stanford, California 94305, United States
| | - Marcel D. Baer
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Blvd., Richland, Washington 99352, United States
| | - Thomas E. Markland
- Department of Chemistry, Stanford University, 333 Campus Drive, Stanford, California 94305, United States
| | - Christopher J. Mundy
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Blvd., Richland, Washington 99352, United States
- Affiliate Professor, Department of Chemical Engineering, University of Washington, Seattle, United States
| |
Collapse
|
13
|
Roy S, Wu L, Goverapet Srinivasan S, Stack AG, Navrotsky A, Bryantsev VS. Hydration structure and water exchange kinetics at xenotime-water interfaces: implications for rare earth minerals separation. Phys Chem Chem Phys 2020; 22:7719-7727. [PMID: 32215419 DOI: 10.1039/d0cp00087f] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Hydration of surface ions gives rise to structural heterogeneity and variable exchange kinetics of water at complex mineral-water interfaces. Here, we employ ab initio molecular dynamics (AIMD) simulations and water adsorption calorimetry to examine the aqueous interfaces of xenotime, a phosphate mineral that contains predominantly Y3+ and heavy rare earth elements. Consistent with natural crystal morphology, xenotime is predicted to have a tetragonal prismatic shape, dominated by the {100} surface. Hydration of this surface induces multilayer interfacial water structures with distinct OH orientations, which agrees with recent crystal truncation rod measurements. The exchange kinetics between two adjacent water layers exhibits a wide range of underlying timescales (5-180 picoseconds), dictated by ion-water electrostatics. Adsorption of a bidentate hydroxamate ligand reveals that {100} xenotime surface can only accommodate monodentate coordination with water exchange kinetics strongly depending on specific ligand orientation, prompting us to reconsider traditional strategies for selective separation of rare-earth minerals.
Collapse
Affiliation(s)
- Santanu Roy
- Chemical Sciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Rd., Oak Ridge, TN 37830, USA.
| | - Lili Wu
- Peter A. Rock Thermochemistry Laboratory and NEAT ORU, University of California Davis, 1 Shields Avenue, Davis, CA 95616, USA
| | | | - Andrew G Stack
- Chemical Sciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Rd., Oak Ridge, TN 37830, USA.
| | - Alexandra Navrotsky
- Peter A. Rock Thermochemistry Laboratory and NEAT ORU, University of California Davis, 1 Shields Avenue, Davis, CA 95616, USA
| | - Vyacheslav S Bryantsev
- Chemical Sciences Division, Oak Ridge National Laboratory, 1 Bethel Valley Rd., Oak Ridge, TN 37830, USA.
| |
Collapse
|
14
|
Roy S, Wu F, Wang H, Ivanov AS, Sharma S, Halstenberg P, Gill SK, Milinda Abeykoon AM, Kwon G, Topsakal M, Layne B, Sasaki K, Zhang Y, Mahurin SM, Dai S, Margulis CJ, Maginn EJ, Bryantsev VS. Structure and dynamics of the molten alkali-chloride salts from an X-ray, simulation, and rate theory perspective. Phys Chem Chem Phys 2020; 22:22900-22917. [DOI: 10.1039/d0cp03672b] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Molten salts are of great interest as alternative solvents, electrolytes, and heat transfer fluids in many emerging technologies.
Collapse
Affiliation(s)
- Santanu Roy
- Chemical Sciences Division
- Oak Ridge National Laboratory
- Oak Ridge
- USA
| | - Fei Wu
- Department of Chemistry
- The University of Iowa
- USA
| | - Haimeng Wang
- Department of Chemical and Biomolecular Engineering
- University of Notre Dame
- USA
| | | | | | - Phillip Halstenberg
- Chemical Sciences Division
- Oak Ridge National Laboratory
- Oak Ridge
- USA
- Department of Chemistry
| | - Simerjeet K Gill
- Nuclear Science and Technology Department
- Brookhaven National Lab
- USA
| | | | - Gihan Kwon
- National Synchrotron Light Source II (NSLS-II)
- Brookhaven National Lab
- USA
| | - Mehmet Topsakal
- Nuclear Science and Technology Department
- Brookhaven National Lab
- USA
| | - Bobby Layne
- Chemistry Division
- Brookhaven National Lab
- USA
| | | | - Yong Zhang
- Department of Chemical and Biomolecular Engineering
- University of Notre Dame
- USA
| | | | - Sheng Dai
- Chemical Sciences Division
- Oak Ridge National Laboratory
- Oak Ridge
- USA
- Department of Chemistry
| | | | - Edward J. Maginn
- Department of Chemical and Biomolecular Engineering
- University of Notre Dame
- USA
| | | |
Collapse
|
15
|
Yappert R, Kamat K, Peters B. The overdamped transmission coefficient: Recovering the true mean first passage time from an inaccurate reaction coordinate. J Chem Phys 2019; 151:184108. [DOI: 10.1063/1.5117237] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Affiliation(s)
- Ryan Yappert
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Kartik Kamat
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
| | - Baron Peters
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
16
|
Brotzakis ZF, Bolhuis PG. Approximating free energy and committor landscapes in standard transition path sampling using virtual interface exchange. J Chem Phys 2019; 151:174111. [DOI: 10.1063/1.5119252] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Z. Faidon Brotzakis
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Peter G. Bolhuis
- Van’t Hoff Institute for Molecular Sciences, University of Amsterdam, P.O. Box 94157, 1090 GD Amsterdam, The Netherlands
| |
Collapse
|
17
|
Paul S, Nair NN, Vashisth H. Phase space and collective variable based simulation methods for studies of rare events. MOLECULAR SIMULATION 2019. [DOI: 10.1080/08927022.2019.1634268] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Sanjib Paul
- Department of Chemical Engineering, University of New Hampshire, Durham, NH, USA
| | - Nisanth N. Nair
- Department of Chemistry, Indian Institute of Technology, Kanpur, India
| | - Harish Vashisth
- Department of Chemical Engineering, University of New Hampshire, Durham, NH, USA
| |
Collapse
|
18
|
Paul S, Paul TK, Taraphder S. Orthogonal order parameters to model the reaction coordinate of an enzyme catalyzed reaction. J Mol Graph Model 2019; 90:18-32. [PMID: 30959266 DOI: 10.1016/j.jmgm.2019.03.023] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2018] [Revised: 03/01/2019] [Accepted: 03/26/2019] [Indexed: 01/03/2023]
Abstract
The choice of suitable collective variables in formulating an optimal reaction coordinate is a challenging task for activated transitions between a pair of stable states especially when dealing with biochemical changes such as enzyme catalyzed reactions. A detailed benchmarking study is carried out on the choice of collective variables that can distinguish between the stable states unambiguously. We specifically address the issue if these variables may be directly used to model the optimal reaction coordinate, or if it would be better to use their orthogonalized counterparts. The proposed computational scheme is applied to the rate determining intramolecular proton transfer step in the enzyme human carbonic anhydrase II. The optimum reaction coordinate is determined with and without orthogonalization of the collective variables pertinent to a key conformational fluctuation and the actual proton transfer step at the active site of the enzyme. Suitability of the predicted reaction coordinates in different processes is examined in terms of the free energy profile projected along the reaction coordinate, the rate constant of transition and the underlying molecular mechanism of barrier crossing. Our results indicate that a better agreement with earlier simulation and experimental data is obtained when the orthogonalized collective variables are used to model the reaction coordinate.
Collapse
Affiliation(s)
- Sanjib Paul
- Department of Chemistry, Indian Institute of Technology, Kharagpur, 721302, India
| | - Tanmoy Kumar Paul
- Department of Chemistry, Indian Institute of Technology, Kharagpur, 721302, India
| | - Srabani Taraphder
- Department of Chemistry, Indian Institute of Technology, Kharagpur, 721302, India.
| |
Collapse
|
19
|
Burgin T, Mayes HB. Mechanism of oligosaccharide synthesis via a mutant GH29 fucosidase. REACT CHEM ENG 2019. [DOI: 10.1039/c8re00240a] [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/21/2022]
Abstract
First unbiased transition path sampling study of a glycosynthase enzyme reveals single-step mechanism with oxocarbenium-like transition state.
Collapse
Affiliation(s)
- Tucker Burgin
- University of Michigan Department of Chemical Engineering
- Ann Arbor
- USA
| | - Heather B. Mayes
- University of Michigan Department of Chemical Engineering
- Ann Arbor
- USA
| |
Collapse
|
20
|
Roy S, Bryantsev VS. Finding Order in the Disordered Hydration Shell of Rapidly Exchanging Water Molecules around the Heaviest Alkali Cs+ and Fr+. J Phys Chem B 2018; 122:12067-12076. [DOI: 10.1021/acs.jpcb.8b08414] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Santanu Roy
- Chemical Science Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee 37830, United States
| | - Vyacheslav S. Bryantsev
- Chemical Science Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee 37830, United States
| |
Collapse
|
21
|
Li W. Equipartition terms in transition path ensemble: Insights from molecular dynamics simulations of alanine dipeptide. J Chem Phys 2018; 148:084105. [PMID: 29495774 DOI: 10.1063/1.5010408] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Transition path ensemble consists of reactive trajectories and possesses all the information necessary for the understanding of the mechanism and dynamics of important condensed phase processes. However, quantitative description of the properties of the transition path ensemble is far from being established. Here, with numerical calculations on a model system, the equipartition terms defined in thermal equilibrium were for the first time estimated in the transition path ensemble. It was not surprising to observe that the energy was not equally distributed among all the coordinates. However, the energies distributed on a pair of conjugated coordinates remained equal. Higher energies were observed to be distributed on several coordinates, which are highly coupled to the reaction coordinate, while the rest were almost equally distributed. In addition, the ensemble-averaged energy on each coordinate as a function of time was also quantified. These quantitative analyses on energy distributions provided new insights into the transition path ensemble.
Collapse
Affiliation(s)
- Wenjin Li
- Institute for Advanced Study, Shenzhen University, Shenzhen, China
| |
Collapse
|
22
|
Paul S, Paul TK, Taraphder S. Reaction Coordinate, Free Energy, and Rate of Intramolecular Proton Transfer in Human Carbonic Anhydrase II. J Phys Chem B 2018; 122:2851-2866. [DOI: 10.1021/acs.jpcb.7b10713] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sanjib Paul
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Tanmoy Kumar Paul
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| | - Srabani Taraphder
- Department of Chemistry, Indian Institute of Technology Kharagpur, Kharagpur 721302, India
| |
Collapse
|
23
|
Abstract
Desolvation barriers are present for solute-solvent exchange events, such as ligand binding to an enzyme active site, during protein folding, and at battery electrodes. For solution-grown crystals, desolvation at kink sites can be the rate-limiting step for growth. However, desolvation and the associated kinetic barriers are poorly understood. In this work, we use rare-event simulation techniques to investigate attachment/detachment events at kink sites of a NaCl crystal in water. We elucidate the desolvation mechanism and present an optimized reaction coordinate, which involves one solute collective variable and one solvent collective variable. The attachment/detachment pathways for Na+ and Cl- are qualitatively similar, with quantitative differences that we attribute to different ion sizes and solvent coordination. The attachment barriers primarily result from kink site desolvation, while detachment barriers largely result from breaking ion-crystal bonds. We compute ion detachment rates from kink sites and compare with results from an independent study. We anticipate that the reaction coordinate and desolvation mechanism identified in this work may be applicable to other alkali halides.
Collapse
|
24
|
Zinovjev K, Tuñón I. Reaction coordinates and transition states in enzymatic catalysis. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1329] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Kirill Zinovjev
- Departament de Química FísicaUniversitat de València Valencia Spain
| | - Iñaki Tuñón
- Departament de Química FísicaUniversitat de València Valencia Spain
| |
Collapse
|
25
|
McGibbon RT, Husic BE, Pande VS. Identification of simple reaction coordinates from complex dynamics. J Chem Phys 2017; 146:044109. [PMID: 28147508 PMCID: PMC5272828 DOI: 10.1063/1.4974306] [Citation(s) in RCA: 55] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2016] [Accepted: 01/05/2017] [Indexed: 11/14/2022] Open
Abstract
Reaction coordinates are widely used throughout chemical physics to model and understand complex chemical transformations. We introduce a definition of the natural reaction coordinate, suitable for condensed phase and biomolecular systems, as a maximally predictive one-dimensional projection. We then show that this criterion is uniquely satisfied by a dominant eigenfunction of an integral operator associated with the ensemble dynamics. We present a new sparse estimator for these eigenfunctions which can search through a large candidate pool of structural order parameters and build simple, interpretable approximations that employ only a small number of these order parameters. Example applications with a small molecule's rotational dynamics and simulations of protein conformational change and folding show that this approach can filter through statistical noise to identify simple reaction coordinates from complex dynamics.
Collapse
Affiliation(s)
- Robert T McGibbon
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Brooke E Husic
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Vijay S Pande
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
26
|
Li W, Ma A. A benchmark for reaction coordinates in the transition path ensemble. J Chem Phys 2016; 144:134104. [PMID: 27059559 DOI: 10.1063/1.4945337] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The molecular mechanism of a reaction is embedded in its transition path ensemble, the complete collection of reactive trajectories. Utilizing the information in the transition path ensemble alone, we developed a novel metric, which we termed the emergent potential energy, for distinguishing reaction coordinates from the bath modes. The emergent potential energy can be understood as the average energy cost for making a displacement of a coordinate in the transition path ensemble. Where displacing a bath mode invokes essentially no cost, it costs significantly to move the reaction coordinate. Based on some general assumptions of the behaviors of reaction and bath coordinates in the transition path ensemble, we proved theoretically with statistical mechanics that the emergent potential energy could serve as a benchmark of reaction coordinates and demonstrated its effectiveness by applying it to a prototypical system of biomolecular dynamics. Using the emergent potential energy as guidance, we developed a committor-free and intuition-independent method for identifying reaction coordinates in complex systems. We expect this method to be applicable to a wide range of reaction processes in complex biomolecular systems.
Collapse
Affiliation(s)
- Wenjin Li
- Department of Bioengineering, the University of Illinois at Chicago, 851 South Morgan Street, Chicago, Illinois 60607, USA
| | - Ao Ma
- Department of Bioengineering, the University of Illinois at Chicago, 851 South Morgan Street, Chicago, Illinois 60607, USA
| |
Collapse
|
27
|
Lupi L, Peters B, Molinero V. Pre-ordering of interfacial water in the pathway of heterogeneous ice nucleation does not lead to a two-step crystallization mechanism. J Chem Phys 2016; 145:211910. [DOI: 10.1063/1.4961652] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Affiliation(s)
- Laura Lupi
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, USA
| | - Baron Peters
- Department of Chemical Engineering and Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, USA
| | - Valeria Molinero
- Department of Chemistry, The University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112-0850, USA
| |
Collapse
|
28
|
van Erp TS, Moqadam M, Riccardi E, Lervik A. Analyzing Complex Reaction Mechanisms Using Path Sampling. J Chem Theory Comput 2016; 12:5398-5410. [DOI: 10.1021/acs.jctc.6b00642] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Titus S. van Erp
- Department of Chemistry,
Faculty of Natural Sciences and Technology, NTNU, Norwegian University of Science and Technology, 7941 Trondheim, Norway
| | - Mahmoud Moqadam
- Department of Chemistry,
Faculty of Natural Sciences and Technology, NTNU, Norwegian University of Science and Technology, 7941 Trondheim, Norway
| | - Enrico Riccardi
- Department of Chemistry,
Faculty of Natural Sciences and Technology, NTNU, Norwegian University of Science and Technology, 7941 Trondheim, Norway
| | - Anders Lervik
- Department of Chemistry,
Faculty of Natural Sciences and Technology, NTNU, Norwegian University of Science and Technology, 7941 Trondheim, Norway
| |
Collapse
|
29
|
Banushkina PV, Krivov SV. Optimal reaction coordinates. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2016. [DOI: 10.1002/wcms.1276] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Polina V. Banushkina
- Astbury Center for Structural Molecular Biology; Faculty of Biological Sciences, University of Leeds; Leeds UK
| | - Sergei V. Krivov
- Astbury Center for Structural Molecular Biology; Faculty of Biological Sciences, University of Leeds; Leeds UK
| |
Collapse
|
30
|
Affiliation(s)
- M. C. Sherman
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - S. A. Corcelli
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| |
Collapse
|
31
|
Affiliation(s)
- Baron Peters
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106;
| |
Collapse
|
32
|
Pluhařová E, Baer MD, Schenter GK, Jungwirth P, Mundy CJ. Dependence of the Rate of LiF Ion-Pairing on the Description of Molecular Interaction. J Phys Chem B 2015; 120:1749-58. [DOI: 10.1021/acs.jpcb.5b09344] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Eva Pluhařová
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nam. 2, 16610 Prague 6, Czech Republic
| | - Marcel D. Baer
- Physical
Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Gregory K. Schenter
- Physical
Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Pavel Jungwirth
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nam. 2, 16610 Prague 6, Czech Republic
| | - Christopher J. Mundy
- Physical
Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| |
Collapse
|
33
|
Li W, Ma A. Reducing the cost of evaluating the committor by a fitting procedure. J Chem Phys 2015; 143:174103. [PMID: 26547154 PMCID: PMC4636499 DOI: 10.1063/1.4934782] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2015] [Accepted: 10/14/2015] [Indexed: 11/14/2022] Open
Abstract
Correct identification of reaction coordinates in complex systems is essential for understanding the mechanisms of their reaction dynamics. Existing methods for identifying reaction coordinates typically require knowledge of the committor--the probability of a given configuration to reach the product basin. The high computational cost of evaluating committors has limited applications of methods for identifying reaction coordinates. We proposed a fitting procedure that can reduce the cost of evaluating committors by an order of magnitude or more. The method only requires evaluating the committors of a few configurations in a transition path by the standard and costly shooting procedure. The committors of the other configurations are then estimated with great accuracy by a sigmoid function derived from fitting the few numerically evaluated committors. The method has been systematically tested on a model system of a Brownian particle moving in a one-dimensional double-well potential, and a small biomolecular system--the isomerization of alanine dipeptide in vacuum and in explicit water.
Collapse
Affiliation(s)
- Wenjin Li
- Department of Bioengineering, The University of Illinois at Chicago, 851 South Morgan St., Chicago, Illinois 60607, USA
| | - Ao Ma
- Department of Bioengineering, The University of Illinois at Chicago, 851 South Morgan St., Chicago, Illinois 60607, USA
| |
Collapse
|
34
|
Zinovjev K, Tuñón I. Transition state ensemble optimization for reactions of arbitrary complexity. J Chem Phys 2015; 143:134111. [DOI: 10.1063/1.4931596] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Affiliation(s)
- Kirill Zinovjev
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| | - Iñaki Tuñón
- Departament de Química Física, Universitat de València, 46100 Burjassot, Spain
| |
Collapse
|
35
|
Mullen RG, Shea JE, Peters B. Easy Transition Path Sampling Methods: Flexible-Length Aimless Shooting and Permutation Shooting. J Chem Theory Comput 2015; 11:2421-8. [DOI: 10.1021/acs.jctc.5b00032] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ryan Gotchy Mullen
- Department of Chemical Engineering, ‡Department of Chemistry & Biochemistry, §Department of Physics, University of California, Santa Barbara, California 93106, United States
| | - Joan-Emma Shea
- Department of Chemical Engineering, ‡Department of Chemistry & Biochemistry, §Department of Physics, University of California, Santa Barbara, California 93106, United States
| | - Baron Peters
- Department of Chemical Engineering, ‡Department of Chemistry & Biochemistry, §Department of Physics, University of California, Santa Barbara, California 93106, United States
| |
Collapse
|
36
|
Affiliation(s)
- Baron Peters
- Department
of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
- Department
of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States
| |
Collapse
|
37
|
Rosta E, Hummer G. Free Energies from Dynamic Weighted Histogram Analysis Using Unbiased Markov State Model. J Chem Theory Comput 2014; 11:276-85. [DOI: 10.1021/ct500719p] [Citation(s) in RCA: 96] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Edina Rosta
- Department
of Chemistry, King’s College London, London SE1 1DB, United Kingdom
| | - Gerhard Hummer
- Department
of Theoretical Biophysics, Max Planck Institute of Biophysics, 60438 Frankfurt am Main, Germany
| |
Collapse
|
38
|
Zinovjev K, Tuñón I. Exploring chemical reactivity of complex systems with path-based coordinates: Role of the distance metric. J Comput Chem 2014; 35:1672-81. [DOI: 10.1002/jcc.23673] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2014] [Revised: 05/16/2014] [Accepted: 06/18/2014] [Indexed: 01/06/2023]
Affiliation(s)
- Kirill Zinovjev
- Departament de Química Física; Universitat de València; 46100 Burjassot Spain
| | - Iñaki Tuñón
- Departament de Química Física; Universitat de València; 46100 Burjassot Spain
| |
Collapse
|
39
|
Annapureddy HVR, Dang LX. Understanding the rates and molecular mechanism of water-exchange around aqueous ions using molecular simulations. J Phys Chem B 2014; 118:8917-27. [PMID: 24911526 DOI: 10.1021/jp502922c] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Solvation processes occurring around aqueous ions are of fundamental importance in physics, chemistry, and biology. Over the past few decades, several experimental and theoretical studies were devoted to understanding ion solvation and the processes involved in it. In this article, we present a summary of our recent efforts that, through computer simulations, focused on providing a comprehensive understanding of solvent-exchange processes around aqueous ions. To accomplish these activities, we have looked at the mechanistic properties associated with the water-exchange process, such as potentials of mean force, time-dependent transmission coefficients, and the corresponding rate constants using transition state theory, the reactive flux method, and Grote-Hynes treatments of the dynamic response of the solvent.
Collapse
Affiliation(s)
- Harsha V R Annapureddy
- Physical Sciences Division, Pacific Northwest National Laboratory , Richland, Washington 93352, United States
| | | |
Collapse
|
40
|
Agarwal V, Peters B. Solute Precipitate Nucleation: A Review of Theory and Simulation Advances. ADVANCES IN CHEMICAL PHYSICS 2014. [DOI: 10.1002/9781118755815.ch03] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
41
|
Daru J, Stirling A. Divided Saddle Theory: A New Idea for Rate Constant Calculation. J Chem Theory Comput 2014; 10:1121-7. [DOI: 10.1021/ct400970y] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- János Daru
- Institute
of Organic Chemistry, Research Centre for Natural Sciences of the Hungarian Academy of Sciences, Budapest 1025, Hungary
- Institute
of Chemistry, Loránd Eötvös University, Budapest 1518, Hungary
| | - András Stirling
- Institute
of Organic Chemistry, Research Centre for Natural Sciences of the Hungarian Academy of Sciences, Budapest 1025, Hungary
| |
Collapse
|
42
|
Mullen RG, Shea JE, Peters B. Transmission Coefficients, Committors, and Solvent Coordinates in Ion-Pair Dissociation. J Chem Theory Comput 2014; 10:659-67. [DOI: 10.1021/ct4009798] [Citation(s) in RCA: 89] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ryan Gotchy Mullen
- Department of Chemical Engineering, ‡Department of Chemistry & Biochemistry, §Department of Physics, University of California, Santa Barbara, California 93106
| | - Joan-Emma Shea
- Department of Chemical Engineering, ‡Department of Chemistry & Biochemistry, §Department of Physics, University of California, Santa Barbara, California 93106
| | - Baron Peters
- Department of Chemical Engineering, ‡Department of Chemistry & Biochemistry, §Department of Physics, University of California, Santa Barbara, California 93106
| |
Collapse
|
43
|
Mullen RG, Shea JE, Peters B. Communication: An existence test for dividing surfaces without recrossing. J Chem Phys 2014; 140:041104. [DOI: 10.1063/1.4862504] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
44
|
Li W, Ma A. Recent developments in methods for identifying reaction coordinates. MOLECULAR SIMULATION 2014; 40:784-793. [PMID: 25197161 PMCID: PMC4152980 DOI: 10.1080/08927022.2014.907898] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
In the study of rare events in complex systems with many degrees of freedom, a key element is to identify the reaction coordinates of a given process. Over recent years, a number of methods and protocols have been developed to extract the reaction coordinates based on limited information from molecular dynamics simulations. In this review, we provide a brief survey over a number of major methods developed in the past decade, some of which are discussed in greater detail, to provide an overview of the problems that are partially solved and challenges that still remain. A particular emphasis has been placed on methods for identifying reaction coordinates that are related to the committor.
Collapse
Affiliation(s)
- Wenjin Li
- Department of Bioengineering, The University of Illinois at Chicago, 851 South Morgan Street, Chicago, IL 60607, USA
| | - Ao Ma
- Department of Bioengineering, The University of Illinois at Chicago, 851 South Morgan Street, Chicago, IL 60607, USA
| |
Collapse
|