51
|
Nie X, Ren X, Ji X, Chen Y, Janik MJ, Guo X, Song C. Mechanistic Insight into Propylene Epoxidation with H 2O 2 over Titanium Silicalite-1: Effects of Zeolite Confinement and Solvent. J Phys Chem B 2019; 123:7410-7423. [PMID: 31387353 DOI: 10.1021/acs.jpcb.9b04439] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Density functional theory (DFT) calculations were performed to investigate the effects of zeolite confinement and solvent on propylene epoxidation with H2O2 over the titanium silicalite-1 (TS-1) catalyst. The 144T and 143T cluster models containing typical 10MR channels of TS-1 were constructed to represent the tripodal(2I) and Ti/defect sites. It was found that the confinement of the zeolite pore channel not only impacts the adsorption stability of guest molecules but also alters reaction barriers, as compared to the results obtained based on small cluster models. When dispersion corrections were considered, an enhancement of the adsorption stability of guest molecules was observed because of the important contribution from van der Waals interactions, especially for propylene adsorption. An explicit protic methanol molecule was introduced into the catalytic system to probe the influence of the solvent on propylene epoxidation, based on which a significant enhancement of CH3OH-H2O2 co-adsorption was obtained owing to H-bond formation. More importantly, the energy barrier for H2O2 dissociation was largely reduced by ∼13 kcal/mol because of the participation of the methanol in the H-transfer process and the formation of H-bond network, resulting in an alteration of the rate-limiting step. By comparison, adding an aprotic acetonitrile solvent did not have substantial effect on reaction path and kinetics. The calculation results clearly demonstrate the important role of the protic methanol solvent, which not only strengthens the adsorption of guest molecules but also promotes the kinetics for propylene epoxidation with H2O2 over TS-1 catalyst.
Collapse
|
52
|
Metal-organic structures with formate and sulfate anions: Synthesis, crystallographic studies and hydrogen storage by PM7 and ONIOM. J Mol Struct 2019. [DOI: 10.1016/j.molstruc.2019.03.102] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
53
|
Sanusi ZK, Lawal MM, Govender T, Maguire GEM, Honarparvar B, Kruger HG. Theoretical Model for HIV-1 PR That Accounts for Substrate Recognition and Preferential Cleavage of Natural Substrates. J Phys Chem B 2019; 123:6389-6400. [DOI: 10.1021/acs.jpcb.9b02207] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Zainab K. Sanusi
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4041, South Africa
| | - Monsurat M. Lawal
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4041, South Africa
| | | | - Glenn E. M. Maguire
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4041, South Africa
- School of Chemistry and Physics, University of KwaZulu-Natal, Durban 4041, South Africa
| | - Bahareh Honarparvar
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4041, South Africa
| | - Hendrik G. Kruger
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4041, South Africa
| |
Collapse
|
54
|
Abstract
Complex chemical systems present challenges to electronic structure theory stemming from large system sizes, subtle interactions, coupled dynamical time scales, and electronically nonadiabatic effects. New methods are needed to perform reliable, rigorous, and affordable electronic structure calculations for simulating the properties and dynamics of such systems. This Account reviews projection-based quantum embedding for electronic structure, which provides a formally exact method for density functional theory (DFT) embedding. The method also provides a rigorous and accurate approach for describing a small part of a chemical system at the level of a correlated wavefunction (WF) method while the remainder of the system is described at the level of DFT. A key advantage of projection-based embedding is that it can be formulated in terms of an extremely simple level-shift projection operator, which eliminates the need for any optimized effective potential calculation or kinetic energy functional approximation while simultaneously ensuring that no extra programming is needed to perform WF-in-DFT embedding with an arbitrary WF method. The current work presents the theoretical underpinnings of projection-based embedding, describes use of the method for combining wavefunction and density functional theories, and discusses technical refinements that have improved the applicability and robustness of the method. Applications of projection-based WF-in-DFT embedding are also reviewed, with particular focus on recent work on transition-metal catalysis, enzyme reactivity, and battery electrolyte decomposition. In particular, we review the application of projection-based embedding for the prediction of electrochemical potentials and reaction pathways in a Co-centered hydrogen evolution catalyst. Projection-based WF-in-DFT calculations are shown to provide quantitative accuracy while greatly reducing the computational cost compared with a reference coupled cluster calculation on the full system. Additionally, projection-based WF-in-DFT embedding is used to study the mechanism of citrate synthase; it is shown that projection-based WF-in-DFT largely eliminates the sensitivity of the potential energy landscape to the employed DFT exchange-correlation functional. Finally, we demonstrate the use of projection-based WF-in-DFT to study electron transfer reactions associated with battery electrolyte decomposition. Projection-based WF-in-DFT embedding is used to calculate the oxidation potentials of neat ethylene carbonate (EC), neat dimethyl carbonate (DMC), and 1:1 mixtures of EC and DMC in order to overcome qualitative inaccuracies in the electron densities and ionization energies obtained from conventional DFT methods. By further embedding the WF-in-DFT description in a molecular mechanics point-charge environment, this work enables an explicit description of the solvent and ensemble averaging of the solvent configurations. Looking forward, we anticipate continued refinement of the projection-based embedding methodology as well as its increasingly widespread application in diverse areas of chemistry, biology, and materials science.
Collapse
Affiliation(s)
- Sebastian J. R. Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Matthew Welborn
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Frederick R. Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Thomas F. Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
55
|
Quesne MG, Silveri F, de Leeuw NH, Catlow CRA. Advances in Sustainable Catalysis: A Computational Perspective. Front Chem 2019; 7:182. [PMID: 31032245 PMCID: PMC6473102 DOI: 10.3389/fchem.2019.00182] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2019] [Accepted: 03/07/2019] [Indexed: 11/13/2022] Open
Abstract
The enormous challenge of moving our societies to a more sustainable future offers several exciting opportunities for computational chemists. The first principles approach to "catalysis by design" will enable new and much greener chemical routes to produce vital fuels and fine chemicals. This prospective outlines a wide variety of case studies to underscore how the use of theoretical techniques, from QM/MM to unrestricted DFT and periodic boundary conditions, can be applied to biocatalysis and to both homogeneous and heterogenous catalysts of all sizes and morphologies to provide invaluable insights into the reaction mechanisms they catalyze.
Collapse
|
56
|
Debnath S, Sengupta A, Jose KVJ, Raghavachari K. Fragment-Based Approaches for Supramolecular Interaction Energies: Applications to Foldamers and Their Complexes with Anions. J Chem Theory Comput 2018; 14:6226-6239. [DOI: 10.1021/acs.jctc.8b00525] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Sibali Debnath
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Arkajyoti Sengupta
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - K. V. Jovan Jose
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
57
|
Welborn M, Manby FR, Miller TF. Even-handed subsystem selection in projection-based embedding. J Chem Phys 2018; 149:144101. [DOI: 10.1063/1.5050533] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Affiliation(s)
- Matthew Welborn
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Frederick R. Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Thomas F. Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
58
|
Do Better Quality Embedding Potentials Accelerate the Convergence of QM/MM Models? The Case of Solvated Acid Clusters. Molecules 2018; 23:molecules23102466. [PMID: 30261616 PMCID: PMC6222911 DOI: 10.3390/molecules23102466] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2018] [Revised: 09/22/2018] [Accepted: 09/24/2018] [Indexed: 12/28/2022] Open
Abstract
This study examines whether the use of more accurate embedding potentials improves the convergence of quantum mechanics/molecular mechanics (QM/MM) models with respect to the size of the QM region. In conjunction with density functional theory calculations using the ωB97X-D functional, various embedding potentials including the TIP3P water model, the effective fragment potential (EFP), and semi-empirical methods (PM6, PM7, and DFTB) were used to simulate the deprotonation energies of solvated acid clusters. The calculations were performed on solvated neutral (HA) and cationic (HB⁺) acids clusters containing 160 and 480 water molecules using configurations sampled from molecular dynamics simulations. Consistently, the ωB97X-D/EFP model performed the best when using a minimal QM region size. The performance for the other potentials appears to be highly sensitive to the charge character of the acid/base pair. Neutral acids display the expected trend that semi-empirical methods generally perform better than TIP3P; however, an opposite trend was observed for the cationic acids. Additionally, electronic embedding provided an improvement over mechanical embedding for the cationic systems, but not the neutral acids. For the best performing ωB97X-D/EFP model, a QM region containing about 6% of the total number of solvent molecules is needed to approach within 10 kJ mol-1 of the pure QM result if the QM region was chosen based on the distance from the reaction centre.
Collapse
|
59
|
Sameera WMC, Maseras F. Expanding the Range of Force Fields Available for ONIOM Calculations: The SICTWO Interface. J Chem Inf Model 2018; 58:1828-1835. [PMID: 30137980 DOI: 10.1021/acs.jcim.8b00332] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
The ONIOM scheme is one of the most popular QM/MM approaches, but its extended application has been so far hindered by the limited availability of force fields in most practical implementations. This paper describes a simple software code to overcome this limitation, and its application to three representative chemical problems. The "Shell Interface for Combining Tinker With ONIOM" (SICTWO) program gives access to all force fields available in the Tinker molecular mechanics program from a Gaussian09 or Gaussian16 calculation. The first application presented is the geometry optimization of dithienobicyclo-[4.4.1]-undeca-3,8-diene-11-one ethylene glycol ketal. A variety of force fields were tested for the MM part, and the molecular structure using the ONIOM(B3LYP:OPLS-AA) method shows good agreement with the experimental results. The second problem is related with enantioselectivity of the vanadium-catalyzed asymmetric oxidation of 1,2-bis( tert-butyl)-disulfide. ONIOM(B3LYP:MM3) results, obtained with SICTWO, are consistent with those of a previous IMOMM(B3LYP:MM3) study. In the third application, we have combined AMOEBA09 polarizable force field with ONIOM to study a water molecule binding on water ice. Calculated ONIOM(ωB97X-D:AMOEBA09) binding energies are consistent with the ωB97X-D results. These studies show SICTWO as an easy-to-use tool that can further expand the use of the multiscale ONIOM method within computational chemistry and computational biology.
Collapse
Affiliation(s)
- W M C Sameera
- Institute of Chemical Research of Catalonia (ICIQ) , The Barcelona Institute of Science and Technology , Avgda. Països Catalans 16 , 43007 , Tarragona , Spain.,Institute of Low Temperature Science , Hokkaido University , Sapporo , Hokkaido 060-0819 , Japan
| | - Feliu Maseras
- Institute of Chemical Research of Catalonia (ICIQ) , The Barcelona Institute of Science and Technology , Avgda. Països Catalans 16 , 43007 , Tarragona , Spain.,Departament de Química , Universitat Autònoma de Barcelona , 08193 Bellaterra , Spain
| |
Collapse
|
60
|
S Fernandes H, Ramos MJ, M F S A Cerqueira N. molUP: A VMD plugin to handle QM and ONIOM calculations using the gaussian software. J Comput Chem 2018; 39:1344-1353. [PMID: 29464735 DOI: 10.1002/jcc.25189] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2017] [Revised: 01/12/2018] [Accepted: 01/30/2018] [Indexed: 11/08/2022]
Abstract
The notable advances obtained by computational (bio)chemistry provided its widespread use in many areas of science, in particular, in the study of reaction mechanisms. These studies involve a huge number of complex calculations, which are often carried out using the Gaussian suite of programs. The preparation of input files and the analysis of the output files are not easy tasks and often involve laborious and complex steps. Taking this into account, we developed molUP: a VMD plugin that offers a complete set of tools that enhance the preparation of QM and ONIOM (QM/MM, QM/QM, and QM/QM/MM) calculations. The starting structures for these calculations can be imported from different chemical formats. A set of tools is available to help the user to examine or modify any geometry parameter. This includes the definition of layers in ONIOM calculations, choosing fixed atoms during geometry optimizations, the recalculation or adjustment of the atomic charges, performing SCANs or IRC calculations, etc. molUP also extracts the geometries from the output files as well as the energies of each of them. All of these tasks are performed in an interactive GUI that is extremely helpful for the user. MolUP was developed to be easy to handle by inexperienced users, but simultaneously to be a fast and flexible graphical interface to allow the advanced users to take full advantage of this plugin. The program is available, free of charges, for macOS, Linux, and Windows at the PortoBioComp page https://www.fc.up.pt/PortoBioComp/database/doku.php?id=molup. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Henrique S Fernandes
- UCIBIO@REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências da Universidade do Porto, Rua do Campo Alegre s/n, Porto, 4169-007, Portugal
| | - Maria J Ramos
- UCIBIO@REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências da Universidade do Porto, Rua do Campo Alegre s/n, Porto, 4169-007, Portugal
| | - Nuno M F S A Cerqueira
- UCIBIO@REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências da Universidade do Porto, Rua do Campo Alegre s/n, Porto, 4169-007, Portugal
| |
Collapse
|
61
|
Sanusi ZK, Govender T, Maguire GEM, Maseko SB, Lin J, Kruger HG, Honarparvar B. An insight to the molecular interactions of the FDA approved HIV PR drugs against L38L↑N↑L PR mutant. J Comput Aided Mol Des 2018; 32:459-471. [PMID: 29397520 DOI: 10.1007/s10822-018-0099-9] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2017] [Accepted: 01/16/2018] [Indexed: 01/12/2023]
Abstract
The aspartate protease of the human immune deficiency type-1 virus (HIV-1) has become a crucial antiviral target in which many useful antiretroviral inhibitors have been developed. However, it seems the emergence of new HIV-1 PR mutations enhances drug resistance, hence, the available FDA approved drugs show less activity towards the protease. A mutation and insertion designated L38L↑N↑L PR was recently reported from subtype of C-SA HIV-1. An integrated two-layered ONIOM (QM:MM) method was employed in this study to examine the binding affinities of the nine HIV PR inhibitors against this mutant. The computed binding free energies as well as experimental data revealed a reduced inhibitory activity towards the L38L↑N↑L PR in comparison with subtype C-SA HIV-1 PR. This observation suggests that the insertion and mutations significantly affect the binding affinities or characteristics of the HIV PIs and/or parent PR. The same trend for the computational binding free energies was observed for eight of the nine inhibitors with respect to the experimental binding free energies. The outcome of this study shows that ONIOM method can be used as a reliable computational approach to rationalize lead compounds against specific targets. The nature of the intermolecular interactions in terms of the host-guest hydrogen bond interactions is discussed using the atoms in molecules (AIM) analysis. Natural bond orbital analysis was also used to determine the extent of charge transfer between the QM region of the L38L↑N↑L PR enzyme and FDA approved drugs. AIM analysis showed that the interaction between the QM region of the L38L↑N↑L PR and FDA approved drugs are electrostatic dominant, the bond stability computed from the NBO analysis supports the results from the AIM application. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide information that will aid in the design of much improved HIV-1 PR antiviral drugs.
Collapse
Affiliation(s)
- Zainab K Sanusi
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Thavendran Govender
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Glenn E M Maguire
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa.,School of Chemistry and Physics, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Sibusiso B Maseko
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Johnson Lin
- School of Life Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Hendrik G Kruger
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa.
| | - Bahareh Honarparvar
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa.
| |
Collapse
|
62
|
Sharma A, Delile S, Jabri M, Adamo C, Fave C, Marchal D, Perrier A. Interaction of osmium(ii) redox probes with DNA: insights from theory. Phys Chem Chem Phys 2018; 18:30029-30039. [PMID: 27774536 DOI: 10.1039/c6cp05105g] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
In the course of developing ultrasensitive and quantitative electrochemical point-of-care analytical tools for genetic detection of infectious diseases, osmium(ii) metallointercalators were revealed to be suitable and efficient redox probes to monitor the in vitro DNA amplification [Defever etal, Anal. Chem., 2011, 83, 1815-1821]. In this work, we thus propose a complete computational protocol in order to evaluate the affinity between Os(ii) complexes with double-stranded DNA. This protocol is based on molecular dynamics, with the parametrization of the GAFF force field for the Os(ii) complexes presenting an octahedral environment with polypyridine ligands, and QM/QM' calculations to evaluate the binding energy. For three Os(ii) probes and different binding sites, molecular dynamics simulations and interaction energies calculated at the QM/QM' level are successively discussed and compared to experimental data in order to identify the most stable binding sites. The computational protocol we propose should then be used to design more efficient Os(ii) metallointercalators.
Collapse
Affiliation(s)
- Ashwani Sharma
- Chimie ParisTech, PSL Research University, CNRS, Institut de Recherche de Chimie Paris (IRCP), F-75005 Paris, France and Laboratoire d'Electrochimie Moléculaire, UMR 7591 CNRS, Université Paris Diderot, Sorbonne Paris Cité, 15 rue J-A de Baif, F-75205 Paris Cedex 13, France.
| | - Sebastien Delile
- Laboratoire d'Electrochimie Moléculaire, UMR 7591 CNRS, Université Paris Diderot, Sorbonne Paris Cité, 15 rue J-A de Baif, F-75205 Paris Cedex 13, France.
| | - Mohamed Jabri
- Chimie ParisTech, PSL Research University, CNRS, Institut de Recherche de Chimie Paris (IRCP), F-75005 Paris, France and E-pôle de génoinformatique, Institut Jacques Monod, UMR7592, CNRS, Université Paris Diderot, Sorbonne Paris Cité, F-75013 Paris, France
| | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institut de Recherche de Chimie Paris (IRCP), F-75005 Paris, France and Institut Universitaire de France, 103 Boulevard Saint Michel, F-75005 Paris, France
| | - Claire Fave
- Laboratoire d'Electrochimie Moléculaire, UMR 7591 CNRS, Université Paris Diderot, Sorbonne Paris Cité, 15 rue J-A de Baif, F-75205 Paris Cedex 13, France.
| | - Damien Marchal
- Laboratoire d'Electrochimie Moléculaire, UMR 7591 CNRS, Université Paris Diderot, Sorbonne Paris Cité, 15 rue J-A de Baif, F-75205 Paris Cedex 13, France.
| | - Aurélie Perrier
- Chimie ParisTech, PSL Research University, CNRS, Institut de Recherche de Chimie Paris (IRCP), F-75005 Paris, France and Université Paris Diderot, Sorbonne Paris Cité, 5 rue Thomas Mann, F-75205 Paris Cedex 13, France.
| |
Collapse
|
63
|
Mann D, Güldenhaupt J, Schartner J, Gerwert K, Kötting C. The protonation states of GTP and GppNHp in Ras proteins. J Biol Chem 2018; 293:3871-3879. [PMID: 29382720 DOI: 10.1074/jbc.ra117.001110] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2017] [Revised: 01/29/2018] [Indexed: 01/09/2023] Open
Abstract
The small GTPase Ras transmits signals in a variety of cellular signaling pathways, most prominently in cell proliferation. GTP hydrolysis in the active center of Ras acts as a prototype for many GTPases and is the key to the understanding of several diseases, including cancer. Therefore, Ras has been the focus of intense research over the last decades. A recent neutron diffraction crystal structure of Ras indicated a protonated γ-guanylyl imidodiphosphate (γ-GppNHp) group, which has put the protonation state of GTP in question. A possible protonation of GTP was not considered in previously published mechanistic studies. To determine the detailed prehydrolysis state of Ras, we calculated infrared and NMR spectra from quantum mechanics/molecular mechanics (QM/MM) simulations and compared them with those from previous studies. Furthermore, we measured infrared spectra of GTP and several GTP analogs bound to lipidated Ras on a membrane system under near-native conditions. Our findings unify results from previous studies and indicate a structural model confirming the hypothesis that γ-GTP is fully deprotonated in the prehydrolysis state of Ras.
Collapse
Affiliation(s)
- Daniel Mann
- From the Department of Biophysics, Ruhr University Bochum, 44780 Bochum, Germany and
| | - Jörn Güldenhaupt
- From the Department of Biophysics, Ruhr University Bochum, 44780 Bochum, Germany and
| | - Jonas Schartner
- From the Department of Biophysics, Ruhr University Bochum, 44780 Bochum, Germany and
| | - Klaus Gerwert
- From the Department of Biophysics, Ruhr University Bochum, 44780 Bochum, Germany and .,Max-Planck-Gesellschaft-Chinese Academy of Sciences (MPG-CAS) Partner Institute for Computational Biology (PICB), Shanghai 200031, China
| | - Carsten Kötting
- From the Department of Biophysics, Ruhr University Bochum, 44780 Bochum, Germany and
| |
Collapse
|
64
|
De S, C. H. R, Thamleena A. H, Joseph A, Ben A, V. U. K. Roles of different amino-acid residues towards binding and selective transport of K+ through KcsA K+-ion channel. Phys Chem Chem Phys 2018; 20:17517-17529. [DOI: 10.1039/c8cp01282b] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Each amino acid in the selectivity filter plays a distinct role towards binding and transport of K+ ion through KcsA.
Collapse
Affiliation(s)
- Susmita De
- Department of Applied Chemistry
- Cochin University of Science and Technology
- Trikakkara
- Kochi
- India – 682 022
| | - Rinsha C. H.
- Theoretical and Computational Chemistry Laboratory
- Department of Chemistry
- National Institute of Technology Calicut
- Kozhikode
- India – 673 601
| | - Hanna Thamleena A.
- Theoretical and Computational Chemistry Laboratory
- Department of Chemistry
- National Institute of Technology Calicut
- Kozhikode
- India – 673 601
| | - Annu Joseph
- Theoretical and Computational Chemistry Laboratory
- Department of Chemistry
- National Institute of Technology Calicut
- Kozhikode
- India – 673 601
| | - Anju Ben
- Theoretical and Computational Chemistry Laboratory
- Department of Chemistry
- National Institute of Technology Calicut
- Kozhikode
- India – 673 601
| | - Krishnapriya V. U.
- Theoretical and Computational Chemistry Laboratory
- Department of Chemistry
- National Institute of Technology Calicut
- Kozhikode
- India – 673 601
| |
Collapse
|
65
|
Gozem S, Luk HL, Schapiro I, Olivucci M. Theory and Simulation of the Ultrafast Double-Bond Isomerization of Biological Chromophores. Chem Rev 2017; 117:13502-13565. [DOI: 10.1021/acs.chemrev.7b00177] [Citation(s) in RCA: 175] [Impact Index Per Article: 25.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Samer Gozem
- Department
of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States
| | - Hoi Ling Luk
- Chemistry
Department, Bowling Green State University, Overman Hall, Bowling Green, Ohio 43403, United States
| | - Igor Schapiro
- Fritz
Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Massimo Olivucci
- Chemistry
Department, Bowling Green State University, Overman Hall, Bowling Green, Ohio 43403, United States
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro
2, 53100 Siena, Italy
| |
Collapse
|
66
|
Libisch F, Marsman M, Burgdörfer J, Kresse G. Embedding for bulk systems using localized atomic orbitals. J Chem Phys 2017; 147:034110. [DOI: 10.1063/1.4993795] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- F. Libisch
- Institute for Theoretical Physics, Vienna University of Technology, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria, EU
| | - M. Marsman
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Sensengasse 8/12, A-1090 Vienna, Austria, EU
| | - J. Burgdörfer
- Institute for Theoretical Physics, Vienna University of Technology, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria, EU
| | - G. Kresse
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Sensengasse 8/12, A-1090 Vienna, Austria, EU
| |
Collapse
|
67
|
Sanusi ZK, Govender T, Maguire GEM, Maseko SB, Lin J, Kruger HG, Honarparvar B. Investigation of the binding free energies of FDA approved drugs against subtype B and C-SA HIV PR: ONIOM approach. J Mol Graph Model 2017; 76:77-85. [PMID: 28711760 DOI: 10.1016/j.jmgm.2017.06.026] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2017] [Revised: 06/08/2017] [Accepted: 06/09/2017] [Indexed: 01/15/2023]
Abstract
Human immune virus subtype C is the most widely spread HIV subtype in Sub-Sahara Africa and South Africa. A profound structural insight on finding potential lead compounds is therefore necessary for drug discovery. The focus of this study is to rationalize the nine Food and Drugs Administration (FDA) HIV antiviral drugs complexed to subtype B and C-SA PR using ONIOM approach. To achieve this, an integrated two-layered ONIOM model was used to optimize the geometrics of the FDA approved HIV-1 PR inhibitors for subtype B. In our hybrid ONIOM model, the HIV-1 PR inhibitors as well as the ASP 25/25' catalytic active residues were treated at high level quantum mechanics (QM) theory using B3LYP/6-31G(d), and the remaining HIV PR residues were considered using the AMBER force field. The experimental binding energies of the PR inhibitors were compared to the ONIOM calculated results. The theoretical binding free energies (?Gbind) for subtype B follow a similar trend to the experimental results, with one exemption. The computational model was less suitable for C-SA PR. Analysis of the results provided valuable information about the shortcomings of this approach. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide much improved binding energies for complex enzyme drug interactions.
Collapse
Affiliation(s)
- Z K Sanusi
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4001, South Africa
| | - T Govender
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4001, South Africa
| | - G E M Maguire
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4001, South Africa; School of Chemistry and Physics, University of KwaZulu-Natal, 4001 Durban, South Africa
| | - S B Maseko
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4001, South Africa
| | - J Lin
- School of Life Sciences, University of KwaZulu-Natal, Durban 4001, South Africa
| | - H G Kruger
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4001, South Africa.
| | - B Honarparvar
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4001, South Africa.
| |
Collapse
|
68
|
Elucidation of Single Hydrogen Bonds in GTPases via Experimental and Theoretical Infrared Spectroscopy. Biophys J 2017; 112:66-77. [PMID: 28076817 PMCID: PMC5232353 DOI: 10.1016/j.bpj.2016.11.3195] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2016] [Revised: 11/04/2016] [Accepted: 11/28/2016] [Indexed: 11/28/2022] Open
Abstract
Time-resolved Fourier transform infrared (FTIR) spectroscopy is a powerful tool to elucidate label-free the reaction mechanisms of proteins. After assignment of the absorption bands to individual groups of the protein, the order of events during the reaction mechanism can be monitored and rate constants can be obtained. Additionally, structural information is encoded into infrared spectra and can be decoded by combining the experimental data with biomolecular simulations. We have determined recently the infrared vibrations of GTP and guanosine diphosphate (GDP) bound to Gαi1, a ubiquitous GTPase. These vibrations are highly sensitive for the environment of the phosphate groups and thereby for the binding mode the GTPase adopts to enable fast hydrolysis of GTP. In this study we calculated these infrared vibrations from biomolecular simulations to transfer the spectral information into a computational model that provides structural information far beyond crystal structure resolution. Conformational ensembles were generated using 15 snapshots of several 100 ns molecular-mechanics/molecular-dynamics (MM-MD) simulations, followed by quantum-mechanics/molecular-mechanics (QM/MM) minimization and normal mode analysis. In comparison with other approaches, no time-consuming QM/MM-MD simulation was necessary. We carefully benchmarked the simulation systems by deletion of single hydrogen bonds between the GTPase and GTP through several Gαi1 point mutants. The missing hydrogen bonds lead to blue-shifts of the corresponding absorption bands. These band shifts for α-GTP (Gαi1-T48A), γ-GTP (Gαi1-R178S), and for both β-GTP/γ-GTP (Gαi1-K46A, Gαi1-D200E) were found in agreement in the experimental and the theoretical spectra. We applied our approach to open questions regarding Gαi1: we show that the GDP state of Gαi1 carries a Mg2+, which is not found in x-ray structures. Further, the catalytic role of K46, a central residue of the P-loop, and the protonation state of the GTP are elucidated.
Collapse
|
69
|
Beckett D, Krukau A, Raghavachari K. Charge redistribution in QM:QM ONIOM model systems: a constrained density functional theory approach. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1333643] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Daniel Beckett
- Department of Chemistry, Indiana University, Bloomington, IN, USA
| | | | | |
Collapse
|
70
|
Farberow CA, Cheah S, Kim S, Miller JT, Gallagher JR, Hensley JE, Schaidle JA, Ruddy DA. Exploring Low-Temperature Dehydrogenation at Ionic Cu Sites in Beta Zeolite To Enable Alkane Recycle in Dimethyl Ether Homologation. ACS Catal 2017. [DOI: 10.1021/acscatal.6b03582] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Carrie A. Farberow
- National
Bioenergy Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Singfoong Cheah
- National
Bioenergy Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Seonah Kim
- National
Bioenergy Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Jeffrey T. Miller
- School
of Chemical Engineering, Purdue University, West Lafayette, Indiana 47907, United States
| | - James R. Gallagher
- Chemical
Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Jesse E. Hensley
- National
Bioenergy Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
- Department
of Chemistry, Colorado School of Mines, Golden, Colorado 80401, United States
| | - Joshua A. Schaidle
- National
Bioenergy Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Daniel A. Ruddy
- National
Bioenergy Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
- Department
of Chemistry, Colorado School of Mines, Golden, Colorado 80401, United States
| |
Collapse
|
71
|
Pennifold RCR, Bennie SJ, Miller TF, Manby FR. Correcting density-driven errors in projection-based embedding. J Chem Phys 2017; 146:084113. [PMID: 28249446 DOI: 10.1063/1.4974929] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Projection-based embedding provides a simple and numerically robust framework for multiscale wavefunction-in-density-functional-theory (WF-in-DFT) calculations. The approach works well when the approximate DFT is sufficiently accurate to describe the energetics of the low-level subsystem and the coupling between subsystems. It is also necessary that the low-level DFT produces a qualitatively reasonable description of the total density, and in this work, we study model systems where delocalization error prevents this from being the case. We find substantial errors in embedding calculations on open-shell doublet systems in which self-interaction errors cause spurious delocalization of the singly occupied orbital. We propose a solution to this error by evaluating the DFT energy using a more accurate self-consistent density, such as that of Hartree-Fock (HF) theory. These so-called WF-in-(HF-DFT) calculations show excellent convergence towards full-system wavefunction calculations.
Collapse
Affiliation(s)
- Robert C R Pennifold
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Simon J Bennie
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Thomas F Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Frederick R Manby
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
72
|
Salavati-fard T, Caratzoulas S, Lobo RF, Doren DJ. Catalysis of the Diels–Alder Reaction of Furan and Methyl Acrylate in Lewis Acidic Zeolites. ACS Catal 2017. [DOI: 10.1021/acscatal.6b02682] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Taha Salavati-fard
- Department
of Physics and Astronomy, ‡Catalysis Center for Energy Innovation
(CCEI), Department of Chemical and Biomolecular Engineering, and §Department of
Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, United States
| | - Stavros Caratzoulas
- Department
of Physics and Astronomy, ‡Catalysis Center for Energy Innovation
(CCEI), Department of Chemical and Biomolecular Engineering, and §Department of
Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, United States
| | - Raul F. Lobo
- Department
of Physics and Astronomy, ‡Catalysis Center for Energy Innovation
(CCEI), Department of Chemical and Biomolecular Engineering, and §Department of
Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, United States
| | - Douglas J. Doren
- Department
of Physics and Astronomy, ‡Catalysis Center for Energy Innovation
(CCEI), Department of Chemical and Biomolecular Engineering, and §Department of
Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
73
|
Jose KVJ, Raghavachari K. Fragment-Based Approach for the Evaluation of NMR Chemical Shifts for Large Biomolecules Incorporating the Effects of the Solvent Environment. J Chem Theory Comput 2017; 13:1147-1158. [DOI: 10.1021/acs.jctc.6b00922] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- K. V. Jovan Jose
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
74
|
Biancardi A, Barnes J, Caricato M. Point charge embedding for ONIOM excited states calculations. J Chem Phys 2017; 145:224109. [PMID: 27984901 DOI: 10.1063/1.4972000] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Hybrid quantum mechanical methods can assist in the interpretation and prediction of the electronic spectra of large molecular structures. In this work, we study the performance of the ONIOM (Our own N-layered Integrated molecular Orbital molecular Mechanics) hybrid method for the calculation of transition energies and oscillator strengths by embedding the core region in a field of fixed point charges. These charges introduce polarization effects from the substituent groups to the core region. We test various charge definitions, with particular attention to the issue of overpolarization near the boundary between layers. To minimize this issue, we fit the charges on the electrostatic potential of the entire structure in the presence of the link atoms used to cap dangling bonds. We propose two constrained fitting strategies: one that produces an average set of charges common to both model system calculations, EE(L1), and one that produces two separate sets of embedding charges, EE(L2). The results from our tests show that indeed electronic embedding with constrained-fitted charges tends to improve the performance of ONIOM compared to non-embedded calculations. However, the EE(L2) charges work best for transition energies, and the EE(L1) charges work best for oscillator strengths. This may be an indication that fixed point charges do not have enough flexibility to adapt to each system, and other effects (e.g., polarization of the embedding field) may be necessary.
Collapse
Affiliation(s)
- Alessandro Biancardi
- Department of Chemistry, University of Kansas, 1251 Wescoe Hall Drive, Lawrence, Kansas 66045, USA
| | - Jeremy Barnes
- Department of Chemistry, University of Kansas, 1251 Wescoe Hall Drive, Lawrence, Kansas 66045, USA
| | - Marco Caricato
- Department of Chemistry, University of Kansas, 1251 Wescoe Hall Drive, Lawrence, Kansas 66045, USA
| |
Collapse
|
75
|
Pirillo J, Mazzone G, Russo N, Bertini L. Photophysical Properties of S, Se and Te-Substituted Deoxyguanosines: Insight into Their Ability To Act as Chemotherapeutic Agents. J Chem Inf Model 2017; 57:234-242. [DOI: 10.1021/acs.jcim.6b00486] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jenny Pirillo
- Dipartimento
di Chimica e Tecnologie Chimiche, University of Calabria, Via P. Bucci, 87036 Rende, Italy
| | - Gloria Mazzone
- Dipartimento
di Chimica e Tecnologie Chimiche, University of Calabria, Via P. Bucci, 87036 Rende, Italy
| | - Nino Russo
- Dipartimento
di Chimica e Tecnologie Chimiche, University of Calabria, Via P. Bucci, 87036 Rende, Italy
| | - Luca Bertini
- Department
of Biotechnologies and Biosciences, University of Milano-Bicocca, Piazza
della Scienza 2, I-20126 Milano, Italy
| |
Collapse
|
76
|
Cerezo J, Petrone A, Ferrer FJA, Donati G, Santoro F, Improta R, Rega N. Electronic spectroscopy of a solvatochromic dye in water: comparison of static cluster/implicit and dynamical/explicit solvent models on structures and energies. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-2009-3] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
77
|
Lu J, Hu L, Cheng J, Fang D, Wang C, Yu K, Jiang H, Cui Q, Xu Y, Luo C. A computational investigation on the substrate preference of ten-eleven-translocation 2 (TET2). Phys Chem Chem Phys 2016; 18:4728-38. [PMID: 26799843 DOI: 10.1039/c5cp07266b] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
TET proteins iteratively convert 5-methylcytosine (5mC) into 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC) in a Fe(ii)/α-ketoglutarate-dependent manner. Our previous biochemical studies revealed that TET proteins are more active on 5mC than on 5hmC and 5fC. However, the source of the substrate preference of TET proteins still remains largely elusive. Here, we investigated the substrate binding and catalytic mechanisms of oxidation reactions mediated by TET2 on different substrates through computational approaches. In accordance with previous experimental reports, our computational results suggest that TET2 can bind to different substrates with comparable binding affinities and the hydrogen abstraction step in the catalytic cycle acts as the rate-limiting step. Further structural characterization of the intermediate structures revealed that the 5-substitution groups on 5hmC and 5fC adopt an unfavorable orientation for hydrogen abstraction, which leads to a higher energy barrier for 5hmC and 5fC (compared to 5mC) and thus a lower catalytic efficiency. In summary, our mechanical insights demonstrate that substrate preference is the intrinsic property of TET proteins and our theoretical calculation results can guide further dry-lab or wet-lab studies on the catalytic mechanism of TET proteins as well as other Fe(ii)/α-ketoglutarate (KG)-dependent dioxygenases.
Collapse
Affiliation(s)
- Junyan Lu
- Drug Discovery and Design Centre, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China.
| | - Lulu Hu
- Fudan University Shanghai Cancer Centre, Institute of Biomedical Sciences, Shanghai Medical College of Fudan University, Shanghai 200032, China. and State Key Laboratory of Genetic Engineering, Collaborative Innovation Centre of Genetics and Development, School of Life Sciences, Fudan University, Shanghai 200433, China
| | - Jingdong Cheng
- Fudan University Shanghai Cancer Centre, Institute of Biomedical Sciences, Shanghai Medical College of Fudan University, Shanghai 200032, China.
| | - Dong Fang
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison, 1101 University Ave, Madison, WI 53706, USA
| | - Chen Wang
- Drug Discovery and Design Centre, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China.
| | - Kunqian Yu
- Drug Discovery and Design Centre, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China.
| | - Hualiang Jiang
- Drug Discovery and Design Centre, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China.
| | - Qiang Cui
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison, 1101 University Ave, Madison, WI 53706, USA
| | - Yanhui Xu
- Fudan University Shanghai Cancer Centre, Institute of Biomedical Sciences, Shanghai Medical College of Fudan University, Shanghai 200032, China. and State Key Laboratory of Genetic Engineering, Collaborative Innovation Centre of Genetics and Development, School of Life Sciences, Fudan University, Shanghai 200433, China
| | - Cheng Luo
- Drug Discovery and Design Centre, State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China.
| |
Collapse
|
78
|
Sharma D, Sameera WMC, Andersson S, Nyman G, Paterson MJ. Computational Study of the Interactions between Benzene and Crystalline Ice Ih: Ground and Excited States. Chemphyschem 2016; 17:4079-4089. [DOI: 10.1002/cphc.201600660] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2016] [Indexed: 01/06/2023]
Affiliation(s)
- Divya Sharma
- Institute of Chemical Sciences; School of Engineering and Physical Sciences; Heriot Watt University Edinburgh EH14 4AS United Kingdom
| | - W. M. C. Sameera
- Fukui Institute for Fundamental Chemistry; Kyoto University; Kyoto 606-8103 Japan
- University of Gothenburg; Department of Chemistry and Molecular Biology; Kemigården 4 412 96 Gothenburg Sweden
| | - Stefan Andersson
- University of Gothenburg; Department of Chemistry and Molecular Biology; Kemigården 4 412 96 Gothenburg Sweden
- SINTEF Materials and Chemistry; P.O. Box 4760 7465 Trondheim Norway
| | - Gunnar Nyman
- University of Gothenburg; Department of Chemistry and Molecular Biology; Kemigården 4 412 96 Gothenburg Sweden
| | - Martin J. Paterson
- Institute of Chemical Sciences; School of Engineering and Physical Sciences; Heriot Watt University Edinburgh EH14 4AS United Kingdom
| |
Collapse
|
79
|
Zanetti-Polzi L, Corni S, Daidone I, Amadei A. Extending the essential dynamics analysis to investigate molecular properties: application to the redox potential of proteins. Phys Chem Chem Phys 2016; 18:18450-9. [PMID: 27339768 DOI: 10.1039/c6cp03394f] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Here, a methodology is proposed to investigate the collective fluctuation modes of an arbitrary set of observables, maximally contributing to the fluctuation of another functionally relevant observable. The methodology, based on the analysis of fully classical molecular dynamics (MD) simulations, exploits the essential dynamics (ED) method, originally developed to analyse the collective motions in proteins. We apply this methodology to identify the residues that are more relevant for determining the reduction potential (E(0)) of a redox-active protein. To this aim, the fluctuation modes of the single-residue electrostatic potentials mostly contributing to the fluctuations of the total electrostatic potential (the main determinant of E(0)) are investigated for wild-type azurin and two of its mutants with a higher E(0). By comparing the results here obtained with a previous study on the same systems [Zanetti-Polzi et al., Org. Biomol. Chem., 2015, 13, 11003] we show that the proposed methodology is able to identify the key sites that determine E(0). This information can be used for a general deeper understanding of the molecular mechanisms on the basis of the redox properties of the proteins under investigation, as well as for the rational design of mutants with a higher or lower E(0). From the results of the present analysis we propose a new azurin mutant that, according to our calculations, shows a further increase of E(0).
Collapse
Affiliation(s)
- Laura Zanetti-Polzi
- Department of Physical and Chemical Sciences, University of L'Aquila, via Vetoio (Coppito 1), 67010, L'Aquila, Italy.
| | | | | | | |
Collapse
|
80
|
Samudio BM, Couch V, Stuchebrukhov AA. Monte Carlo Simulations of Glu-242 in Cytochrome c Oxidase. J Phys Chem B 2016; 120:2095-105. [PMID: 26865374 DOI: 10.1021/acs.jpcb.5b10998] [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/26/2023]
Abstract
Monte Carlo (MC) simulations of conformational changes and protonation of Glu-242, a key residue that shuttles protons in cytochrome c oxidase (CcO), are reported. Previous studies suggest that this residue may play a role of the valve of the enzyme proton pump. Here we examine how sensitive the results of simulations are to the computational method used. We applied both molecular mechanic (MM) and hybrid quantum mechanic:molecular mechanic (QM:MM) methods and find that the results are qualitatively different. The results indicate that the mechanism for proton gating in CcO is still an open issue.
Collapse
Affiliation(s)
- Benjamin M Samudio
- Department of Chemistry , University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| | - Vernon Couch
- Department of Chemistry , University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| | - Alexei A Stuchebrukhov
- Department of Chemistry , University of California, Davis, One Shields Avenue, Davis, California 95616, United States
| |
Collapse
|
81
|
Peterson C, Penchoff D, Wilson A. Prediction of Thermochemical Properties Across the Periodic Table. ANNUAL REPORTS IN COMPUTATIONAL CHEMISTRY 2016. [DOI: 10.1016/bs.arcc.2016.04.001] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
82
|
Zanetti-Polzi L, Corni S. A dynamical approach to non-adiabatic electron transfers at the bio-inorganic interface. Phys Chem Chem Phys 2016; 18:10538-49. [DOI: 10.1039/c6cp00044d] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
A methodology is proposed to investigate the role of the energy fluctuations, determined by the dynamical evolution of a system, and the role of non-adiabaticity in affecting the kinetic rate of electron transfer reactions at the bio-inorganic interface.
Collapse
|
83
|
Shift-and-invert parallel spectral transformation eigensolver: Massively parallel performance for density-functional based tight-binding. J Comput Chem 2015; 37:448-59. [DOI: 10.1002/jcc.24254] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2015] [Revised: 10/15/2015] [Accepted: 10/25/2015] [Indexed: 01/12/2023]
|
84
|
Ferguson GA, Cheng L, Bu L, Kim S, Robichaud DJ, Nimlos MR, Curtiss LA, Beckham GT. Carbocation Stability in H-ZSM5 at High Temperature. J Phys Chem A 2015; 119:11397-405. [PMID: 26501585 DOI: 10.1021/acs.jpca.5b07025] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Glen A. Ferguson
- National
Bioenergy Center, National Renewable Energy Laboratory, 15013 Denver
West Parkway, Golden, Colorado 80401, United States
| | - Lei Cheng
- Material
Science Division, Argonne National Laboratory, 9700 South Cass Avenue B109, Lemont, Illinois 60439, United States
| | - Lintao Bu
- National
Bioenergy Center, National Renewable Energy Laboratory, 15013 Denver
West Parkway, Golden, Colorado 80401, United States
| | - Seonah Kim
- National
Bioenergy Center, National Renewable Energy Laboratory, 15013 Denver
West Parkway, Golden, Colorado 80401, United States
| | - David J. Robichaud
- National
Bioenergy Center, National Renewable Energy Laboratory, 15013 Denver
West Parkway, Golden, Colorado 80401, United States
| | - Mark R. Nimlos
- National
Bioenergy Center, National Renewable Energy Laboratory, 15013 Denver
West Parkway, Golden, Colorado 80401, United States
| | - Larry A. Curtiss
- Material
Science Division, Argonne National Laboratory, 9700 South Cass Avenue B109, Lemont, Illinois 60439, United States
| | - Gregg T. Beckham
- National
Bioenergy Center, National Renewable Energy Laboratory, 15013 Denver
West Parkway, Golden, Colorado 80401, United States
| |
Collapse
|
85
|
Encapsulation of monomers, homodimers and heterodimers of amides and carboxylic acids in three non-covalent assemblies. Struct Chem 2015. [DOI: 10.1007/s11224-015-0682-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
86
|
Zanetti-Polzi L, Bortolotti CA, Daidone I, Aschi M, Amadei A, Corni S. A few key residues determine the high redox potential shift in azurin mutants. Org Biomol Chem 2015; 13:11003-13. [PMID: 26381463 DOI: 10.1039/c5ob01819f] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The wide range of variability of the reduction potential (E(0)) of blue-copper proteins has been the subject of a large number of studies in the past several years. In particular, a series of azurin mutants have been recently rationally designed tuning E(0) over a very broad range (700 mV) without significantly altering the redox-active site [Marshall et al., Nature, 2009, 462, 113]. This clearly suggests that interactions outside the primary coordination sphere are relevant to determine E(0) in cupredoxins. However, the molecular determinants of the redox potential variability are still undisclosed. Here, by means of atomistic molecular dynamics simulations and hybrid quantum/classical calculations, the mechanisms that determine the E(0) shift of two azurin mutants with high potential shifts are unravelled. The reduction potentials of native azurin and of the mutants are calculated obtaining results in good agreement with the experiments. The analysis of the simulations reveals that only a small number of residues (including non-mutated ones) are relevant in determining the experimentally observed E(0) variation via site-specific, but diverse, mechanisms. These findings open the path to the rational design of new azurin mutants with different E(0).
Collapse
Affiliation(s)
- Laura Zanetti-Polzi
- Center S3, CNR-Institute of Nanoscience, Via Campi 213/A, 41125, Modena, Italy.
| | | | | | | | | | | |
Collapse
|
87
|
Jovan Jose K, Raghavachari K. Molecules-in-molecules fragment-based method for the evaluation of Raman spectra of large molecules. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1074744] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
88
|
Revealing the active role of the terminal CNO moiety in the photochemical oxaziridine conversion process of some chemopreventive retinylnitrones through hybrid QM:QM and QM:MM ONIOM calculations. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2015.05.005] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
89
|
Tzeli D, Petsalakis ID, Theodorakopoulos G, Rebek J. Reversible encapsulation in a covalent capsule. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2015.05.018] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
90
|
Byrd JN, Jindal N, Molt RW, Bartlett RJ, Sanders BA, Lotrich VF. Molecular cluster perturbation theory. I. Formalism. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1036145] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
91
|
Chung LW, Sameera WMC, Ramozzi R, Page AJ, Hatanaka M, Petrova GP, Harris TV, Li X, Ke Z, Liu F, Li HB, Ding L, Morokuma K. The ONIOM Method and Its Applications. Chem Rev 2015; 115:5678-796. [PMID: 25853797 DOI: 10.1021/cr5004419] [Citation(s) in RCA: 738] [Impact Index Per Article: 82.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Lung Wa Chung
- †Department of Chemistry, South University of Science and Technology of China, Shenzhen 518055, China
| | - W M C Sameera
- ‡Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo, Kyoto 606-8103, Japan
| | - Romain Ramozzi
- ‡Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo, Kyoto 606-8103, Japan
| | - Alister J Page
- §Newcastle Institute for Energy and Resources, The University of Newcastle, Callaghan 2308, Australia
| | - Miho Hatanaka
- ‡Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo, Kyoto 606-8103, Japan
| | - Galina P Petrova
- ∥Faculty of Chemistry and Pharmacy, University of Sofia, Bulgaria Boulevard James Bourchier 1, 1164 Sofia, Bulgaria
| | - Travis V Harris
- ‡Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo, Kyoto 606-8103, Japan.,⊥Department of Chemistry, State University of New York at Oswego, Oswego, New York 13126, United States
| | - Xin Li
- #State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Zhuofeng Ke
- ∇School of Chemistry and Chemical Engineering, Sun Yat-sen University, Guangzhou 510275, China
| | - Fengyi Liu
- ○Key Laboratory of Macromolecular Science of Shaanxi Province, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an, Shaanxi 710119, China
| | - Hai-Bei Li
- ■School of Ocean, Shandong University, Weihai 264209, China
| | - Lina Ding
- ▲School of Pharmaceutical Sciences, Zhengzhou University, 100 Kexue Avenue, Zhengzhou, Henan 450001, China
| | - Keiji Morokuma
- ‡Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo, Kyoto 606-8103, Japan
| |
Collapse
|
92
|
Jose KVJ, Raghavachari K. Evaluation of Energy Gradients and Infrared Vibrational Spectra through Molecules-in-Molecules Fragment-Based Approach. J Chem Theory Comput 2015; 11:950-61. [DOI: 10.1021/ct501026m] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- K. V. Jovan Jose
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
93
|
Saini P, Chattopadhyay A. A comprehensive spectroscopic investigation of α-(2-naphthyl)-N-methylnitrone: a computational study on photochemical nitrone–oxaziridine conversion and thermal E–Z isomerization processes. RSC Adv 2015. [DOI: 10.1039/c4ra16375c] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Abstract
CASSCF and 2-layer hybrid ONIOM-based computational studies on α-(2-naphthyl)-N-methylnitrone have proposed its photochemical oxaziridine formation and thermal E–Z isomerization mechanisms.
Collapse
Affiliation(s)
- Praveen Saini
- Department of Chemistry
- Birla Institute of Technology and Science (BITS)
- Goa
- India
| | - Anjan Chattopadhyay
- Department of Chemistry
- Birla Institute of Technology and Science (BITS)
- Goa
- India
| |
Collapse
|
94
|
Bera A, Maroo S, Bhattacharya A. Electronically nonadiabatic decomposition mechanisms of clusters of zinc and dimethylnitramine. Chem Phys 2015. [DOI: 10.1016/j.chemphys.2014.11.003] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
95
|
Jurinovich S, Viani L, Curutchet C, Mennucci B. Limits and potentials of quantum chemical methods in modelling photosynthetic antennae. Phys Chem Chem Phys 2015; 17:30783-92. [DOI: 10.1039/c5cp00986c] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
A critical overview of quantum chemical approaches to simulate the light-harvesting process in photosynthetic antennae is presented together with a perspective on the developments that need to be introduced to reach a quantitative predictive power.
Collapse
Affiliation(s)
- Sandro Jurinovich
- Dipartimento di Chimica e Chimica Industriale
- University of Pisa
- 56124 Pisa
- Italy
| | - Lucas Viani
- Dipartimento di Chimica e Chimica Industriale
- University of Pisa
- 56124 Pisa
- Italy
- Institute for Fluid Dynamics
| | - Carles Curutchet
- Departament de Fisicoquímica
- Facultat de Farmàcia
- Universitat de Barcelona
- 08028 Barcelona
- Spain
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale
- University of Pisa
- 56124 Pisa
- Italy
| |
Collapse
|
96
|
Acevedo O. Simulating chemical reactions in ionic liquids using QM/MM methodology. J Phys Chem A 2014; 118:11653-66. [PMID: 25329366 DOI: 10.1021/jp507967z] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
The use of ionic liquids as a reaction medium for chemical reactions has dramatically increased in recent years due in large part to the numerous reported advances in catalysis and organic synthesis. In some extreme cases, ionic liquids have been shown to induce mechanistic changes relative to conventional solvents. Despite the large interest in the solvents, a clear understanding of the molecular factors behind their chemical impact is largely unknown. This feature article reviews our efforts developing and applying mixed quantum and molecular mechanical (QM/MM) methodology to elucidate the microscopic details of how these solvents operate to enhance rates and alter mechanisms for industrially and academically important reactions, e.g., Diels-Alder, Kemp eliminations, nucleophilic aromatic substitutions, and β-eliminations. Explicit solvent representation provided the medium dependence of the activation barriers and atomic-level characterization of the solute-solvent interactions responsible for the experimentally observed "ionic liquid effects". Technical advances are also discussed, including a linear-scaling pairwise electrostatic interaction alternative to Ewald sums, an efficient polynomial fitting method for modeling proton transfers, and the development of a custom ionic liquid OPLS-AA force field.
Collapse
Affiliation(s)
- Orlando Acevedo
- Department of Chemistry and Biochemistry, Auburn University , Auburn, Alabama 36849, United States
| |
Collapse
|
97
|
Khavrutskii IV, Legler PM, Friedlander AM, Wallqvist A. A reaction path study of the catalysis and inhibition of the Bacillus anthracis CapD γ-glutamyl transpeptidase. Biochemistry 2014; 53:6954-67. [PMID: 25334088 DOI: 10.1021/bi500623c] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
The CapD enzyme of Bacillus anthracis is a γ-glutamyl transpeptidase from the N-terminal nucleophile hydrolase superfamily that covalently anchors the poly-γ-D-glutamic acid (pDGA) capsule to the peptidoglycan. The capsule hinders phagocytosis of B. anthracis by host cells and is essential for virulence. The role CapD plays in capsule anchoring and remodeling makes the enzyme a promising target for anthrax medical countermeasures. Although the structure of CapD is known, and a covalent inhibitor, capsidin, has been identified, the mechanisms of CapD catalysis and inhibition are poorly understood. Here, we used a computational approach to map out the reaction steps involved in CapD catalysis and inhibition. We found that the rate-limiting step of either CapD catalysis or inhibition was a concerted asynchronous formation of the tetrahedral intermediate with a barrier of 22-23 kcal/mol. However, the mechanisms of these reactions differed for the two amides. The formation of the tetrahedral intermediate with pDGA was substrate-assisted with two proton transfers. In contrast, capsidin formed the tetrahedral intermediate in a conventional way with one proton transfer. Interestingly, capsidin coupled a conformational change in the catalytic residue of the tetrahedral intermediate to stretching of the scissile amide bond. Furthermore, capsidin took advantage of iminol-amide tautomerism of its diacetamide moiety to convert the tetrahedral intermediate to the acetylated CapD. As evidence of the promiscuous nature of CapD, the enzyme cleaved the amide bond of capsidin by attacking it on the opposite side compared to pDGA.
Collapse
Affiliation(s)
- Ilja V Khavrutskii
- Department of Defense Biotechnology High Performance Computing Software Applications Institute, Telemedicine and Advanced Technology Research Center, United States Army Medical Research and Materiel Command , Fort Detrick, Maryland 21702, United States
| | | | | | | |
Collapse
|
98
|
Libisch F, Huang C, Carter EA. Embedded correlated wavefunction schemes: theory and applications. Acc Chem Res 2014; 47:2768-75. [PMID: 24873211 DOI: 10.1021/ar500086h] [Citation(s) in RCA: 149] [Impact Index Per Article: 14.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Conspectus Ab initio modeling of matter has become a pillar of chemical research: with ever-increasing computational power, simulations can be used to accurately predict, for example, chemical reaction rates, electronic and mechanical properties of materials, and dynamical properties of liquids. Many competing quantum mechanical methods have been developed over the years that vary in computational cost, accuracy, and scalability: density functional theory (DFT), the workhorse of solid-state electronic structure calculations, features a good compromise between accuracy and speed. However, approximate exchange-correlation functionals limit DFT's ability to treat certain phenomena or states of matter, such as charge-transfer processes or strongly correlated materials. Furthermore, conventional DFT is purely a ground-state theory: electronic excitations are beyond its scope. Excitations in molecules are routinely calculated using time-dependent DFT linear response; however applications to condensed matter are still limited. By contrast, many-electron wavefunction methods aim for a very accurate treatment of electronic exchange and correlation. Unfortunately, the associated computational cost renders treatment of more than a handful of heavy atoms challenging. On the other side of the accuracy spectrum, parametrized approaches like tight-binding can treat millions of atoms. In view of the different (dis-)advantages of each method, the simulation of complex systems seems to force a compromise: one is limited to the most accurate method that can still handle the problem size. For many interesting problems, however, compromise proves insufficient. A possible solution is to break up the system into manageable subsystems that may be treated by different computational methods. The interaction between subsystems may be handled by an embedding formalism. In this Account, we review embedded correlated wavefunction (CW) approaches and some applications. We first discuss our density functional embedding theory, which is formally exact. We show how to determine the embedding potential, which replaces the interaction between subsystems, at the DFT level. CW calculations are performed using a fixed embedding potential, that is, a non-self-consistent embedding scheme. We demonstrate this embedding theory for two challenging electron transfer phenomena: (1) initial oxidation of an aluminum surface and (2) hot-electron-mediated dissociation of hydrogen molecules on a gold surface. In both cases, the interaction between gas molecules and metal surfaces were treated by sophisticated CW techniques, with the remainder of the extended metal surface being treated by DFT. Our embedding approach overcomes the limitations of conventional Kohn-Sham DFT in describing charge transfer, multiconfigurational character, and excited states. From these embedding simulations, we gained important insights into fundamental processes that are crucial aspects of fuel cell catalysis (i.e., O2 reduction at metal surfaces) and plasmon-mediated photocatalysis by metal nanoparticles. Moreover, our findings agree very well with experimental observations, while offering new views into the chemistry. We finally discuss our recently formulated potential-functional embedding theory that provides a seamless, first-principles way to include back-action onto the environment from the embedded region.
Collapse
Affiliation(s)
- Florian Libisch
- Institute
for Theoretical Physics, Vienna University of Technology, 1040 Vienna, Austria
| | - Chen Huang
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Emily A. Carter
- Department
of Mechanical and Aerospace Engineering, Program in Applied and Computational
Mathematics, and Andlinger Center for Energy and the Environment, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
99
|
Zanetti-Polzi L, Daidone I, Bortolotti CA, Corni S. Surface Packing Determines the Redox Potential Shift of Cytochrome c Adsorbed on Gold. J Am Chem Soc 2014; 136:12929-37. [DOI: 10.1021/ja505251a] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Affiliation(s)
- Laura Zanetti-Polzi
- Center
S3, CNR NANO, Institute of Nanoscience, Via Campi 213/A, 41125, Modena, Italy
| | - Isabella Daidone
- Department
of Physical and Chemical Sciences, University of L’Aquila, via
Vetoio (Coppito 1), 67010, L’Aquila, Italy
| | - Carlo Augusto Bortolotti
- Center
S3, CNR NANO, Institute of Nanoscience, Via Campi 213/A, 41125, Modena, Italy
- Department
of Life Sciences, University of Modena and Reggio Emilia, Via Campi
183, 41125, Modena, Italy
| | - Stefano Corni
- Center
S3, CNR NANO, Institute of Nanoscience, Via Campi 213/A, 41125, Modena, Italy
| |
Collapse
|
100
|
|