1
|
Gorantla K, Krishnan A, Waheed SO, Varghese A, DiCastri I, LaRouche C, Paik M, Fields GB, Karabencheva-Christova TG. Novel Insights into the Catalytic Mechanism of Collagenolysis by Zn(II)-Dependent Matrix Metalloproteinase-1. Biochemistry 2024; 63:1925-1940. [PMID: 38963231 PMCID: PMC11309001 DOI: 10.1021/acs.biochem.4c00076] [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/13/2024] [Revised: 06/14/2024] [Accepted: 06/24/2024] [Indexed: 07/05/2024]
Abstract
Collagen hydrolysis, catalyzed by Zn(II)-dependent matrix metalloproteinases (MMPs), is a critical physiological process. Despite previous computational investigations into the catalytic mechanisms of MMP-mediated collagenolysis, a significant knowledge gap in understanding remains regarding the influence of conformational sampling and entropic contributions at physiological temperature on enzymatic collagenolysis. In our comprehensive multilevel computational study, employing quantum mechanics/molecular mechanics (QM/MM) metadynamics (MetD) simulations, we aimed to bridge this gap and provide valuable insights into the catalytic mechanism of MMP-1. Specifically, we compared the full enzyme-substrate complex in solution, clusters in solution, and gas-phase to elucidate insights into MMP-1-catalyzed collagenolysis. Our findings reveal significant differences in the catalytic mechanism when considering thermal effects and the dynamic evolution of the system, contrasting with conventional static potential energy surface QM/MM reaction path studies. Notably, we observed a significant stabilization of the critical tetrahedral intermediate, attributed to contributions from conformational flexibility and entropy. Moreover, we found that protonation of the scissile bond nitrogen occurs via proton transfer from a Zn(II)-coordinated hydroxide rather than from a solvent water molecule. Following C-N bond cleavage, the C-terminus remains coordinated to the catalytic Zn(II), while the N-terminus forms a hydrogen bond with a solvent water molecule. Subsequently, the release of the C-terminus is facilitated by the coordination of a water molecule. Our study underscores the pivotal role of protein conformational dynamics at physiological temperature in stabilizing the transition state of the rate-limiting step and key intermediates, compared to the corresponding reaction in solution. These fundamental insights into the mechanism of collagen degradation provide valuable guidance for the development of MMP-1-specific inhibitors.
Collapse
Affiliation(s)
- Koteswara
Rao Gorantla
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Anandhu Krishnan
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Sodiq O. Waheed
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Ann Varghese
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Isabella DiCastri
- Department
of Chemical Engineering, Michigan Technological
University, Houghton, Michigan 49931, United States
| | - Ciara LaRouche
- Department
of Chemical Engineering, Michigan Technological
University, Houghton, Michigan 49931, United States
| | - Meredith Paik
- Department
of Chemistry, Michigan Technological University, Houghton, Michigan 49931, United States
| | - Gregg B. Fields
- Department
of Chemistry and Biochemistry and I-HEALTH, Florida Atlantic University, Jupiter, Florida 33458, United States
| | | |
Collapse
|
2
|
Yuan S, Han X, Zhang J, Xie Z, Fan C, Xiao Y, Gao YQ, Yang YI. Generating High-Precision Force Fields for Molecular Dynamics Simulations to Study Chemical Reaction Mechanisms Using Molecular Configuration Transformer. J Phys Chem A 2024; 128:4378-4390. [PMID: 38759697 DOI: 10.1021/acs.jpca.4c01267] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/19/2024]
Abstract
Theoretical studies on chemical reaction mechanisms have been crucial in organic chemistry. Traditionally, calculating the manually constructed molecular conformations of transition states for chemical reactions using quantum chemical calculations is the most commonly used method. However, this way is heavily dependent on individual experience and chemical intuition. In our previous study, we proposed a research paradigm that used enhanced sampling in molecular dynamics simulations to study chemical reactions. This approach can directly simulate the entire process of a chemical reaction. However, the computational speed limited the use of high-precision potential energy functions for simulations. To address this issue, we presented a scheme for training high-precision force fields for molecular modeling using a previously developed graph-neural-network-based molecular model, molecular configuration transformer. This potential energy function allowed for highly accurate simulations at a low computational cost, leading to more precise calculations of the mechanism of chemical reactions. We applied this approach to study a Claisen rearrangement reaction and a carbonyl insertion reaction catalyzed by manganese.
Collapse
Affiliation(s)
- Sihao Yuan
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518132, China
| | - Xu Han
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Jun Zhang
- Changping Laboratory, Beijing 102200, China
| | - Zhaoxin Xie
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518132, China
| | - Cheng Fan
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518132, China
| | - Yunlong Xiao
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Yi Qin Gao
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518132, China
- Changping Laboratory, Beijing 102200, China
| | - Yi Isaac Yang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518132, China
| |
Collapse
|
3
|
Han X, Zhang J, Yang YI, Zhang Z, Yang L, Gao YQ. Enhanced Sampling Simulation Reveals How Solvent Influences Chirogenesis of the Intra-Molecular Diels-Alder Reaction. J Chem Theory Comput 2022; 18:4318-4326. [PMID: 35666128 DOI: 10.1021/acs.jctc.2c00233] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
The timescale involved in chemical reactions is quite often beyond that of normal molecular dynamics simulations. Here, we combine metadynamics with selective integrated tempering sampling to simulate an intra-molecular Diels-Alder reaction in explicit solvents. Based on a one-dimensional collective variable obtained from harmonic linear discriminant analysis, four chiral isomers of products were observed in the simulation. Analyses of reactive trajectories showed that this reaction follows a concerted mechanism in all four solvents. In addition, the hydrogen bond between the reactant and water solvent plays an important role in the water-accelerated reaction mechanism. The dynamics of chirality formation varies significantly with solvents. The chirality of products forms significantly before the transition state, especially in ionic liquid.
Collapse
Affiliation(s)
- Xu Han
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Jun Zhang
- Changping Laboratory, Beijing 102206, China
| | - Yi Isaac Yang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen, Guangdong 518132, China
| | - Zhen Zhang
- School of Physics and Technology, Tangshan Normal University, Tangshan 063000, China
| | - Lijiang Yang
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Yi Qin Gao
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China.,Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen, Guangdong 518132, China
| |
Collapse
|
4
|
Huang Y, Xia Y, Yang L, Wei J, Yang YI, Gao YQ. SPONGE
: A
GPU‐Accelerated
Molecular Dynamics Package with Enhanced Sampling and
AI‐Driven
Algorithms. CHINESE J CHEM 2022. [DOI: 10.1002/cjoc.202100456] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Yu‐Peng Huang
- College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
- Beijing National Laboratory for Molecular Sciences Peking University Beijing 100871 China
- Biomedical Pioneering Innovation Center Peking University Beijing 100871 China
| | - Yijie Xia
- College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
- Beijing National Laboratory for Molecular Sciences Peking University Beijing 100871 China
- Biomedical Pioneering Innovation Center Peking University Beijing 100871 China
| | - Lijiang Yang
- College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
- Beijing National Laboratory for Molecular Sciences Peking University Beijing 100871 China
- Biomedical Pioneering Innovation Center Peking University Beijing 100871 China
- Beijing Advanced Innovation Center for Genomics Peking University Beijing 100871 China
| | - Jiachen Wei
- State Key Laboratory of Nonlinear Mechanics and Beijing Key Laboratory of Engineered Construction and Mechanobiology, Institute of Mechanics Chinese Academy of Sciences Beijing 100190 China
- Shenzhen Bay Laboratory, Gaoke Innovation Center, Guangqiao Road, Guangming District Shenzhen Guangdong 518132 China
| | - Yi Isaac Yang
- Shenzhen Bay Laboratory, Gaoke Innovation Center, Guangqiao Road, Guangming District Shenzhen Guangdong 518132 China
| | - Yi Qin Gao
- College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
- Beijing National Laboratory for Molecular Sciences Peking University Beijing 100871 China
- Biomedical Pioneering Innovation Center Peking University Beijing 100871 China
- Beijing Advanced Innovation Center for Genomics Peking University Beijing 100871 China
- Shenzhen Bay Laboratory, Gaoke Innovation Center, Guangqiao Road, Guangming District Shenzhen Guangdong 518132 China
| |
Collapse
|
5
|
Kitanosono T, Kobayashi S. Synthetic Organic "Aquachemistry" that Relies on Neither Cosolvents nor Surfactants. ACS CENTRAL SCIENCE 2021; 7:739-747. [PMID: 34079894 PMCID: PMC8161484 DOI: 10.1021/acscentsci.1c00045] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/12/2021] [Indexed: 06/12/2023]
Abstract
There is a growing awareness of the underlying power of catalytic reactions in water that is not limited to innate sustainability alone. Some Type III reactions are catalytically accelerated without dissolution of reactants and are occasionally highly selective, as shown by comparison with the corresponding reactions run in organic solvents or under solvent-free conditions. Such catalysts are highly diversified, including hydrophilic, lipophilic, and even solid catalysts. In this Outlook, we highlight the impressive characteristics of illustrative catalysis that is exerted despite the immiscibility of the substrates and reveal the intrinsic benefits of these enigmatic reactions for synthetic organic chemistry, albeit with many details remaining unclear. We hope that this brief introduction to the expanding field of synthetic organic "aquachemistry" will inspire organic chemists to use the platform to invent new transformations.
Collapse
Affiliation(s)
- Taku Kitanosono
- Department of Chemistry,
School of Science, The University of Tokyo, Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
| | - Shu Kobayashi
- Department of Chemistry,
School of Science, The University of Tokyo, Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
| |
Collapse
|
6
|
Zhang J, Lei YK, Zhang Z, Han X, Li M, Yang L, Yang YI, Gao YQ. Deep reinforcement learning of transition states. Phys Chem Chem Phys 2021; 23:6888-6895. [PMID: 33729229 DOI: 10.1039/d0cp06184k] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Combining reinforcement learning (RL) and molecular dynamics (MD) simulations, we propose a machine-learning approach, called RL‡, to automatically unravel chemical reaction mechanisms. In RL‡, locating the transition state of a chemical reaction is formulated as a game, and two functions are optimized, one for value estimation and the other for policy making, to iteratively improve our chance of winning this game. Both functions can be approximated by deep neural networks. By virtue of RL‡, one can directly interpret the reaction mechanism according to the value function. Meanwhile, the policy function allows efficient sampling of the transition path ensemble, which can be further used to analyze reaction dynamics and kinetics. Through multiple experiments, we show that RL‡ can be trained tabula rasa hence allowing us to reveal chemical reaction mechanisms with minimal subjective biases.
Collapse
Affiliation(s)
- Jun Zhang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, 518055 Shenzhen, China.
| | | | | | | | | | | | | | | |
Collapse
|
7
|
Rahnamoun A, Kaymak MC, Manathunga M, Götz AW, van Duin ACT, Merz KM, Aktulga HM. ReaxFF/AMBER-A Framework for Hybrid Reactive/Nonreactive Force Field Molecular Dynamics Simulations. J Chem Theory Comput 2020; 16:7645-7654. [PMID: 33141581 DOI: 10.1021/acs.jctc.0c00874] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Combined quantum mechanical/molecular mechanical (QM/MM) models using semiempirical and ab initio methods have been extensively reported on over the past few decades. These methods have been shown to be capable of providing unique insights into a range of problems, but they are still limited to relatively short time scales, especially QM/MM models using ab initio methods. An intermediate approach between a QM based model and classical mechanics could help fill this time-scale gap and facilitate the study of a range of interesting problems. Reactive force fields represent the intermediate approach explored in this paper. A widely used reactive model is ReaxFF, which has largely been applied to materials science problems and is generally used as a stand-alone (i.e., the full system is modeled using ReaxFF). We report a hybrid ReaxFF/AMBER molecular dynamics (MD) tool, which introduces ReaxFF capabilities to capture bond breaking and formation within the AMBER MD software package. This tool enables us to study local reactive events in large systems at a fraction of the computational costs of QM/MM models. We describe the implementation of ReaxFF/AMBER, validate this implementation using a benzene molecule solvated in water, and compare its performance against a range of similar approaches. To illustrate the predictive capabilities of ReaxFF/AMBER, we carried out a Claisen rearrangement study in aqueous solution. In a first for ReaxFF, we were able to use AMBER's potential of mean force (PMF) capabilities to perform a PMF study on this organic reaction. The ability to capture local reaction events in large systems using combined ReaxFF/AMBER opens up a range of problems that can be tackled using this model to address both chemical and biological processes.
Collapse
Affiliation(s)
- Ali Rahnamoun
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, 578 S. Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Mehmet Cagri Kaymak
- Department of Computer Science and Engineering, Michigan State University, 428 S. Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Madushanka Manathunga
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, 578 S. Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Andreas W Götz
- San Diego Supercomputer Center, University of California San Diego, 9500 Gilman Drive, La Jolla, California 92093-0505, United States
| | - Adri C T van Duin
- Department of Mechanical Engineering, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Kenneth M Merz
- Department of Chemistry and Department of Biochemistry and Molecular Biology, Michigan State University, 578 S. Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Hasan Metin Aktulga
- Department of Computer Science and Engineering, Michigan State University, 428 S. Shaw Lane, East Lansing, Michigan 48824-1322, United States
| |
Collapse
|
8
|
Kitanosono T, Kobayashi S. Reactions in Water Involving the “On‐Water” Mechanism. Chemistry 2020; 26:9408-9429. [DOI: 10.1002/chem.201905482] [Citation(s) in RCA: 74] [Impact Index Per Article: 18.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2019] [Revised: 02/08/2020] [Indexed: 11/07/2022]
Affiliation(s)
- Taku Kitanosono
- Department of ChemistrySchool of ScienceThe University of Tokyo Hongo Bunkyo-ku Tokyo 113-0033 Japan
| | - Shū Kobayashi
- Department of ChemistrySchool of ScienceThe University of Tokyo Hongo Bunkyo-ku Tokyo 113-0033 Japan
| |
Collapse
|
9
|
Yang YI, Shao Q, Zhang J, Yang L, Gao YQ. Enhanced sampling in molecular dynamics. J Chem Phys 2019; 151:070902. [PMID: 31438687 DOI: 10.1063/1.5109531] [Citation(s) in RCA: 177] [Impact Index Per Article: 35.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Abstract
Although molecular dynamics simulations have become a useful tool in essentially all fields of chemistry, condensed matter physics, materials science, and biology, there is still a large gap between the time scale which can be reached in molecular dynamics simulations and that observed in experiments. To address the problem, many enhanced sampling methods were introduced, which effectively extend the time scale being approached in simulations. In this perspective, we review a variety of enhanced sampling methods. We first discuss collective-variables-based methods including metadynamics and variationally enhanced sampling. Then, collective variable free methods such as parallel tempering and integrated tempering methods are presented. At last, we conclude with a brief introduction of some newly developed combinatory methods. We summarize in this perspective not only the theoretical background and numerical implementation of these methods but also the new challenges and prospects in the field of the enhanced sampling.
Collapse
Affiliation(s)
- Yi Isaac Yang
- Institute of Systems Biology, Shenzhen Bay Laboratory, Shenzhen 518055, Guangdong, China
| | - Qiang Shao
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
| | - Jun Zhang
- Department of Mathematics and Computer Science, Freie Universität Berlin, Arnimallee 6, Berlin 14195, Germany
| | - Lijiang Yang
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Yi Qin Gao
- Institute of Systems Biology, Shenzhen Bay Laboratory, Shenzhen 518055, Guangdong, China
| |
Collapse
|
10
|
Yang X, Zou J, Wang Y, Xue Y, Yang S. Role of Water in the Reaction Mechanism and endo/exo Selectivity of 1,3-Dipolar Cycloadditions Elucidated by Quantum Chemistry and Machine Learning. Chemistry 2019; 25:8289-8303. [PMID: 30887586 DOI: 10.1002/chem.201900617] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2019] [Indexed: 02/05/2023]
Abstract
Asymmetric 1,3-dipolar cycloadditions of azomethine ylides with activated olefins are among the most important and versatile methods for the synthesis of enantioenriched pyrroline and pyrrolidine derivatives. Despite both theoretical and practical importance, the role of water molecules in the reactivity and endo/exo selectivity remains unclear. To explore how water accelerates the reactions and improves the endo/exo selectivity of the cycloadditions of 1,3-dipole phthalazinium-2-dicyanomethanide (1) and two dipolarophiles, an ab initio-quality neural network potential that overcomes the computational bottleneck of explicitly considering water molecules was used. It is demonstrated that not only the nature of both the dipolarophile and the 1,3-dipole, but also the solvent medium, can perturb or even alter the reaction mechanism. An extreme case was found for the reaction of 1,3-dipole 1 with methyl vinyl ketone, in which the reaction mechanism changes from a concerted to a stepwise mode on going from MeCN to H2 O as solvent, with formation of a zwitterionic intermediate that is a very shallow minimum on the energy surface. Thus, high stereocontrol can still be expected despite the stepwise nature of the mechanism. The results indicate that water can induce global polarization along the reaction coordinate and highlight the role of microsolvation effects and bulk-phase effects in reproducing the experimentally observed aqueous acceleration and enhanced endo/exo selectivity.
Collapse
Affiliation(s)
- Xin Yang
- State Key Laboratory of Biotherapy and Cancer Center, West China Hospital, Sichuan University, Chengdu, Sichuan, 610041, P.R. China
| | - Jun Zou
- State Key Laboratory of Biotherapy and Cancer Center, West China Hospital, Sichuan University, Chengdu, Sichuan, 610041, P.R. China
| | - Yifei Wang
- State Key Laboratory of Biotherapy and Cancer Center, West China Hospital, Sichuan University, Chengdu, Sichuan, 610041, P.R. China
| | - Ying Xue
- College of Chemistry, Key Lab of Green Chemistry and Technology in Ministry of Education, Sichuan University, Chengdu, Sichuan, 610041, P.R. China
| | - Shengyong Yang
- State Key Laboratory of Biotherapy and Cancer Center, West China Hospital, Sichuan University, Chengdu, Sichuan, 610041, P.R. China
| |
Collapse
|
11
|
Xie L, Cheng H, Fang D, Chen ZN, Yang M. Enhanced QM/MM sampling for free energy calculation of chemical reactions: A case study of double proton transfer. J Chem Phys 2019; 150:044111. [PMID: 30709281 DOI: 10.1063/1.5072779] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Free energy calculations for chemical reactions with a steep energy barrier require well defined reaction coordinates (RCs). However, when multiple parallel channels exist along selected RC, the application of conventional enhanced samplings is difficult to generate correct sampling within limited simulation time and thus cannot give correct prediction about the favorable pathways, the relative stability of multiple products or intermediates. Here, we implement the selective integrated tempering sampling (SITS) method with quantum mechanical and molecular mechanical (QM/MM) potential to investigate the chemical reactions in solution. The combined SITS-QM/MM scheme is used to identify possible reaction paths, intermediate and product states, and the free energy profiles for the different reaction paths. Two double proton transfer reactions were studied to validate the implemented method and simulation protocol, from which the independent and correlated proton transfer processes are identified in two representative systems, respectively. This protocol can be generalized to various kinds of chemical reactions for both academic studies and industry applications, such as in exploration and optimization of potential reactions in DNA encoded compound library and halogen or deuterium substitution of the hit discovery and lead optimization stages of drug design via providing a better understanding of the reaction mechanism along the designed chemical reaction pathways.
Collapse
Affiliation(s)
- Liangxu Xie
- Institute of Bioinformatics and Medical Engineering, School of Electrical and Information Engineering, Jiangsu University of Technology, Changzhou 21300, China
| | - Huimin Cheng
- XtalPi Inc. (Shenzhen Jingtai Technology Co., Ltd.), Times Science & Tech Mansion E. 20F, 7028 Shennan Ave, Futian District, Shenzhen, China
| | - Dong Fang
- XtalPi Inc. (Shenzhen Jingtai Technology Co., Ltd.), Times Science & Tech Mansion E. 20F, 7028 Shennan Ave, Futian District, Shenzhen, China
| | - Zhe-Ning Chen
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, 155 Yangqiao Road West, Fuzhou 350002, China
| | - Mingjun Yang
- XtalPi Inc. (Shenzhen Jingtai Technology Co., Ltd.), Times Science & Tech Mansion E. 20F, 7028 Shennan Ave, Futian District, Shenzhen, China
| |
Collapse
|
12
|
Chang YL, Sasaki T, Ribas-Ariño J, Machida M, Shiga M. Understanding Competition of Polyalcohol Dehydration Reactions in Hot Water. J Phys Chem B 2019; 123:1662-1671. [DOI: 10.1021/acs.jpcb.8b11615] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Yong Lik Chang
- Department of Complexity Science and Engineering, Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5, Kashiwanoha, Kashiwa 277-8561, Japan
| | - Takehiko Sasaki
- Department of Complexity Science and Engineering, Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5, Kashiwanoha, Kashiwa 277-8561, Japan
| | - Jordi Ribas-Ariño
- Departament de Química-Física i CERQT, Universitat de Barcelona, Diagonal, 645, 08028 Barcelona, Spain
| | - Masahiko Machida
- Center for Computational Science and e-Systems, Japan Atomic Energy Agency, 178-4-4, Wakashiba, Kashiwa, Chiba 277-0871, Japan
| | - Motoyuki Shiga
- Center for Computational Science and e-Systems, Japan Atomic Energy Agency, 178-4-4, Wakashiba, Kashiwa, Chiba 277-0871, Japan
| |
Collapse
|
13
|
Brickel S, Meuwly M. Molecular Determinants for Rate Acceleration in the Claisen Rearrangement Reaction. J Phys Chem B 2018; 123:448-456. [DOI: 10.1021/acs.jpcb.8b11059] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Sebastian Brickel
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, Basel CH-4056, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, Basel CH-4056, Switzerland
| |
Collapse
|
14
|
Wu J, Shen L, Yang W. Internal force corrections with machine learning for quantum mechanics/molecular mechanics simulations. J Chem Phys 2018; 147:161732. [PMID: 29096448 DOI: 10.1063/1.5006882] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Ab initio quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulation is a useful tool to calculate thermodynamic properties such as potential of mean force for chemical reactions but intensely time consuming. In this paper, we developed a new method using the internal force correction for low-level semiempirical QM/MM molecular dynamics samplings with a predefined reaction coordinate. As a correction term, the internal force was predicted with a machine learning scheme, which provides a sophisticated force field, and added to the atomic forces on the reaction coordinate related atoms at each integration step. We applied this method to two reactions in aqueous solution and reproduced potentials of mean force at the ab initio QM/MM level. The saving in computational cost is about 2 orders of magnitude. The present work reveals great potentials for machine learning in QM/MM simulations to study complex chemical processes.
Collapse
Affiliation(s)
- Jingheng Wu
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Lin Shen
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| |
Collapse
|
15
|
Yan F, Liu X, Zhang S, Su J, Zhang Q, Chen J. Computational revelation of binding mechanisms of inhibitors to endocellular protein tyrosine phosphatase 1B using molecular dynamics simulations. J Biomol Struct Dyn 2017; 36:3636-3650. [DOI: 10.1080/07391102.2017.1394221] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Fangfang Yan
- School of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| | - Jing Su
- School of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, 250357, China
| |
Collapse
|
16
|
Zhang J, Zhang Z, Yang YI, Liu S, Yang L, Gao YQ. Rich Dynamics Underlying Solution Reactions Revealed by Sampling and Data Mining of Reactive Trajectories. ACS CENTRAL SCIENCE 2017; 3:407-414. [PMID: 28573202 PMCID: PMC5445542 DOI: 10.1021/acscentsci.7b00037] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/23/2017] [Indexed: 05/15/2023]
Abstract
Efficient sampling in both configuration and trajectory spaces, combined with mechanism analyses via data mining, allows a systematic investigation of the thermodynamics, kinetics, and molecular-detailed dynamics of chemical reactions in solution. Through a Bayesian learning algorithm, the reaction coordinate(s) of a (retro-)Claisen rearrangement in bulk water was variationally optimized. The bond formation/breakage was found to couple with intramolecular charge separation and dipole change, and significant dynamic solvent effects manifest, leading to the "in-water" acceleration of Claisen rearrangement. In addition, the vibrational modes of the reactant and the solvation states are significantly coupled to the reaction dynamics, leading to heterogeneous and oscillatory reaction paths. The calculated reaction rate is well interpreted by the Kramers' theory with a diffusion term accounting for solvent-solute interactions. These findings demonstrated that the reaction mechanisms can be complicated in homogeneous solutions since the solvent-solute interactions can profoundly influence the reaction dynamics and the energy transfer process.
Collapse
Affiliation(s)
- Jun Zhang
- Institute of Theoretical and Computational Chemistry, College of
Chemistry and Molecular Engineering and Biodynamic Optical Imaging Center, Peking University, Beijing 100871, China
| | - Zhen Zhang
- Institute of Theoretical and Computational Chemistry, College of
Chemistry and Molecular Engineering and Biodynamic Optical Imaging Center, Peking University, Beijing 100871, China
| | - Yi Isaac Yang
- Institute of Theoretical and Computational Chemistry, College of
Chemistry and Molecular Engineering and Biodynamic Optical Imaging Center, Peking University, Beijing 100871, China
| | - Sirui Liu
- Institute of Theoretical and Computational Chemistry, College of
Chemistry and Molecular Engineering and Biodynamic Optical Imaging Center, Peking University, Beijing 100871, China
| | - Lijiang Yang
- Institute of Theoretical and Computational Chemistry, College of
Chemistry and Molecular Engineering and Biodynamic Optical Imaging Center, Peking University, Beijing 100871, China
| | - Yi Qin Gao
- Institute of Theoretical and Computational Chemistry, College of
Chemistry and Molecular Engineering and Biodynamic Optical Imaging Center, Peking University, Beijing 100871, China
| |
Collapse
|
17
|
Shen L, Xie L, Yang M. Thermodynamic properties of solvated peptides from selective integrated tempering sampling with a new weighting factor estimation algorithm. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1292009] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Lin Shen
- Department of Chemistry, The University of Hong Kong, Hong Kong, P. R. China
| | - Liangxu Xie
- Department of Chemistry, The University of Hong Kong, Hong Kong, P. R. China
| | - Mingjun Yang
- Department of Chemistry, The University of Hong Kong, Hong Kong, P. R. China
| |
Collapse
|
18
|
Shen L, Wu J, Yang W. Multiscale Quantum Mechanics/Molecular Mechanics Simulations with Neural Networks. J Chem Theory Comput 2016; 12:4934-4946. [PMID: 27552235 PMCID: PMC6209101 DOI: 10.1021/acs.jctc.6b00663] [Citation(s) in RCA: 73] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Molecular dynamics simulation with multiscale quantum mechanics/molecular mechanics (QM/MM) methods is a very powerful tool for understanding the mechanism of chemical and biological processes in solution or enzymes. However, its computational cost can be too high for many biochemical systems because of the large number of ab initio QM calculations. Semiempirical QM/MM simulations have much higher efficiency. Its accuracy can be improved with a correction to reach the ab initio QM/MM level. The computational cost on the ab initio calculation for the correction determines the efficiency. In this paper we developed a neural network method for QM/MM calculation as an extension of the neural-network representation reported by Behler and Parrinello. With this approach, the potential energy of any configuration along the reaction path for a given QM/MM system can be predicted at the ab initio QM/MM level based on the semiempirical QM/MM simulations. We further applied this method to three reactions in water to calculate the free energy changes. The free-energy profile obtained from the semiempirical QM/MM simulation is corrected to the ab initio QM/MM level with the potential energies predicted with the constructed neural network. The results are in excellent accordance with the reference data that are obtained from the ab initio QM/MM molecular dynamics simulation or corrected with direct ab initio QM/MM potential energies. Compared with the correction using direct ab initio QM/MM potential energies, our method shows a speed-up of 1 or 2 orders of magnitude. It demonstrates that the neural network method combined with the semiempirical QM/MM calculation can be an efficient and reliable strategy for chemical reaction simulations.
Collapse
Affiliation(s)
- Lin Shen
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Jingheng Wu
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
- School of Chemistry and Chemical Engineering, Sun Yat-sen University, Guangzhou 510275, P.R. China
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
19
|
Yang L, Zhang J, Che X, Gao YQ. Simulation Studies of Protein and Small Molecule Interactions and Reaction. Methods Enzymol 2016; 578:169-212. [PMID: 27497167 DOI: 10.1016/bs.mie.2016.05.031] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Computational studies of protein and small molecule (protein-ligand/enzyme-substrate) interactions become more and more important in biological science and drug discovery. Computer modeling can provide molecular details of the processes such as conformational change, binding, and transportation of small molecules/proteins, which are not easily to be captured in experiments. In this chapter, we discussed simulation studies of both protein and small molecules from three aspects: conformation sampling, transportations of small molecules in enzymes, and enzymatic reactions involving small molecules. Both methodology developments and examples of simulation studies in this field were presented.
Collapse
Affiliation(s)
- L Yang
- College of Chemistry and Molecular Engineering, Beijing National Laboratory for Molecular Sciences, Peking University, Beijing, PR China; Biodynamic Optical Imaging Center (BIOPIC), Peking University, Beijing, PR China
| | - J Zhang
- College of Chemistry and Molecular Engineering, Beijing National Laboratory for Molecular Sciences, Peking University, Beijing, PR China; Biodynamic Optical Imaging Center (BIOPIC), Peking University, Beijing, PR China
| | - X Che
- College of Chemistry and Molecular Engineering, Beijing National Laboratory for Molecular Sciences, Peking University, Beijing, PR China; Biodynamic Optical Imaging Center (BIOPIC), Peking University, Beijing, PR China
| | - Y Q Gao
- College of Chemistry and Molecular Engineering, Beijing National Laboratory for Molecular Sciences, Peking University, Beijing, PR China; Biodynamic Optical Imaging Center (BIOPIC), Peking University, Beijing, PR China.
| |
Collapse
|
20
|
Ruiz-Barragan S, Ribas Ariño J, Shiga M. The reaction mechanism of polyalcohol dehydration in hot pressurized water. Phys Chem Chem Phys 2016; 18:32438-32447. [DOI: 10.1039/c6cp05695d] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
The use of high-temperature liquid water (HTW) as a reaction medium is a very promising technology in the field of green chemistry.
Collapse
Affiliation(s)
- Sergi Ruiz-Barragan
- CCSE
- Japan Atomic Energy Agency
- Kashiwa
- Japan
- Department of Theoretical and Computational Molecular Science
| | - Jordi Ribas Ariño
- Departament de Química-Física i CERQT
- Universitat de Barcelona
- 08028-Barcelona
- Spain
| | | |
Collapse
|