1
|
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
|
2
|
Kitanosono T, Hashidoko A, Yamashita Y, Kobayashi S. 2-Methoxyethyl Nitrite as a Reagent for Chemoselective On-Water Nitration. Chem Asian J 2022; 17:e202200457. [PMID: 35612572 DOI: 10.1002/asia.202200457] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2022] [Revised: 05/23/2022] [Indexed: 11/11/2022]
Abstract
An on-water approach has been developed that allows a nitration of tyrosines and phenols under mild conditions. We envisioned that the assembly of tyrosine/tyrosyl radical couples with interfacial water molecules would realize a biomimetic stacking hydrogen atom transfer (HAT) transition state to facilitate the electron-transfer process. The optimal organic nitrite, 2-methoxyethyl nitrite, resulted in rapid coupling of the tyrosyl radicals with •NO 2 at the oil-water interface to afford the nitrated phenols. Many characteristics found in our on-water strategy are distinct from other complementary systems that include radical nitration. These enticing roles of water in the reaction process introduce new avenues to explore in the design of synthetic organic chemistry systems.
Collapse
Affiliation(s)
- Taku Kitanosono
- The University of Tokyo: Tokyo Daigaku, Department of Chemistry, JAPAN
| | - Airu Hashidoko
- The University of Tokyo: Tokyo Daigaku, Department of Chemistry, JAPAN
| | | | - Shu Kobayashi
- The University of Tokyo, Department of Chemistry, School of Science, 7-3-1 Hongo, Bunkyo-ku, 113-0033, Tokyo, JAPAN
| |
Collapse
|
3
|
BF3–Catalyzed Diels–Alder Reaction between Butadiene and Methyl Acrylate in Aqueous Solution—An URVA and Local Vibrational Mode Study. Catalysts 2022. [DOI: 10.3390/catal12040415] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
In this study we investigate the Diels–Alder reaction between methyl acrylate and butadiene, which is catalyzed by BF3 Lewis acid in explicit water solution, using URVA and Local Mode Analysis as major tools complemented with NBO, electron density and ring puckering analyses. We considered four different starting orientations of methyl acrylate and butadiene, which led to 16 DA reactions in total. In order to isolate the catalytic effects of the BF3 catalyst and those of the water environment and exploring how these effects are synchronized, we systematically compared the non-catalyzed reaction in gas phase and aqueous solution with the catalyzed reaction in gas phase and aqueous solution. Gas phase studies were performed at the B3LYP/6-311+G(2d,p) level of theory and studies in aqueous solution were performed utilizing a QM/MM approach at the B3LYP/6-311+G(2d,p)/AMBER level of theory. The URVA results revealed reaction path curvature profiles with an overall similar pattern for all 16 reactions showing the same sequence of CC single bond formation for all of them. In contrast to the parent DA reaction with symmetric substrates causing a synchronous bond formation process, here, first the new CC single bond on the CH2 side of methyl acrylate is formed followed by the CC bond at the ester side. As for the parent DA reaction, both bond formation events occur after the TS, i.e., they do not contribute to the energy barrier. What determines the barrier is the preparation process for CC bond formation, including the approach diene and dienophile, CC bond length changes and, in particular, rehybridization of the carbon atoms involved in the formation of the cyclohexene ring. This process is modified by both the BF3 catalyst and the water environment, where both work in a hand-in-hand fashion leading to the lowest energy barrier of 9.06 kcal/mol found for the catalyzed reaction R1 in aqueous solution compared to the highest energy barrier of 20.68 kcal/mol found for the non-catalyzed reaction R1 in the gas phase. The major effect of the BF3 catalyst is the increased mutual polarization and the increased charge transfer between methyl acrylate and butadiene, facilitating the approach of diene and dienophile and the pyramidalization of the CC atoms involved in the ring formation, which leads to a lowering of the activation energy. The catalytic effect of water solution is threefold. The polar environment leads also to increased polarization and charge transfer between the reacting species, similar as in the case of the BF3 catalyst, although to a smaller extend. More important is the formation of hydrogen bonds with the reaction complex, which are stronger for the TS than for the reactant, thus stabilizing the TS which leads to a further reduction of the activation energy. As shown by the ring puckering analysis, the third effect of water is space confinement of the reacting partners, conserving the boat form of the six-member ring from the entrance to the exit reaction channel. In summary, URVA combined with LMA has led to a clearer picture on how both BF3 catalyst and aqueous environment in a synchronized effort lower the reaction barrier. These new insights will serve to further fine-tune the DA reaction of methyl acrylate and butadiene and DA reactions in general.
Collapse
|
4
|
Unke OT, Koner D, Patra S, Käser S, Meuwly M. High-dimensional potential energy surfaces for molecular simulations: from empiricism to machine learning. MACHINE LEARNING-SCIENCE AND TECHNOLOGY 2020. [DOI: 10.1088/2632-2153/ab5922] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
|
5
|
Pestana LR, Hao H, Head-Gordon T. Diels-Alder Reactions in Water Are Determined by Microsolvation. NANO LETTERS 2020; 20:606-611. [PMID: 31771330 DOI: 10.1021/acs.nanolett.9b04369] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Nanoconfined aqueous environments and the recent advent of accelerated chemistry in microdroplets are increasingly being investigated for catalysis. The mechanisms underlying the enhanced reactivity in alternate solvent environments and whether the enhanced reactivity due to nanoconfinement is a universal phenomenon are not fully understood. Here, we use ab initio molecular dynamics simulations to characterize the free energy of a retro-Diels-Alder reaction in bulk water at very different densities and in water nanoconfined by parallel graphene sheets. We find that the broadly different global solvation environments accelerate the reactions to a similar degree with respect to the gas-phase reaction, with activation free energies that do not differ by more than kbT from each other. The reason for the same acceleration factor in the extremely different solvation environments is that it is the microsolvation of the dienophile's carbonyl group that governs the transition-state stabilization and mechanism, which is not significantly disrupted by either the lower density in bulk water or the strong nanoconfinement conditions used here. Our results also suggest that significant acceleration of Diels-Alder reactions in microdroplets or on-water conditions cannot arise from local microsolvation when water is present but instead must come from highly altered reaction environments that drastically change the reaction mechanisms.
Collapse
Affiliation(s)
- Luis Ruiz Pestana
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
- Pitzer Center for Theoretical Chemistry, Departments of Chemistry, Bioengineering, and Chemical and Biomolecular Engineering , University of California , Berkeley , California 94720 , United States
| | - Hongxia Hao
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
- Pitzer Center for Theoretical Chemistry, Departments of Chemistry, Bioengineering, and Chemical and Biomolecular Engineering , University of California , Berkeley , California 94720 , United States
| | - Teresa Head-Gordon
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
- Pitzer Center for Theoretical Chemistry, Departments of Chemistry, Bioengineering, and Chemical and Biomolecular Engineering , University of California , Berkeley , California 94720 , United States
| |
Collapse
|
6
|
Lodola A, Callegari D, Scalvini L, Rivara S, Mor M. Design and SAR Analysis of Covalent Inhibitors Driven by Hybrid QM/MM Simulations. Methods Mol Biol 2020; 2114:307-337. [PMID: 32016901 DOI: 10.1007/978-1-0716-0282-9_19] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Quantum mechanics/molecular mechanics (QM/MM) hybrid technique is emerging as a reliable computational method to investigate and characterize chemical reactions occurring in enzymes. From a drug discovery perspective, a thorough understanding of enzyme catalysis appears pivotal to assist the design of inhibitors able to covalently bind one of the residues belonging to the enzyme catalytic machinery. Thanks to the current advances in computer power, and the availability of more efficient algorithms for QM-based simulations, the use of QM/MM methodology is becoming a viable option in the field of covalent inhibitor design. In the present review, we summarized our experience in the field of QM/MM simulations applied to drug design problems which involved the optimization of agents working on two well-known drug targets, namely fatty acid amide hydrolase (FAAH) and epidermal growth factor receptor (EGFR). In this context, QM/MM simulations gave valuable information in terms of geometry (i.e., of transition states and metastable intermediates) and reaction energetics that allowed to correctly predict inhibitor binding orientation and substituent effect on enzyme inhibition. What is more, enzyme reaction modelling with QM/MM provided insights that were translated into the synthesis of new covalent inhibitor featured by a unique combination of intrinsic reactivity, on-target activity, and selectivity.
Collapse
Affiliation(s)
- Alessio Lodola
- Drug Design and Discovery Group, Department of Food and Drug, University of Parma, Parma, Italy.
| | - Donatella Callegari
- Drug Design and Discovery Group, Department of Food and Drug, University of Parma, Parma, Italy
| | - Laura Scalvini
- Drug Design and Discovery Group, Department of Food and Drug, University of Parma, Parma, Italy
| | - Silvia Rivara
- Drug Design and Discovery Group, Department of Food and Drug, University of Parma, Parma, Italy
| | - Marco Mor
- Drug Design and Discovery Group, Department of Food and Drug, University of Parma, Parma, Italy
| |
Collapse
|
7
|
Nguyen TH, Ngo V, Castro Zerba JP, Noskov S, Minh DDL. Nonequilibrium path-ensemble averages for symmetric protocols. J Chem Phys 2019; 151:194103. [DOI: 10.1063/1.5121306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Trung Hai Nguyen
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City, Vietnam
- Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Van Ngo
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - João Paulo Castro Zerba
- Department of Chemistry, Illinois Institute of Technology, Chicago, Illinois 60616, USA
- Institute of Physics “Gleb Wataghin”, University of Campinas, Campinas, SP 13083-859, Brazil
| | - Sergei Noskov
- SYN: Centre for Molecular Simulation, Department of Biological Sciences, University of Calgary, Calgary, Alberta T2N 1N4, Canada
| | - David D. L. Minh
- Department of Chemistry, Illinois Institute of Technology, Chicago, Illinois 60616, USA
| |
Collapse
|
8
|
Rivero U, Unke OT, Meuwly M, Willitsch S. Reactive atomistic simulations of Diels-Alder reactions: The importance of molecular rotations. J Chem Phys 2019; 151:104301. [DOI: 10.1063/1.5114981] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Uxía Rivero
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland
| | - Oliver T. Unke
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland
| | - Stefan Willitsch
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland
| |
Collapse
|
9
|
Li P, Liu F, Shao Y, Mei Y. Computational Insights into Endo/Exo Selectivity of the Diels–Alder Reaction in Explicit Solvent at Ab Initio Quantum Mechanical/Molecular Mechanical Level. J Phys Chem B 2019; 123:5131-5138. [DOI: 10.1021/acs.jpcb.9b01989] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Pengfei Li
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Fengjiao Liu
- College of Chemistry and Chemical Engineering, Shanghai University of Engineering Science, Shanghai 201620, China
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, 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
|
Affiliation(s)
- Belinda L. Slakman
- Department of Chemical Engineering Northeastern University Boston Massachusetts
| | - Richard H. West
- Department of Chemical Engineering Northeastern University Boston Massachusetts
| |
Collapse
|
12
|
Karabencheva-Christova TG, Torras J, Mulholland AJ, Lodola A, Christov CZ. Mechanistic Insights into the Reaction of Chlorination of Tryptophan Catalyzed by Tryptophan 7-Halogenase. Sci Rep 2017; 7:17395. [PMID: 29234124 PMCID: PMC5727139 DOI: 10.1038/s41598-017-17789-x] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2017] [Accepted: 11/30/2017] [Indexed: 12/14/2022] Open
Abstract
Tryptophan 7-halogenase catalyzes chlorination of free tryptophan to 7-chlorotryptophan, which is the first step in the antibiotic pyrrolnitrin biosynthesis. Many biologically and pharmaceutically active natural products contain chlorine and thus, an understanding of the mechanism of its introduction into organic molecules is important. Whilst enzyme-catalyzed chlorination is accomplished with ease, it remains a difficult task for the chemists. Therefore, utilizing enzymes in the synthesis of chlorinated organic compounds is important, and providing atomistic mechanistic insights about the reaction mechanism of tryptophan 7-halogenase is vital and timely. In this work, we examined a mechanism for the reaction of tryptophan chlorination, performed by tryptophan 7-halogenase, by calculating potential energy and free energy surfaces using two different Combined Quantum Mechanical/Molecular Mechanical (QM/MM) methods both employing Density Functional Theory (DFT) for the QM region. Both computational strategies agree on the nature of the rate-limiting step and provided close results for the reaction barriers of the two reaction steps. The calculations for both the potential energy and the free energy profiles showed very similar geometric features and hydrogen bonding interactions for the characterized stationary points.
Collapse
Affiliation(s)
- Tatyana G Karabencheva-Christova
- Department of Chemistry, Michigan Technological University, Houghton, 49931, MI, USA.
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock's Close, Bristol, BS8 1TS, UK.
| | - Juan Torras
- Department of Chemical Engineering, Escola d'Enginyeria de Barcelona Est (EEBE), Universitat Politècnica de Catalunya, C. Eduard Maristany 10-14, 08019, Barcelona, Spain.
| | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock's Close, Bristol, BS8 1TS, UK
| | - Alessio Lodola
- Pharmacy Department, Università di Parma, V. le P.G Usberti 27/A, Campus Universitario, 431124, Parma, Italy
| | - Christo Z Christov
- Department of Chemistry, Michigan Technological University, Houghton, 49931, MI, USA
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock's Close, Bristol, BS8 1TS, UK
| |
Collapse
|
13
|
Supady A, Hecht S, Baldauf C. About Underappreciated Yet Active Conformations of Thiourea Organocatalysts. Org Lett 2017; 19:4199-4202. [PMID: 28786673 DOI: 10.1021/acs.orglett.7b01782] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Conformational dynamics can define the function of organocatalysts. While the accepted mechanism of Schreiner's catalyst features a double hydrogen bond to the substrate that only forms with the anti-anti conformation of its central thiourea group, our electronic-structure theory study reveals that binding of the model substrate methyl vinyl ketone prefers syn-anti conformations. We find a new mechanism featuring π stacking interactions and highlight the need for extensive structure searches for flexible molecules, especially when aiming for structure-based design of catalytic activity.
Collapse
Affiliation(s)
- Adriana Supady
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , 14195 Berlin, Germany
| | - Stefan Hecht
- Department of Chemistry, Humboldt-Universität zu Berlin , 12489 Berlin, Germany
| | - Carsten Baldauf
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , 14195 Berlin, Germany
| |
Collapse
|