51
|
Yang Z, Twidale RM, Gervasoni S, Suardíaz R, Colenso CK, Lang EJM, Spencer J, Mulholland AJ. Multiscale Workflow for Modeling Ligand Complexes of Zinc Metalloproteins. J Chem Inf Model 2021; 61:5658-5672. [PMID: 34748329 DOI: 10.1021/acs.jcim.1c01109] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Zinc metalloproteins are ubiquitous, with protein zinc centers of structural and functional importance, involved in interactions with ligands and substrates and often of pharmacological interest. Biomolecular simulations are increasingly prominent in investigations of protein structure, dynamics, ligand interactions, and catalysis, but zinc poses a particular challenge, in part because of its versatile, flexible coordination. A computational workflow generating reliable models of ligand complexes of biological zinc centers would find broad application. Here, we evaluate the ability of alternative treatments, using (nonbonded) molecular mechanics (MM) and quantum mechanics/molecular mechanics (QM/MM) at semiempirical (DFTB3) and density functional theory (DFT) levels of theory, to describe the zinc centers of ligand complexes of six metalloenzyme systems differing in coordination geometries, zinc stoichiometries (mono- and dinuclear), and the nature of interacting groups (specifically the presence of zinc-sulfur interactions). MM molecular dynamics (MD) simulations can overfavor octahedral geometries, introducing additional water molecules to the zinc coordination shell, but this can be rectified by subsequent semiempirical (DFTB3) QM/MM MD simulations. B3LYP/MM geometry optimization further improved the accuracy of the description of coordination distances, with the overall effectiveness of the approach depending upon factors, including the presence of zinc-sulfur interactions that are less well described by semiempirical methods. We describe a workflow comprising QM/MM MD using DFTB3 followed by QM/MM geometry optimization using DFT (e.g., B3LYP) that well describes our set of zinc metalloenzyme complexes and is likely to be suitable for creating accurate models of zinc protein complexes when structural information is more limited.
Collapse
Affiliation(s)
- Zongfan Yang
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TH, U.K.,School of Cellular and Molecular Medicine, University of Bristol, Bristol BS8 1TD, U.K
| | - Rebecca M Twidale
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TH, U.K
| | - Silvia Gervasoni
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TH, U.K.,Department of Pharmaceutical Sciences, University of Milan, Via Mangiagalli, 25, I-20133 Milano, Italy
| | - Reynier Suardíaz
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TH, U.K
| | - Charlotte K Colenso
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TH, U.K.,School of Cellular and Molecular Medicine, University of Bristol, Bristol BS8 1TD, U.K
| | - Eric J M Lang
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TH, U.K
| | - James Spencer
- School of Cellular and Molecular Medicine, University of Bristol, Bristol BS8 1TD, U.K
| | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TH, U.K
| |
Collapse
|
52
|
Lu SY, Zuehlsdorff TJ, Hong H, Aguirre VP, Isborn CM, Shi L. The Influence of Electronic Polarization on Nonlinear Optical Spectroscopy. J Phys Chem B 2021; 125:12214-12227. [PMID: 34726915 DOI: 10.1021/acs.jpcb.1c05914] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The environment surrounding a chromophore can dramatically affect the energy absorption and relaxation process, as manifested in optical spectra. Simulations of nonlinear optical spectroscopy, such as two-dimensional electronic spectroscopy (2DES) and transient absorption (TA), will be influenced by the computational model of the environment. We here compare a fixed point charge molecular mechanics model and a quantum mechanical (QM) model of the environment in computed 2DES and TA spectra of Nile red in water and the chromophore of photoactive yellow protein (PYP) in water and protein environments. In addition to simulating these nonlinear optical spectra, we directly juxtapose the computed excitation energy correlation function to the dynamic Stokes shift function often used to analyze environment dynamics. Overall, we find that for the three systems studied here the mutual electronic polarization provided by the QM environment manifests in broader 2DES signals, as well as a larger reorganization energy and a larger static Stokes shift due to stronger coupling between the chromophore and the environment.
Collapse
Affiliation(s)
- Shao-Yu Lu
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, United States
| | - Hanbo Hong
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Vincent P Aguirre
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Christine M Isborn
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Liang Shi
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
53
|
Modeling Catalysis in Allosteric Enzymes: Capturing Conformational Consequences. Top Catal 2021; 65:165-186. [DOI: 10.1007/s11244-021-01521-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
54
|
Berger MB, Walker AR, Vázquez-Montelongo EA, Cisneros GA. Computational investigations of selected enzymes from two iron and α-ketoglutarate-dependent families. Phys Chem Chem Phys 2021; 23:22227-22240. [PMID: 34586107 PMCID: PMC8516722 DOI: 10.1039/d1cp03800a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
DNA alkylation is used as the key epigenetic mark in eukaryotes, however, most alkylation in DNA can result in deleterious effects. Therefore, this process needs to be tightly regulated. The enzymes of the AlkB and Ten-Eleven Translocation (TET) families are members of the Fe and alpha-ketoglutarate-dependent superfamily of enzymes that are tasked with dealkylating DNA and RNA in cells. Members of these families span all species and are an integral part of transcriptional regulation. While both families catalyze oxidative dealkylation of various bases, each has specific preference for alkylated base type as well as distinct catalytic mechanisms. This perspective aims to provide an overview of computational work carried out to investigate several members of these enzyme families including AlkB, ALKB Homolog 2, ALKB Homolog 3 and Ten-Eleven Translocate 2. Insights into structural details, mutagenesis studies, reaction path analysis, electronic structure features in the active site, and substrate preferences are presented and discussed.
Collapse
Affiliation(s)
- Madison B Berger
- Department of Chemistry, University of North Texas, Denton, Texas, 76201, USA.
| | - Alice R Walker
- Department of Chemistry, Wayne State University, Detroit, Michigan, 48202, USA
| | | | - G Andrés Cisneros
- Department of Chemistry, University of North Texas, Denton, Texas, 76201, USA.
| |
Collapse
|
55
|
Kang H, Zheng M. Influence of the quantum mechanical region size in QM/MM modelling: A case study of fluoroacetate dehalogenase catalyzed C F bond cleavage. COMPUT THEOR CHEM 2021. [DOI: 10.1016/j.comptc.2021.113399] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
56
|
Prioli S, Kongsted J. Modeling environmental effects in two-photon circular dichroism calculations. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02838-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
57
|
Unke O, Chmiela S, Sauceda HE, Gastegger M, Poltavsky I, Schütt KT, Tkatchenko A, Müller KR. Machine Learning Force Fields. Chem Rev 2021; 121:10142-10186. [PMID: 33705118 PMCID: PMC8391964 DOI: 10.1021/acs.chemrev.0c01111] [Citation(s) in RCA: 397] [Impact Index Per Article: 132.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Indexed: 12/27/2022]
Abstract
In recent years, the use of machine learning (ML) in computational chemistry has enabled numerous advances previously out of reach due to the computational complexity of traditional electronic-structure methods. One of the most promising applications is the construction of ML-based force fields (FFs), with the aim to narrow the gap between the accuracy of ab initio methods and the efficiency of classical FFs. The key idea is to learn the statistical relation between chemical structure and potential energy without relying on a preconceived notion of fixed chemical bonds or knowledge about the relevant interactions. Such universal ML approximations are in principle only limited by the quality and quantity of the reference data used to train them. This review gives an overview of applications of ML-FFs and the chemical insights that can be obtained from them. The core concepts underlying ML-FFs are described in detail, and a step-by-step guide for constructing and testing them from scratch is given. The text concludes with a discussion of the challenges that remain to be overcome by the next generation of ML-FFs.
Collapse
Affiliation(s)
- Oliver
T. Unke
- Machine
Learning Group, Technische Universität
Berlin, 10587 Berlin, Germany
- DFG
Cluster of Excellence “Unifying Systems in Catalysis”
(UniSysCat), Technische Universität Berlin, 10623 Berlin, Germany
| | - Stefan Chmiela
- Machine
Learning Group, Technische Universität
Berlin, 10587 Berlin, Germany
| | - Huziel E. Sauceda
- Machine
Learning Group, Technische Universität
Berlin, 10587 Berlin, Germany
- BASLEARN,
BASF-TU Joint Lab, Technische Universität
Berlin, 10587 Berlin, Germany
| | - Michael Gastegger
- Machine
Learning Group, Technische Universität
Berlin, 10587 Berlin, Germany
- DFG
Cluster of Excellence “Unifying Systems in Catalysis”
(UniSysCat), Technische Universität Berlin, 10623 Berlin, Germany
- BASLEARN,
BASF-TU Joint Lab, Technische Universität
Berlin, 10587 Berlin, Germany
| | - Igor Poltavsky
- Department
of Physics and Materials Science, University
of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Kristof T. Schütt
- Machine
Learning Group, Technische Universität
Berlin, 10587 Berlin, Germany
| | - Alexandre Tkatchenko
- Department
of Physics and Materials Science, University
of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Klaus-Robert Müller
- Machine
Learning Group, Technische Universität
Berlin, 10587 Berlin, Germany
- BIFOLD−Berlin
Institute for the Foundations of Learning and Data, Berlin, Germany
- Department
of Artificial Intelligence, Korea University, Anam-dong, Seongbuk-gu, Seoul 02841, Korea
- Max Planck
Institute for Informatics, Stuhlsatzenhausweg, 66123 Saarbrücken, Germany
- Google
Research, Brain Team, Berlin, Germany
| |
Collapse
|
58
|
Chen J, Kato J, Harper JB, Shao Y, Ho J. On the Accuracy of QM/MM Models: A Systematic Study of Intramolecular Proton Transfer Reactions of Amino Acids in Water. J Phys Chem B 2021; 125:9304-9316. [PMID: 34355564 DOI: 10.1021/acs.jpcb.1c04876] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
This work presents a systematic assessment of QM/QM' and QM/MM models with respect to direct QM calculations for the tautomerization (neutral to zwitterion) reactions of amino acids (glycine, alanine, valine, aspartate, and neutral and protonated histidine) solvated in a 160 water cluster. The effect of varying QM region size and choice of embedding potentials, including fixed-charge and polarizable molecular mechanics force fields (TIP3P and EFP) and various semiempirical QM methods (PM7, GFN2-xTB, DFTBA, DFTB3, HF-3c, and PBEh-3c), on the accuracy of the models was examined. A surprising finding was that molecular mechanics force fields outperformed many of the semiempirical methods. Generally, the errors in the QM/QM' and QM/MM models converge slowly with respect to the QM region size, requiring 50 or more waters to be included in the QM region before the error in the model falls below 1 kcal mol-1 of its pure QM result. Different QM region selection schemes were also compared, and it was found that selection based on Natural Population Analysis (NPA) atomic charges significantly reduced the error in the QM/QM' and QM/MM models particularly if a low-quality embedding potential was used. It is envisaged that these results will be useful for the development of future hybrid QM models.
Collapse
Affiliation(s)
- Junbo Chen
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| | - Jin Kato
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| | - Jason B Harper
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Junming Ho
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| |
Collapse
|
59
|
Yagi K, Sugita Y. Anharmonic Vibrational Calculations Based on Group-Localized Coordinates: Applications to Internal Water Molecules in Bacteriorhodopsin. J Chem Theory Comput 2021; 17:5007-5020. [PMID: 34296615 PMCID: PMC10986902 DOI: 10.1021/acs.jctc.1c00060] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
An efficient anharmonic vibrational method is developed exploiting the locality of molecular vibration. Vibrational coordinates localized to a group of atoms are employed to divide the potential energy surface (PES) of a system into intra- and inter-group contributions. Then, the vibrational Schrödinger equation is solved based on a PES, in which the inter-group coupling is truncated at the harmonic level while accounting for the intra-group anharmonicity. The method is applied to a pentagonal hydrogen bond network (HBN) composed of internal water molecules and charged residues in a membrane protein, bacteriorhodopsin. The PES is calculated by the quantum mechanics/molecular mechanics (QM/MM) calculation at the level of B3LYP-D3/aug-cc-pVDZ. The infrared (IR) spectrum is computed using a set of coordinates localized to each water molecule and amino acid residue by second-order vibrational quasi-degenerate perturbation theory (VQDPT2). Benchmark calculations show that the proposed method yields the N-D/O-D stretching frequencies with an error of 7 cm-1 at the cost reduced by more than five times. In contrast, the harmonic approximation results in a severe error of 150 cm-1. Furthermore, the size of QM regions is carefully assessed to find that the QM regions should include not only the pentagonal HBN itself but also its HB partners. VQDPT2 calculations starting from transient structures obtained by molecular dynamics simulations have shown that the structural sampling has a significant impact on the calculated IR spectrum. The incorporation of anharmonicity, sufficiently large QM regions, and structural samplings are of essential importance to reproduce the experimental IR spectrum. The computational spectrum paves the way for decoding the IR signal of strong HBNs and helps elucidate their functional roles in biomolecules.
Collapse
Affiliation(s)
- Kiyoshi Yagi
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Yuji Sugita
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, 1-6-5 Minatojima-Minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
60
|
Harder, better, faster, stronger: Large-scale QM and QM/MM for predictive modeling in enzymes and proteins. Curr Opin Struct Biol 2021; 72:9-17. [PMID: 34388673 DOI: 10.1016/j.sbi.2021.07.004] [Citation(s) in RCA: 36] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2021] [Revised: 06/25/2021] [Accepted: 07/05/2021] [Indexed: 11/23/2022]
Abstract
Computational prediction of enzyme mechanism and protein function requires accurate physics-based models and suitable sampling. We discuss recent advances in large-scale quantum mechanical (QM) modeling of biochemical systems that have reduced the cost of high-accuracy models. Tradeoffs between sampling and accuracy have motivated modeling with molecular mechanics (MM) in a multiscale QM/MM or iterative approach. Limitations to both conventional density-functional theory and classical MM force fields remain for describing noncovalent interactions in comparison to experiment or wavefunction theory. Because predictions of enzyme action (i.e. electrostatics), free energy barriers, and mechanisms are sensitive to the protocol and embedding method in QM/MM, convergence tests and systematic methods for quantifying QM-level interactions are a needed, active area of development.
Collapse
|
61
|
Summers TJ, Cheng Q, Palma MA, Pham DT, Kelso DK, Webster CE, DeYonker NJ. Cheminformatic quantum mechanical enzyme model design: A catechol-O-methyltransferase case study. Biophys J 2021; 120:3577-3587. [PMID: 34358526 DOI: 10.1016/j.bpj.2021.07.029] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2021] [Revised: 05/26/2021] [Accepted: 07/29/2021] [Indexed: 10/20/2022] Open
Abstract
To accurately simulate the inner workings of an enzyme active site with quantum mechanics (QM), not only must the reactive species be included in the model but also important surrounding residues, solvent, or coenzymes involved in crafting the microenvironment. Our lab has been developing the Residue Interaction Network Residue Selector (RINRUS) toolkit to utilize interatomic contact network information for automated, rational residue selection and QM-cluster model generation. Starting from an x-ray crystal structure of catechol-O-methyltransferase, RINRUS was used to construct a series of QM-cluster models. The reactant, product, and transition state of the methyl transfer reaction were computed for a total of 550 models, and the resulting free energies of activation and reaction were used to evaluate model convergence. RINRUS-designed models with only 200-300 atoms are shown to converge. RINRUS will serve as a cornerstone for improved and automated cheminformatics-based enzyme model design.
Collapse
Affiliation(s)
- Thomas J Summers
- Department of Chemistry, The University of Memphis, Memphis, Tennessee
| | - Qianyi Cheng
- Department of Chemistry, The University of Memphis, Memphis, Tennessee
| | - Manuel A Palma
- Department of Chemistry, The University of Memphis, Memphis, Tennessee
| | - Diem-Trang Pham
- Department of Chemistry, The University of Memphis, Memphis, Tennessee; Department of Computer Science, The University of Memphis, Memphis, Tennessee
| | - Dudley K Kelso
- Department of Chemistry, The University of Memphis, Memphis, Tennessee
| | - Charles Edwin Webster
- Department of Chemistry, Mississippi State University, Mississippi State, Mississippi
| | - Nathan J DeYonker
- Department of Chemistry, The University of Memphis, Memphis, Tennessee.
| |
Collapse
|
62
|
Hix MA, Leddin EM, Cisneros GA. Combining Evolutionary Conservation and Quantum Topological Analyses To Determine Quantum Mechanics Subsystems for Biomolecular Quantum Mechanics/Molecular Mechanics Simulations. J Chem Theory Comput 2021; 17:4524-4537. [PMID: 34087064 PMCID: PMC8477969 DOI: 10.1021/acs.jctc.1c00313] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Selection of residues and other molecular fragments for inclusion in the quantum mechanics (QM) region for QM/molecular mechanics (MM) simulations is an important step for these calculations. Here, we present an approach that combines protein sequence/structure evolution and electron localization function (ELF) analyses. The combination of these two analyses allows the determination of whether a residue needs to be included in the QM subsystem or can be represented by the MM environment. We have applied this approach on two systems previously investigated by QM/MM simulations, 4-oxalocrotonate tautomerase (4OT) and ten-eleven translocation-2 (TET2), that provide examples where fragments may or may not need to be included in the QM subsystem. Subsequently, we present the use of this approach to determine the appropriate QM subsystem to calculate the minimum energy path (MEP) for the reaction catalyzed by human DNA polymerase λ (Polλ) with a third cation in the active site. Our results suggest that the combination of protein evolutionary and ELF analyses provides insights into residue/molecular fragment selection for QM/MM simulations.
Collapse
Affiliation(s)
- Mark A Hix
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - Emmett M Leddin
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - G Andrés Cisneros
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| |
Collapse
|
63
|
Yang Z, Kulik HJ. Protein Dynamics and Substrate Protonation States Mediate the Catalytic Action of trans-4-Hydroxy-l-Proline Dehydratase. J Phys Chem B 2021; 125:7774-7784. [PMID: 34236200 DOI: 10.1021/acs.jpcb.1c05320] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The enzyme trans-4-hydroxy-l-proline (Hyp) dehydratase (HypD) is among the most abundant glycyl radical enzymes (GREs) in the healthy human gut microbiome and is considered a promising antibiotic target for the prominent antibiotic-resistant pathogen Clostridium difficile. Although an enzymatic mechanism has been proposed, the role of the greater HypD protein environment in mediating radical reactivity is not well understood. To fill this gap in understanding, we investigate HypD across multiple time- and length-scales using electronic structure modeling and classical molecular dynamics. We observe that the Hyp substrate protonation state significantly alters both its enzyme-free reactivity and its dynamics within the enzyme active site. Accurate coupled-cluster modeling suggests the deprotonated form of Hyp to be the most reactive protonation state for C5-Hpro-S activation. In the protein environment, hydrophobic interactions modulate the positioning of the Cys434 radical to enhance the reactivity of C5-Hpro-S abstraction. Long-time dynamics reveal that changing Hyp protonation states triggers the switching of a Leu643-gated water tunnel, a functional feature that has not yet been observed for members of the GRE superfamily.
Collapse
Affiliation(s)
- Zhongyue Yang
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
64
|
Abstract
Machine learning (ML) techniques applied to chemical reactions have a long history. The present contribution discusses applications ranging from small molecule reaction dynamics to computational platforms for reaction planning. ML-based techniques can be particularly relevant for problems involving both computation and experiments. For one, Bayesian inference is a powerful approach to develop models consistent with knowledge from experiments. Second, ML-based methods can also be used to handle problems that are formally intractable using conventional approaches, such as exhaustive characterization of state-to-state information in reactive collisions. Finally, the explicit simulation of reactive networks as they occur in combustion has become possible using machine-learned neural network potentials. This review provides an overview of the questions that can and have been addressed using machine learning techniques, and an outlook discusses challenges in this diverse and stimulating field. It is concluded that ML applied to chemistry problems as practiced and conceived today has the potential to transform the way with which the field approaches problems involving chemical reactions, in both research and academic teaching.
Collapse
Affiliation(s)
- Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland.,Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
65
|
Brunken C, Reiher M. Automated Construction of Quantum–Classical Hybrid Models. J Chem Theory Comput 2021; 17:3797-3813. [DOI: 10.1021/acs.jctc.1c00178] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Christoph Brunken
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
66
|
Zuehlsdorff TJ, Shedge SV, Lu SY, Hong H, Aguirre VP, Shi L, Isborn CM. Vibronic and Environmental Effects in Simulations of Optical Spectroscopy. Annu Rev Phys Chem 2021; 72:165-188. [DOI: 10.1146/annurev-physchem-090419-051350] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Including both environmental and vibronic effects is important for accurate simulation of optical spectra, but combining these effects remains computationally challenging. We outline two approaches that consider both the explicit atomistic environment and the vibronic transitions. Both phenomena are responsible for spectral shapes in linear spectroscopy and the electronic evolution measured in nonlinear spectroscopy. The first approach utilizes snapshots of chromophore-environment configurations for which chromophore normal modes are determined. We outline various approximations for this static approach that assumes harmonic potentials and ignores dynamic system-environment coupling. The second approach obtains excitation energies for a series of time-correlated snapshots. This dynamic approach relies on the accurate truncation of the cumulant expansion but treats the dynamics of the chromophore and the environment on equal footing. Both approaches show significant potential for making strides toward more accurate optical spectroscopy simulations of complex condensed phase systems.
Collapse
Affiliation(s)
- Tim J. Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | - Sapana V. Shedge
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Shao-Yu Lu
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Hanbo Hong
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Vincent P. Aguirre
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Liang Shi
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Christine M. Isborn
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| |
Collapse
|
67
|
Wang JN, Liu W, Li P, Mo Y, Hu W, Zheng J, Pan X, Shao Y, Mei Y. Accelerated Computation of Free Energy Profile at Ab Initio Quantum Mechanical/Molecular Mechanics Accuracy via a Semiempirical Reference Potential. 4. Adaptive QM/MM. J Chem Theory Comput 2021; 17:1318-1325. [PMID: 33593057 PMCID: PMC8335528 DOI: 10.1021/acs.jctc.0c01149] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Although quantum mechanical/molecular mechanics (QM/MM) methods are now routinely applied to the studies of chemical reactions in condensed phases and enzymatic reactions, they may experience technical difficulties when the reactive region is varying over time. For instance, when the solvent molecules are directly participating in the reaction, the exchange of water molecules between the QM and MM regions may occur on a time scale comparable to the reaction time. To cope with this situation, several adaptive QM/MM schemes have been proposed. However, these methods either add significantly to the computational cost or introduce artificial restraints to the system. In this work, we developed a novel adaptive QM/MM scheme and applied it to the study of a nucleophilic addition reaction. In this scheme, the configuration sampling was performed with a small QM region (without solvent molecules), and the thermodynamic properties under another potential energy function with a larger QM region (with a certain number of solvent molecules and/or different levels of QM theory) are computed via extrapolation using the reference-potential method. Our simulation results show that this adaptive QM/MM scheme is numerically stable, at least for the case studied in this work. Furthermore, this method also offers an inexpensive way to examine the convergence of the QM/MM calculation with respect to the size of the QM region.
Collapse
Affiliation(s)
- Jia-Ning Wang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Wei Liu
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Pengfei Li
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Yan Mo
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| | - Wenxin Hu
- The Computer Center, School of Data Science & Engineering, East China Normal University, Shanghai 200062, China
| | - Jun Zheng
- The Computer Center, School of Data Science & Engineering, East China Normal University, Shanghai 200062, China
| | - Xiaoliang Pan
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - 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 Electronic Science, East China Normal University, Shanghai 200062, China
- 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
|
68
|
Voice AT, Tresadern G, Twidale RM, van Vlijmen H, Mulholland AJ. Mechanism of covalent binding of ibrutinib to Bruton's tyrosine kinase revealed by QM/MM calculations. Chem Sci 2021; 12:5511-5516. [PMID: 33995994 PMCID: PMC8097726 DOI: 10.1039/d0sc06122k] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Ibrutinib is the first covalent inhibitor of Bruton's tyrosine kinase (BTK) to be used in the treatment of B-cell cancers. Understanding the mechanism of covalent inhibition will aid in the design of safer and more selective covalent inhibitors that target BTK. The mechanism of covalent inhibition in BTK has been uncertain because there is no appropriate residue nearby that can act as a base to deprotonate the cysteine thiol prior to covalent bond formation. We investigate several mechanisms of covalent modification of C481 in BTK by ibrutinib using combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics reaction simulations. The lowest energy pathway involves direct proton transfer from C481 to the acrylamide warhead in ibrutinib, followed by covalent bond formation to form an enol intermediate. There is a subsequent rate-limiting keto-enol tautomerisation step (ΔG ‡ = 10.5 kcal mol-1) to reach the inactivated BTK/ibrutinib complex. Our results represent the first mechanistic study of BTK inactivation by ibrutinib to consider multiple mechanistic pathways. These findings should aid in the design of covalent drugs that target BTK and other similar targets.
Collapse
Affiliation(s)
- Angus T Voice
- Centre for Computational Chemistry, School of Chemistry, University of Bristol Cantock's Close Bristol BS8 1TS UK
| | - Gary Tresadern
- Computational Chemistry, Janssen Research & Development, Janssen Pharmaceutica N. V. Turnhoutseweg 30 B-2340 Beerse Belgium
| | - Rebecca M Twidale
- Centre for Computational Chemistry, School of Chemistry, University of Bristol Cantock's Close Bristol BS8 1TS UK
| | - Herman van Vlijmen
- Computational Chemistry, Janssen Research & Development, Janssen Pharmaceutica N. V. Turnhoutseweg 30 B-2340 Beerse Belgium
| | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol Cantock's Close Bristol BS8 1TS UK
| |
Collapse
|
69
|
Abstract
QM/MM simulations have become an indispensable tool in many chemical and biochemical investigations. Considering the tremendous degree of success, including recognition by a 2013 Nobel Prize in Chemistry, are there still "burning challenges" in QM/MM methods, especially for biomolecular systems? In this short Perspective, we discuss several issues that we believe greatly impact the robustness and quantitative applicability of QM/MM simulations to many, if not all, biomolecules. We highlight these issues with observations and relevant advances from recent studies in our group and others in the field. Despite such limited scope, we hope the discussions are of general interest and will stimulate additional developments that help push the field forward in meaningful directions.
Collapse
Affiliation(s)
- Qiang Cui
- Departments of Chemistry, Physics, and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Tanmoy Pal
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Luke Xie
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
70
|
Maldonado AM, Hagiwara S, Choi TH, Eckert F, Schwarz K, Sundararaman R, Otani M, Keith JA. Quantifying Uncertainties in Solvation Procedures for Modeling Aqueous Phase Reaction Mechanisms. J Phys Chem A 2021; 125:154-164. [PMID: 33393781 DOI: 10.1021/acs.jpca.0c08961] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Computational quantum chemistry provides fundamental chemical and physical insights into solvated reaction mechanisms across many areas of chemistry, especially in homogeneous and heterogeneous renewable energy catalysis. Such reactions may depend on explicit interactions with ions and solvent molecules that are nontrivial to characterize. Rigorously modeling explicit solvent effects with molecular dynamics usually brings steep computational costs while the performance of continuum solvent models such as polarizable continuum model (PCM), charge-asymmetric nonlocally determined local-electric (CANDLE), conductor-like screening model for real solvents (COSMO-RS), and effective screening medium method with the reference interaction site model (ESM-RISM) are less well understood for reaction mechanisms. Here, we revisit a fundamental aqueous hydride transfer reaction-carbon dioxide (CO2) reduction by sodium borohydride (NaBH4)-as a test case to evaluate how different solvent models perform in aqueous phase charge migrations that would be relevant to renewable energy catalysis mechanisms. For this system, quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations almost exactly reproduced energy profiles from QM simulations, and the Na+ counterion in the QM/MM simulations plays an insignificant role over ensemble averaged trajectories that describe the reaction pathway. However, solvent models used on static calculations gave much more variability in data depending on whether the system was modeled using explicit solvent shells and/or the counterion. We pinpoint this variability due to unphysical descriptions of charge-separated states in the gas phase (i.e., self-interaction errors), and we show that using more accurate hybrid functionals and/or explicit solvent shells lessens these errors. This work closes with recommended procedures for treating solvation in future computational efforts in studying renewable energy catalysis mechanisms.
Collapse
Affiliation(s)
- Alex M Maldonado
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Satoshi Hagiwara
- National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba 305-8568, Japan
| | - Tae Hoon Choi
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Frank Eckert
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Kathleen Schwarz
- Material Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Ravishankar Sundararaman
- Department of Materials Science and Engineering, Rensselaer Polytechnic Institute, Troy, New York 12180, United States
| | - Minoru Otani
- National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba 305-8568, Japan
| | - John A Keith
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| |
Collapse
|
71
|
Gheorghiu A, Coveney PV, Arabi AA. The influence of base pair tautomerism on single point mutations in aqueous DNA. Interface Focus 2020; 10:20190120. [PMID: 33178413 PMCID: PMC7653342 DOI: 10.1098/rsfs.2019.0120] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/24/2020] [Indexed: 12/11/2022] Open
Abstract
The relationship between base pair hydrogen bond proton transfer and the rate of spontaneous single point mutations at ambient temperatures and pressures in aqueous DNA is investigated. By using an ensemble-based multiscale computational modelling method, statistically robust rates of proton transfer for the A:T and G:C base pairs within a solvated DNA dodecamer are calculated. Several different proton transfer pathways are observed within the same base pair. It is shown that, in G:C, the double proton transfer tautomer is preferred, while the single proton transfer process is favoured in A:T. The reported range of rate coefficients for double proton transfer is consistent with recent experimental data. Notwithstanding the approximately 1000 times more common presence of single proton transfer products from A:T, observationally there is bias towards G:C to A:T mutations in a wide range of living organisms. We infer that the double proton transfer reactions between G:C base pairs have a negligible contribution towards this bias for the following reasons: (i) the maximum half-life of the G*:C* tautomer is in the range of picoseconds, which is significantly smaller than the milliseconds it takes for DNA to unwind during replication, (ii) statistically, the majority of G*:C* tautomers revert back to their canonical forms through a barrierless process, and (iii) the thermodynamic instability of the tautomers with respect to the canonical base pairs. Through similar reasoning, we also deduce that proton transfer in the A:T base pair does not contribute to single point mutations in DNA.
Collapse
Affiliation(s)
- A Gheorghiu
- Centre for Computational Science, University College London, London, UK
| | - P V Coveney
- Centre for Computational Science, University College London, London, UK.,Informatics Institute, University of Amsterdam, Amsterdam, The Netherlands
| | - A A Arabi
- Centre for Computational Science, University College London, London, UK.,College of Medicine and Health Sciences, Biochemistry Department, United Arab Emirates University, PO Box 17666, Al Ain, United Arab Emirates
| |
Collapse
|
72
|
Raucci U, Chiariello MG, Rega N. Modeling Excited-State Proton Transfer to Solvent: A Dynamics Study of a Super Photoacid with a Hybrid Implicit/Explicit Solvent Model. J Chem Theory Comput 2020; 16:7033-7043. [PMID: 33112132 PMCID: PMC8016186 DOI: 10.1021/acs.jctc.0c00782] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
![]()
The rapid growth of time-resolved
spectroscopies and the theoretical
advances in ab initio molecular dynamics (AIMD) pave the way to look
at the real-time molecular motion following the electronic excitation.
Here, we exploited the capabilities of AIMD combined with a hybrid
implicit/explicit model of solvation to investigate the ultrafast
excited-state proton transfer (ESPT) reaction of a super photoacid,
known as QCy9, in water solution. QCy9 transfers a proton to a water
solvent molecule within 100 fs upon the electronic excitation in aqueous
solution, and it is the strongest photoacid reported in the literature
so far. Because of the ultrafast kinetics, it has been experimentally
hypothesized that the ESPT escapes the solvent dynamics control (Huppert
et al., J. Photochem. Photobiol. A2014,277, 90). The sampling of the solvent configuration
space on the ground electronic state is the first key step toward
the simulation of the ESPT event. Therefore, several configurations
in the Franck–Condon region, describing an average solvation,
were chosen as starting points for the excited-state dynamics. In
all cases, the excited-state evolution spontaneously leads to the
proton transfer event, whose rate is strongly dependent on the hydrogen
bond network around the proton acceptor solvent molecule. Our study
revealed that the explicit representation at least of three solvation
shells around the proton acceptor molecule is necessary to stabilize
the excess proton. Furthermore, the analysis of the solvent molecule
motions in proximity of the reaction site suggested that even in the
case of the strongest photoacid, the ESPT is actually assisted by
the solvation dynamics of the first and second solvation shells of
the water accepting molecule.
Collapse
Affiliation(s)
- Umberto Raucci
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S.Angelo, via Cintia, I-80126 Napoli, Italy
| | - Maria Gabriella Chiariello
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S.Angelo, via Cintia, I-80126 Napoli, Italy
| | - Nadia Rega
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S.Angelo, via Cintia, I-80126 Napoli, Italy.,CRIB, Centro Interdipartimentale di Ricerca sui Biomateriali, Piazzale Tecchio, I-80125 Napoli, Italy
| |
Collapse
|
73
|
Mai S, González L. Molecular Photochemistry: Recent Developments in Theory. Angew Chem Int Ed Engl 2020; 59:16832-16846. [PMID: 32052547 PMCID: PMC7540682 DOI: 10.1002/anie.201916381] [Citation(s) in RCA: 86] [Impact Index Per Article: 21.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2019] [Revised: 02/12/2020] [Indexed: 12/16/2022]
Abstract
Photochemistry is a fascinating branch of chemistry that is concerned with molecules and light. However, the importance of simulating light-induced processes is reflected also in fields as diverse as biology, material science, and medicine. This Minireview highlights recent progress achieved in theoretical chemistry to calculate electronically excited states of molecules and simulate their photoinduced dynamics, with the aim of reaching experimental accuracy. We focus on emergent methods and give selected examples that illustrate the progress in recent years towards predicting complex electronic structures with strong correlation, calculations on large molecules, describing multichromophoric systems, and simulating non-adiabatic molecular dynamics over long time scales, for molecules in the gas phase or in complex biological environments.
Collapse
Affiliation(s)
- Sebastian Mai
- Photonics InstituteVienna University of TechnologyGusshausstrasse 27–291040ViennaAustria
| | - Leticia González
- Institute of Theoretical ChemistryFaculty of ChemistryUniversity of ViennaWähringer Strasse 171090ViennaAustria
| |
Collapse
|
74
|
Roca M, Williams IH. Transition-State Vibrational Analysis and Isotope Effects for COMT-Catalyzed Methyl Transfer. J Am Chem Soc 2020; 142:15548-15559. [PMID: 32812761 PMCID: PMC7498148 DOI: 10.1021/jacs.0c07344] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Isotopic partition-function ratios (IPFRs) computed for transition structures (TSs) of the methyl-transfer reaction catalyzed by catechol O-methyltransferase and modeled by hybrid QM/MM methods are analyzed. The ability of smaller Hessians to reproduce trends in α-3H3 and 14Cα IPFRs as obtained using the much larger subset QM/MM Hessians from which they are extracted is investigated critically. A 6-atom-extracted Hessian reproduces perfectly the α-T3 IPFR values from the full-subset Hessians of all the TSs but not the α-14CIPFRs. Average AM1/OPLS-AA harmonic frequencies and mean-square amplitudes are presented for the 12 normal modes of the α-CH3 moiety within the active site of several enzymic transition structures, together with QM/MM potential energy scans along each of these modes to assess the degree of anharmonicity. A novel investigation of ponderal effects upon IPFRs suggests that the value for α-14C tends toward a limiting minimum whereas that for α-T3 tends toward a limiting maximum as the mass of the rest of the system increases. The transition vector is dominated by motions of atoms within the donor and acceptor moieties and is very well described as a simple combination of Walden-inversion "umbrella" bending and asymmetric stretching of the SCα and CαO bonds. The contribution of atoms of the protein residues Met40, Tyr68, and Asp141 to the transition vector is extremely small. Average valence force constants for the COMT TS show significant differences from early BEBOVIB estimates which were used in support of the compression hypothesis for catalysis. There is no correlation between TS IPFRs and the nonbonded distances for close contacts between the S atom of SAM and Tyr68 or between any of the H atoms of the transferring methyl group and either Met40 or Asp141.
Collapse
Affiliation(s)
- Maite Roca
- Departament de Química Física i Analítica, Universitat Jaume I, 12071 Castellón, Spain
| | - Ian H Williams
- Department of Chemistry, University of Bath, Bath BA2 7AY, United Kingdom
| |
Collapse
|
75
|
Mancuso JL, Mroz AM, Le KN, Hendon CH. Electronic Structure Modeling of Metal-Organic Frameworks. Chem Rev 2020; 120:8641-8715. [PMID: 32672939 DOI: 10.1021/acs.chemrev.0c00148] [Citation(s) in RCA: 97] [Impact Index Per Article: 24.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Owing to their molecular building blocks, yet highly crystalline nature, metal-organic frameworks (MOFs) sit at the interface between molecule and material. Their diverse structures and compositions enable them to be useful materials as catalysts in heterogeneous reactions, electrical conductors in energy storage and transfer applications, chromophores in photoenabled chemical transformations, and beyond. In all cases, density functional theory (DFT) and higher-level methods for electronic structure determination provide valuable quantitative information about the electronic properties that underpin the functions of these frameworks. However, there are only two general modeling approaches in conventional electronic structure software packages: those that treat materials as extended, periodic solids, and those that treat materials as discrete molecules. Each approach has features and benefits; both have been widely employed to understand the emergent chemistry that arises from the formation of the metal-organic interface. This Review canvases these approaches to date, with emphasis placed on the application of electronic structure theory to explore reactivity and electron transfer using periodic, molecular, and embedded models. This includes (i) computational chemistry considerations such as how functional, k-grid, and other model variables are selected to enable insights into MOF properties, (ii) extended solid models that treat MOFs as materials rather than molecules, (iii) the mechanics of cluster extraction and subsequent chemistry enabled by these molecular models, (iv) catalytic studies using both solids and clusters thereof, and (v) embedded, mixed-method approaches, which simulate a fraction of the material using one level of theory and the remainder of the material using another dissimilar theoretical implementation.
Collapse
Affiliation(s)
- Jenna L Mancuso
- Department of Chemistry and Biochemistry, University of Oregon, Eugene, Oregon 97405, United States
| | - Austin M Mroz
- Department of Chemistry and Biochemistry, University of Oregon, Eugene, Oregon 97405, United States
| | - Khoa N Le
- Department of Chemistry and Biochemistry, University of Oregon, Eugene, Oregon 97405, United States
| | - Christopher H Hendon
- Department of Chemistry and Biochemistry, University of Oregon, Eugene, Oregon 97405, United States
| |
Collapse
|
76
|
Mai S, González L. Molekulare Photochemie: Moderne Entwicklungen in der theoretischen Chemie. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.201916381] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Sebastian Mai
- Institut für Photonik Technische Universität Wien Gußhausstraße 27–29 1040 Wien Österreich
| | - Leticia González
- Institut für theoretische Chemie Fakultät für Chemie Universität Wien Währinger Straße 17 1090 Wien Österreich
| |
Collapse
|
77
|
Otero-de-la-Roza A, DiLabio GA. Improved Basis-Set Incompleteness Potentials for Accurate Density-Functional Theory Calculations in Large Systems. J Chem Theory Comput 2020; 16:4176-4191. [PMID: 32470304 DOI: 10.1021/acs.jctc.0c00102] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The accurate calculation of chemical properties using density-functional theory (DFT) requires the use of a nearly complete basis set. In chemical systems involving hundreds to thousands of atoms, the cost of the calculations place practical limitations on the number of basis functions that can be used. Therefore, in most practical applications of DFT to large systems, there exists a basis-set incompleteness error (BSIE). In this article, we present the next iteration of the basis-set incompleteness potentials (BSIPs), one-electron potentials designed to correct for basis-set incompleteness error. The ultimate goal associated with the development of BSIPs is to allow the calculation of molecular properties using DFT with near-complete-basis-set results at a computational cost that is similar to a small basis set calculation. In this work, we develop BSIPs for 10 atoms in the first and second rows (H, B-F, Si-Cl) and 15 common basis sets of the Pople, Dunning, Karlsruhe, and Huzinaga types. Our new BSIPs are constructed to minimize BSIE in the calculation of reaction energies, barrier heights, noncovalent binding energies, and intermolecular distances. The BSIPs were obtained using a training set of 15 944 data points. The fitting approach employed a regularized linear least-squares method with variable selection (the LASSO method), which results in a much better fit to the training data than our previous BSIPs while, at the same time, reducing the computational cost of BSIP development. The proposed BSIPs are tested on various benchmark sets and demonstrate excellent performance in practice. Our new BSIPs are also transferable; i.e., they can be used to correct BSIE in calculations that employ density functionals other than the one used in the BSIP development (B3LYP). Finally, BSIPs can be used in any quantum chemistry program that have implemented effective-core potentials without changes to the software.
Collapse
Affiliation(s)
- A Otero-de-la-Roza
- Departamento de Quı́mica Fı́sica y Analítica and MALTA-Consolider Team, Facultad de Quı́mica, Universidad de Oviedo, 33006 Oviedo, Spain
| | - Gino A DiLabio
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada.,Faculty of Management, University of British Columbia, Okanagan, 1137 Alumni Avenue, Kelowna, British Columbia V1V 1V7, Canada
| |
Collapse
|
78
|
Zhang J, Balsbaugh JL, Gao S, Ahn NG, Klinman JP. Hydrogen deuterium exchange defines catalytically linked regions of protein flexibility in the catechol O-methyltransferase reaction. Proc Natl Acad Sci U S A 2020; 117:10797-10805. [PMID: 32371482 PMCID: PMC7245127 DOI: 10.1073/pnas.1917219117] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Abstract
Human catechol O-methyltransferase (COMT) has emerged as a model for understanding enzyme-catalyzed methyl transfer from S-adenosylmethionine (AdoMet) to small-molecule catecholate acceptors. Mutation of a single residue (tyrosine 68) behind the methyl-bearing sulfonium of AdoMet was previously shown to impair COMT activity by interfering with methyl donor-acceptor compaction within the activated ground state of the wild type enzyme [J. Zhang, H. J. Kulik, T. J. Martinez, J. P. Klinman, Proc. Natl. Acad. Sci. U.S.A. 112, 7954-7959 (2015)]. This predicts the involvement of spatially defined protein dynamical effects that further tune the donor/acceptor distance and geometry as well as the electrostatics of the reactants. Here, we present a hydrogen/deuterium exchange (HDX)-mass spectrometric study of wild type and mutant COMT, comparing temperature dependences of HDX against corresponding kinetic and cofactor binding parameters. The data show that the impaired Tyr68Ala mutant displays similar breaks in Arrhenius plots of both kinetic and HDX properties that are absent in the wild type enzyme. The spatial resolution of HDX below a break point of 15-20 °C indicates changes in flexibility across ∼40% of the protein structure that is confined primarily to the periphery of the AdoMet binding site. Above 20 °C, Tyr68Ala behaves more like WT in HDX, but its rate and enthalpic barrier remain significantly altered. The impairment of catalysis by Tyr68Ala can be understood in the context of a mutationally induced alteration in protein motions that becomes manifest along and perpendicular to the primary group transfer coordinate.
Collapse
Affiliation(s)
- Jianyu Zhang
- Department of Chemistry, University of California, Berkeley, CA 94720
- The California Institute for Quantitative Biosciences (QB3), University of California, Berkeley, CA 94720
| | - Jeremy L Balsbaugh
- Department of Biochemistry, University of Colorado Boulder, Boulder, CO 80309
- BioFrontiers Institute, University of Colorado Boulder, Boulder, CO 80309
| | - Shuaihua Gao
- Department of Chemistry, University of California, Berkeley, CA 94720
- The California Institute for Quantitative Biosciences (QB3), University of California, Berkeley, CA 94720
| | - Natalie G Ahn
- Department of Biochemistry, University of Colorado Boulder, Boulder, CO 80309;
- BioFrontiers Institute, University of Colorado Boulder, Boulder, CO 80309
| | - Judith P Klinman
- Department of Chemistry, University of California, Berkeley, CA 94720;
- The California Institute for Quantitative Biosciences (QB3), University of California, Berkeley, CA 94720
- Department of Molecular and Cell Biology, University of California, Berkeley, CA 94720
| |
Collapse
|
79
|
Gimenez-Dejoz J, Tsuchiya K, Tateishi A, Motoda Y, Kigawa T, Asano Y, Numata K. Computational study on the polymerization reaction of d-aminopeptidase for the synthesis of d-peptides. RSC Adv 2020; 10:17582-17592. [PMID: 35515590 PMCID: PMC9053604 DOI: 10.1039/d0ra01138j] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2020] [Accepted: 04/28/2020] [Indexed: 02/02/2023] Open
Abstract
Almost all natural proteins are composed exclusively of l-amino acids, and this chirality influences their properties, functions, and selectivity. Proteases can recognize proteins composed of l-amino acids but display lower selectivity for their stereoisomers, d-amino acids. Taking this as an advantage, d-amino acids can be used to develop polypeptides or biobased materials with higher biostability. Chemoenzymatic peptide synthesis is a technique that uses proteases as biocatalysts to synthesize polypeptides, and d-stereospecific proteases can be used to synthesize polypeptides incorporating d-amino acids. However, engineered proteases with modified catalytic activities are required to allow the incorporation of d-amino acids with increased efficiency. To understand the stereospecificity presented by proteases and their involvement in polymerization reactions, we studied d-aminopeptidase. This enzyme displays the ability to efficiently synthesize poly d-alanine-based peptides under mild conditions. To elucidate the mechanisms involved in the unique specificity of d-aminopeptidase, we performed quantum mechanics/molecular mechanics simulations of its polymerization reaction and determined the energy barriers presented by the chiral substrates. The enzyme faces higher activation barriers for the acylation and aminolysis reactions with the l-stereoisomer than with the d-substrate (10.7 and 17.7 kcal mol-1 higher, respectively). The simulation results suggest that changes in the interaction of the substrate with Asn155 influence the stereospecificity of the polymerization reaction.
Collapse
Affiliation(s)
- Joan Gimenez-Dejoz
- Biomacromolecules Research Team, RIKEN Center for Sustainable Resource Science 2-1 Hirosawa Wako-shi Saitama 351-0198 Japan
| | - Kousuke Tsuchiya
- Biomacromolecules Research Team, RIKEN Center for Sustainable Resource Science 2-1 Hirosawa Wako-shi Saitama 351-0198 Japan
| | - Ayaka Tateishi
- Biomacromolecules Research Team, RIKEN Center for Sustainable Resource Science 2-1 Hirosawa Wako-shi Saitama 351-0198 Japan
| | - Yoko Motoda
- Biomacromolecules Research Team, RIKEN Center for Sustainable Resource Science 2-1 Hirosawa Wako-shi Saitama 351-0198 Japan
| | - Takanori Kigawa
- Laboratory for Cellular Structural Biology, RIKEN Center for Biosystems Dynamics Research 1-7-22 Suehiro-cho, Tsurumi Yokohama 230-0045 Japan
| | - Yasuhisa Asano
- Biotechnology Research Center, Department of Biotechnology, Toyama Prefectural University 5180 Kurokawa Imizu Toyama 939-0398 Japan
| | - Keiji Numata
- Biomacromolecules Research Team, RIKEN Center for Sustainable Resource Science 2-1 Hirosawa Wako-shi Saitama 351-0198 Japan
| |
Collapse
|
80
|
Magalhães RP, Fernandes HS, Sousa SF. Modelling Enzymatic Mechanisms with QM/MM Approaches: Current Status and Future Challenges. Isr J Chem 2020. [DOI: 10.1002/ijch.202000014] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Rita P. Magalhães
- UCIBIO@REQUIMTE, BioSIMDepartamento de Biomedicina, Faculdade de Medicina da Universidade do Porto Alameda Professor Hernâni Monteiro 4200-319 Porto Portugal
| | - Henriques S. Fernandes
- UCIBIO@REQUIMTE, BioSIMDepartamento de Biomedicina, Faculdade de Medicina da Universidade do Porto Alameda Professor Hernâni Monteiro 4200-319 Porto Portugal
| | - Sérgio F. Sousa
- UCIBIO@REQUIMTE, BioSIMDepartamento de Biomedicina, Faculdade de Medicina da Universidade do Porto Alameda Professor Hernâni Monteiro 4200-319 Porto Portugal
| |
Collapse
|
81
|
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]
|
82
|
Chen J, Chan B, Shao Y, Ho J. How accurate are approximate quantum chemical methods at modelling solute-solvent interactions in solvated clusters? Phys Chem Chem Phys 2020; 22:3855-3866. [PMID: 32022044 PMCID: PMC7394230 DOI: 10.1039/c9cp06792b] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
In this paper, the performance of a wide range of DFT methods is assessed for the calculation of interaction energies of thermal clusters of a solute in water. Three different charge states (neutral, proton transfer transition state and zwitterion) of glycine were solvated by 1 to 40 water molecules as sampled from molecular dynamics simulations. While some ab initio composite methods that employ insufficiently large basis sets incurred significant errors even for a cluster containing only 5 water molecules relative to the W1X-2 benchmark, the DLPNO-CCSD(T)/CBS and DSD-PBEP86 (triple zeta basis set) levels of theory predicted very accurate interaction energies. These levels of theory were used to benchmark the performance of 16 density functionals from different rungs of Jacob's Ladder. Of the Rung 4 functionals examined, the ωB97M-V and ωB97X-V functionals stood out for predicting absolute interaction energies in 40-water clusters with mean absolute deviations (MAD) ∼4 kJ mol-1. The B3LYP-D3(BJ) functional performed exceptionally well with a MAD ∼1.7 kJ mol-1 and is the overall best performing method. Calculations of relative interaction energies allow for cancellation of systematic errors, including basis set truncation and superposition errors, and the ωB97M-V and B3LYP-D3(BJ) double zeta basis set calculations yielded relative interaction energies that are within ∼3 kJ mol-1 of the benchmark. The ONIOM approximation provides another strategy for accelerating the calculation of accurate absolute interaction energies provided that the calculations have converged with respect to the size of the "high-level-layer".
Collapse
Affiliation(s)
- Junbo Chen
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia.
| | - Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo-Machi 1-14, Nagasaki 852-8521, Japan.
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, USA
| | - Junming Ho
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia.
| |
Collapse
|
83
|
Jones LO, Mosquera MA, Schatz GC, Ratner MA. Embedding Methods for Quantum Chemistry: Applications from Materials to Life Sciences. J Am Chem Soc 2020; 142:3281-3295. [PMID: 31986877 DOI: 10.1021/jacs.9b10780] [Citation(s) in RCA: 58] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Quantum mechanical embedding methods hold the promise to transform not just the way calculations are performed, but to significantly reduce computational costs and improve scaling for macro-molecular systems containing hundreds if not thousands of atoms. The field of embedding has grown increasingly broad with many approaches of different intersecting flavors. In this perspective, we lay out the methods into two streams: QM:MM and QM:QM, showcasing the advantages and disadvantages of both. We provide a review of the literature, the underpinning theories including our contributions, and we highlight current applications with select examples spanning both materials and life sciences. We conclude with prospects and future outlook on embedding, and our view on the use of universal test case scenarios for cross-comparisons of the many available (and future) embedding theories.
Collapse
Affiliation(s)
- Leighton O Jones
- Department of Chemistry , Northwestern University , Evanston , Illinois 60208 , United States
| | - Martín A Mosquera
- Department of Chemistry , Northwestern University , Evanston , Illinois 60208 , United States
| | - George C Schatz
- Department of Chemistry , Northwestern University , Evanston , Illinois 60208 , United States
| | - Mark A Ratner
- Department of Chemistry , Northwestern University , Evanston , Illinois 60208 , United States
| |
Collapse
|
84
|
Dasgupta S, Herbert JM. Using Atomic Confining Potentials for Geometry Optimization and Vibrational Frequency Calculations in Quantum-Chemical Models of Enzyme Active Sites. J Phys Chem B 2020; 124:1137-1147. [PMID: 31986049 DOI: 10.1021/acs.jpcb.9b11060] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Quantum-chemical studies of enzymatic reaction mechanisms sometimes use truncated active-site models as simplified alternatives to mixed quantum mechanics molecular mechanics (QM/MM) procedures. Eliminating the MM degrees of freedom reduces the complexity of the sampling problem, but the trade-off is the need to introduce geometric constraints in order to prevent structural collapse of the model system during geometry optimizations that do not contain a full protein backbone. These constraints may impair the efficiency of the optimization, and care must be taken to avoid artifacts such as imaginary vibrational frequencies. We introduce a simple alternative in which terminal atoms of the model system are placed in soft harmonic confining potentials rather than being rigidly constrained. This modification is simple to implement and straightforward to use in vibrational frequency calculations, unlike iterative constraint-satisfaction algorithms, and allows the optimization to proceed without constraint even though the practical result is to fix the anchor atoms in space. The new approach is more efficient for optimizing minima and transition states, as compared to the use of fixed-atom constraints, and also more robust against unwanted imaginary frequencies. We illustrate the method by application to several enzymatic reaction pathways where entropy makes a significant contribution to the relevant reaction barriers. The use of confining potentials correctly describes reaction paths and facilitates calculation of both vibrational zero-point and finite-temperature entropic corrections to barrier heights.
Collapse
Affiliation(s)
- Saswata Dasgupta
- Department of Chemistry and Biochemistry , The Ohio State University , Columbus , Ohio 43210 , United States
| | - John M Herbert
- Department of Chemistry and Biochemistry , The Ohio State University , Columbus , Ohio 43210 , United States
| |
Collapse
|
85
|
da Costa CHS, Bonatto V, Dos Santos AM, Lameira J, Leitão A, Montanari CA. Evaluating QM/MM Free Energy Surfaces for Ranking Cysteine Protease Covalent Inhibitors. J Chem Inf Model 2020; 60:880-889. [PMID: 31944110 DOI: 10.1021/acs.jcim.9b00847] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
One tactic for cysteine protease inhibition is to form a covalent bond between an electrophilic atom of the inhibitor and the thiol of the catalytic cysteine. In this study, we evaluate the reaction free energy obtained from a hybrid quantum mechanical/molecular mechanical (QM/MM) free energy profile as a predictor of affinity for reversible, covalent inhibitors of rhodesain. We demonstrate that the reaction free energy calculated with the PM6/MM potential is in agreement with the experimental data and suggest that the free energy profile for covalent bond formation in a protein environment may be a useful tool for the inhibitor design.
Collapse
Affiliation(s)
- Clauber H S da Costa
- Laboratório de Planejamento e Desenvolvimento de Fármacos , Universidade Federal do Pará , Rua Augusto Correa S/N , 66075-110 Belém , PA , Brazil
| | - Vinícius Bonatto
- Grupo de Quı́mica Medicinal do Instituto de Quı́mica de São Carlos da , Universidade de São Paulo, NEQUIMED/IQSC/USP , 13566-590 São Carlos , SP , Brazil
| | - Alberto M Dos Santos
- Laboratório de Planejamento e Desenvolvimento de Fármacos , Universidade Federal do Pará , Rua Augusto Correa S/N , 66075-110 Belém , PA , Brazil
| | - Jerônimo Lameira
- Laboratório de Planejamento e Desenvolvimento de Fármacos , Universidade Federal do Pará , Rua Augusto Correa S/N , 66075-110 Belém , PA , Brazil.,Grupo de Quı́mica Medicinal do Instituto de Quı́mica de São Carlos da , Universidade de São Paulo, NEQUIMED/IQSC/USP , 13566-590 São Carlos , SP , Brazil
| | - Andrei Leitão
- Grupo de Quı́mica Medicinal do Instituto de Quı́mica de São Carlos da , Universidade de São Paulo, NEQUIMED/IQSC/USP , 13566-590 São Carlos , SP , Brazil
| | - Carlos A Montanari
- Grupo de Quı́mica Medicinal do Instituto de Quı́mica de São Carlos da , Universidade de São Paulo, NEQUIMED/IQSC/USP , 13566-590 São Carlos , SP , Brazil
| |
Collapse
|
86
|
Brunken C, Reiher M. Self-Parametrizing System-Focused Atomistic Models. J Chem Theory Comput 2020; 16:1646-1665. [DOI: 10.1021/acs.jctc.9b00855] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Christoph Brunken
- Laboratory for Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Markus Reiher
- Laboratory for Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
87
|
Abstract
Quantum mechanics (QM) methods provide a fine description of receptor-ligand interactions and of chemical reactions. Their use in drug design and drug discovery is increasing, especially for complex systems including metal ions in the binding sites, for the design of highly selective inhibitors, for the optimization of bi-specific compounds, to understand enzymatic reactions, and for the study of covalent ligands and prodrugs. They are also used for generating molecular descriptors for predictive QSAR/QSPR models and for the parameterization of force fields. Thanks to the continuous increase of computational power offered by GPUs and to the development of sophisticated algorithms, QM methods are becoming part of the standard tools used in computer-aided drug design (CADD). We present the most used QM methods and software packages, and we discuss recent representative applications in drug design and drug discovery.
Collapse
Affiliation(s)
- Martin Kotev
- Global Research Informatics/Cheminformatics and Drug Design, Evotec (France) SAS, Toulouse, France
| | - Laurie Sarrat
- Global Research Informatics/Cheminformatics and Drug Design, Evotec (France) SAS, Toulouse, France
| | | |
Collapse
|
88
|
Zuehlsdorff TJ, Hong H, Shi L, Isborn CM. Influence of Electronic Polarization on the Spectral Density. J Phys Chem B 2019; 124:531-543. [DOI: 10.1021/acs.jpcb.9b10250] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Tim J. Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Hanbo Hong
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Liang Shi
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Christine M. Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
89
|
von der Esch B, Dietschreit JCB, Peters LDM, Ochsenfeld C. Finding Reactive Configurations: A Machine Learning Approach for Estimating Energy Barriers Applied to Sirtuin 5. J Chem Theory Comput 2019; 15:6660-6667. [PMID: 31765138 DOI: 10.1021/acs.jctc.9b00876] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Sirtuin 5 is a class III histone deacetylase that, unlike its classification, mainly catalyzes desuccinylation and demanoylation reactions. It is an interesting drug target that we use here to test new ideas for calculating reaction pathways of large molecular systems such as enzymes. A major issue with most schemes (e.g., adiabatic mapping) is that the resulting activation barrier height heavily depends on the chosen educt conformation. This makes the selection of the initial structure decisive for the success of the characterization. Here, we apply machine learning to a large number of molecular dynamics frames and potential energy barriers obtained by quantum mechanics/molecular mechanics calculations in order to identify (1) suitable start-conformations for reaction path calculations and (2) structural features relevant for the first step of the desuccinylation reaction catalyzed by Sirtuin 5. The latter generally aids the understanding of reaction mechanisms and important interactions in active centers. Using our novel approach, we found eleven key features that govern the reactivity. We were able to estimate reaction barriers with a mean absolute error of 3.6 kcal/mol and identified reactive configurations.
Collapse
Affiliation(s)
- Beatriz von der Esch
- Chair of Theoretical Chemistry, Department of Chemistry , University of Munich (LMU) , Butenandtstr. 7 , D-81377 München , Germany
| | - Johannes C B Dietschreit
- Chair of Theoretical Chemistry, Department of Chemistry , University of Munich (LMU) , Butenandtstr. 7 , D-81377 München , Germany
| | - Laurens D M Peters
- Chair of Theoretical Chemistry, Department of Chemistry , University of Munich (LMU) , Butenandtstr. 7 , D-81377 München , Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry , University of Munich (LMU) , Butenandtstr. 7 , D-81377 München , Germany
| |
Collapse
|
90
|
Summers TJ, Daniel BP, Cheng Q, DeYonker NJ. Quantifying Inter-Residue Contacts through Interaction Energies. J Chem Inf Model 2019; 59:5034-5044. [DOI: 10.1021/acs.jcim.9b00804] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Thomas J. Summers
- The Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, United States
| | - Baty P. Daniel
- The Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, United States
| | - Qianyi Cheng
- The Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, United States
| | - Nathan J. DeYonker
- The Department of Chemistry, The University of Memphis, 213 Smith Chemistry Building, Memphis, Tennessee 38152-3550, United States
| |
Collapse
|
91
|
Abstract
For decades, there has been debate regarding the origin of the catalytic power of enzymes. In this work, we use the approach of computational chemistry to study the enzyme catechol O-methyltransferase (COMT) and reveal that the two current views on the catalytic mechanism of enzymes, the rate-promoting vibrations and the electric field, may both be viewed as part of the chemical step catalyzed by COMT. However, we show that the rate-promoting vibrations cause the electrostatic effect. This work provides insight into the catalytic mechanism of COMT and resolves a longstanding controversy regarding this enzyme's mechanism.
Collapse
Affiliation(s)
- Xi Chen
- Department of Chemistry and Biochemistry, University of Arizona, 1306 East University Boulevard, Tucson, Arizona 85721, United States
| | - Steven D. Schwartz
- Department of Chemistry and Biochemistry, University of Arizona, 1306 East University Boulevard, Tucson, Arizona 85721, United States
| |
Collapse
|
92
|
Loco D, Lagardère L, Cisneros GA, Scalmani G, Frisch M, Lipparini F, Mennucci B, Piquemal JP. Towards large scale hybrid QM/MM dynamics of complex systems with advanced point dipole polarizable embeddings. Chem Sci 2019; 10:7200-7211. [PMID: 31588288 PMCID: PMC6677116 DOI: 10.1039/c9sc01745c] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2019] [Accepted: 06/08/2019] [Indexed: 12/17/2022] Open
Abstract
In this work, we present a general route to hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) Molecular Dynamics for complex systems using a polarizable embedding. We extend the capabilities of our hybrid framework, combining the Gaussian and Tinker/Tinker-HP packages in the context of the AMOEBA polarizable force field to treat large (bio)systems where the QM and the MM subsystems are covalently bound, adopting pseudopotentials at the boundaries between the two regions. We discuss in detail the implementation and demonstrate the global energy conservation of our QM/MM Born-Oppenheimer molecular dynamics approach using Density Functional Theory. Finally, the approach is assessed on the electronic absorption properties of a 16 500 atom complex encompassing an organic dye embedded in a DNA matrix in solution, extending the hybrid method to a time-dependent Density Functional Theory approach. The results obtained comparing different partitions between the quantum and the classical subsystems also suggest that large QM portions are not necessary if accurate polarizable force fields are used in a variational formulation of the embedding, properly including the QM/MM mutual polarization.
Collapse
Affiliation(s)
- Daniele Loco
- Sorbonne Université , CNRS , Laboratoire de Chimie Théorique, LCT , Paris , France . ;
| | - Louis Lagardère
- Sorbonne Université , CNRS , Institut Parisien de Chimie Physique et Théorique, IP2CT , Paris , France
- Sorbonne Université , Institut des Sciences du Calcul et des Données, ISCD , Paris , France
| | | | | | | | - Filippo Lipparini
- Univerisita di Pisa , Dipartimento di Chimica e ChimicaIndustriale , Pisa , Italy
| | - Benedetta Mennucci
- Univerisita di Pisa , Dipartimento di Chimica e ChimicaIndustriale , Pisa , Italy
| | - Jean-Philip Piquemal
- Sorbonne Université , CNRS , Laboratoire de Chimie Théorique, LCT , Paris , France . ;
- Institut Universitaire de France, IUF , Paris , France
- The University of Texas at Austin , Department of Biomedical Engineering , TX , USA
| |
Collapse
|
93
|
Yang Z, Liu F, Steeves AH, Kulik HJ. Quantum Mechanical Description of Electrostatics Provides a Unified Picture of Catalytic Action Across Methyltransferases. J Phys Chem Lett 2019; 10:3779-3787. [PMID: 31244268 DOI: 10.1021/acs.jpclett.9b01555] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Methyl transferases (MTases) are a well-studied class of enzymes for which competing enzymatic enhancement mechanisms have been suggested, ranging from structural methyl group CH···X hydrogen bonds (HBs) to electrostatic- and charge-transfer-driven stabilization of the transition state (TS). We identified all Class I MTases for which reasonable resolution (<2.0 Å) crystal structures could be used to form catalytically competent ternary complexes for multiscale (i.e., quantum-mechanical/molecular-mechanical or QM/MM) simulation of the SN2 methyl transfer reaction coordinate. The four Class I MTases studied have both distinct functions (e.g., protein repair or biosynthesis) and substrate nucleophiles (i.e., C, N, or O). While CH···X HBs stabilize all reactant complexes, no universal TS stabilization role is found for these interactions in MTases. A consistent picture is instead obtained through analysis of charge transfer and electrostatics, wherein much of cofactor-substrate charge separation is maintained in the TS region, and electrostatic potential is correlated with substrate nucleophilicity (i.e., intrinsic reactivity).
Collapse
Affiliation(s)
- Zhongyue Yang
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Fang Liu
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Adam H Steeves
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Heather J Kulik
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| |
Collapse
|
94
|
Ranaghan KE, Shchepanovska D, Bennie SJ, Lawan N, Macrae SJ, Zurek J, Manby FR, Mulholland AJ. Projector-Based Embedding Eliminates Density Functional Dependence for QM/MM Calculations of Reactions in Enzymes and Solution. J Chem Inf Model 2019; 59:2063-2078. [PMID: 30794388 DOI: 10.1021/acs.jcim.8b00940] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Combined quantum mechanics/molecular mechanics (QM/MM) methods are increasingly widely utilized in studies of reactions in enzymes and other large systems. Here, we apply a range of QM/MM methods to investigate the Claisen rearrangement of chorismate to prephenate, in solution, and in the enzyme chorismate mutase. Using projector-based embedding in a QM/MM framework, we apply treatments up to the CCSD(T) level. We test a range of density functional QM/MM methods and QM region sizes. The results show that the calculated reaction energetics are significantly more sensitive to the choice of density functional than they are to the size of the QM region in these systems. Projector-based embedding of a wave function method in DFT reduced the 13 kcal/mol spread in barrier heights calculated at the DFT/MM level to a spread of just 0.3 kcal/mol, essentially eliminating dependence on the functional. Projector-based embedding of correlated ab initio methods provides a practical method for achieving high accuracy for energy profiles derived from DFT and DFT/MM calculations for reactions in condensed phases.
Collapse
Affiliation(s)
- Kara E Ranaghan
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Bristol , U.K. BS8 1TS
| | - Darya Shchepanovska
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Bristol , U.K. BS8 1TS
| | - Simon J Bennie
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Bristol , U.K. BS8 1TS
| | - Narin Lawan
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Bristol , U.K. BS8 1TS
| | - Stephen J Macrae
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Bristol , U.K. BS8 1TS
| | - Jolanta Zurek
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Bristol , U.K. BS8 1TS
| | - Frederick R Manby
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Bristol , U.K. BS8 1TS
| | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry , University of Bristol , Bristol , U.K. BS8 1TS
| |
Collapse
|
95
|
Czarnota S, Johannissen LO, Baxter NJ, Rummel F, Wilson AL, Cliff MJ, Levy CW, Scrutton NS, Waltho JP, Hay S. Equatorial Active Site Compaction and Electrostatic Reorganization in Catechol- O-methyltransferase. ACS Catal 2019; 9:4394-4401. [PMID: 31080692 PMCID: PMC6503465 DOI: 10.1021/acscatal.9b00174] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2019] [Revised: 03/26/2019] [Indexed: 12/18/2022]
Abstract
Catechol-O-methyltransferase (COMT) is a model S-adenosyl-l-methionine (SAM) dependent methyl transferase, which catalyzes the methylation of catecholamine neurotransmitters such as dopamine in the primary pathway of neurotransmitter deactivation in animals. Despite extensive study, there is no consensus view of the physical basis of catalysis in COMT. Further progress requires experimental data that directly probes active site geometry, protein dynamics and electrostatics, ideally in a range of positions along the reaction coordinate. Here we establish that sinefungin, a fungal-derived inhibitor of SAM-dependent enzymes that possess transition state-like charge on the transferring group, can be used as a transition state analog of COMT when combined with a catechol. X-ray crystal structures and NMR backbone assignments of the ternary complexes of the soluble form of human COMT containing dinitrocatechol, Mg2+ and SAM or sinefungin were determined. Comparison and further analysis with the aid of density functional theory calculations and molecular dynamics simulations provides evidence for active site "compaction", which is driven by electrostatic stabilization between the transferring methyl group and "equatorial" active site residues that are orthogonal to the donor-acceptor (pseudo reaction) coordinate. We propose that upon catecholamine binding and subsequent proton transfer to Lys 144, the enzyme becomes geometrically preorganized, with little further movement along the donor-acceptor coordinate required for methyl transfer. Catalysis is then largely facilitated through stabilization of the developing charge on the transferring methyl group via "equatorial" H-bonding and electrostatic interactions orthogonal to the donor-acceptor coordinate.
Collapse
Affiliation(s)
- Sylwia Czarnota
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
- School
of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom
| | - Linus O. Johannissen
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
| | - Nicola J. Baxter
- Krebs
Institute for Biomolecular Research, Department of Molecular Biology
and Biotechnology, The University of Sheffield, Firth Court, Western Bank, Sheffield, S10 2TN, United Kingdom
| | - Felix Rummel
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
- School
of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom
| | - Alex L. Wilson
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
- School
of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom
| | - Matthew J. Cliff
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
| | - Colin W. Levy
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
| | - Nigel S. Scrutton
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
- School
of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom
| | - Jonathan P. Waltho
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
- School
of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom
- Krebs
Institute for Biomolecular Research, Department of Molecular Biology
and Biotechnology, The University of Sheffield, Firth Court, Western Bank, Sheffield, S10 2TN, United Kingdom
| | - Sam Hay
- Manchester
Institute of Biotechnology, The University
of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
- School
of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, United Kingdom
| |
Collapse
|
96
|
Mehmood R, Qi HW, Steeves AH, Kulik HJ. The Protein’s Role in Substrate Positioning and Reactivity for Biosynthetic Enzyme Complexes: The Case of SyrB2/SyrB1. ACS Catal 2019. [DOI: 10.1021/acscatal.9b00865] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
|
97
|
Cuyàs E, Verdura S, Lozano-Sánchez J, Viciano I, Llorach-Parés L, Nonell-Canals A, Bosch-Barrera J, Brunet J, Segura-Carretero A, Sanchez-Martinez M, Encinar JA, Menendez JA. The extra virgin olive oil phenolic oleacein is a dual substrate-inhibitor of catechol-O-methyltransferase. Food Chem Toxicol 2019; 128:35-45. [PMID: 30935952 DOI: 10.1016/j.fct.2019.03.049] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2019] [Revised: 03/25/2019] [Accepted: 03/26/2019] [Indexed: 12/13/2022]
Abstract
Catechol-containing polyphenols present in coffee and tea, while serving as excellent substrates for catechol-O-methyltransferase (COMT)-catalyzed O-methylation, can also operate as COMT inhibitors. However, little is known about the relationship between COMT and the characteristic phenolics present in extra virgin olive oil (EVOO). We here selected the EVOO dihydroxy-phenol oleacein for a computational study of COMT-driven methylation using classic molecular docking/molecular dynamics simulations and hybrid quantum mechanical/molecular mechanics, which were supported by in vitro activity studies using human COMT. Oleacein could be superimposed onto the catechol-binding site of COMT, maintaining the interactions with the atomic positions involved in methyl transfer from the S-adenosyl-L-methionine cofactor. The transition state structure for the meta-methylation in the O5 position of the oleacein benzenediol moiety was predicted to occur preferentially. Enzyme analysis of the conversion ratio of catechol to O-alkylated guaiacol confirmed the inhibitory effect of oleacein on human COMT, which remained unaltered when tested against the protein version encoded by the functional Val158Met polymorphism of the COMT gene. Our study provides a theoretical determination of how EVOO dihydroxy-phenols can be metabolized via COMT. The ability of oleacein to inhibit COMT adds a new dimension to the physiological and therapeutic utility of EVOO secoiridoids.
Collapse
Affiliation(s)
- Elisabet Cuyàs
- ProCURE (Program Against Cancer Therapeutic Resistance), Metabolism & Cancer Group, Catalan Institute of Oncology, Girona, Spain; Girona Biomedical Research Institute (IDIBGI), Girona, Spain
| | - Sara Verdura
- ProCURE (Program Against Cancer Therapeutic Resistance), Metabolism & Cancer Group, Catalan Institute of Oncology, Girona, Spain; Girona Biomedical Research Institute (IDIBGI), Girona, Spain
| | - Jesús Lozano-Sánchez
- Department of Analytical Chemistry, Faculty of Sciences, University of Granada, Granada, Spain; Research and Development Functional Food Centre (CIDAF), PTS Granada, Granada, Spain
| | | | | | | | - Joaquim Bosch-Barrera
- Girona Biomedical Research Institute (IDIBGI), Girona, Spain; Medical Oncology, Catalan Institute of Oncology (ICO) Dr. Josep Trueta University Hospital, Girona, Spain; Department of Medical Sciences, Medical School University of Girona, Girona, Spain
| | - Joan Brunet
- Medical Oncology, Catalan Institute of Oncology (ICO) Dr. Josep Trueta University Hospital, Girona, Spain; Department of Medical Sciences, Medical School University of Girona, Girona, Spain; Hereditary Cancer Programme, Catalan Institute of Oncology (ICO), Bellvitge Institute for Biomedical Research (IDIBELL) L'Hospitalet del Llobregat, Barcelona, Spain; Hereditary Cancer Programme, Catalan Institute of Oncology (ICO) Girona Biomedical Research Institute (IDIBGI), Girona, Spain
| | - Antonio Segura-Carretero
- Department of Analytical Chemistry, Faculty of Sciences, University of Granada, Granada, Spain; Research and Development Functional Food Centre (CIDAF), PTS Granada, Granada, Spain
| | | | - José Antonio Encinar
- Institute of Research, Development and Innovation in Biotechnology of Elche (IDiBE) and Molecular and Cell Biology Institute (IBMC), Miguel Hernández University (UMH), Elche, Spain.
| | - Javier A Menendez
- ProCURE (Program Against Cancer Therapeutic Resistance), Metabolism & Cancer Group, Catalan Institute of Oncology, Girona, Spain; Girona Biomedical Research Institute (IDIBGI), Girona, Spain.
| |
Collapse
|
98
|
Qi HW, Kulik HJ. Evaluating Unexpectedly Short Non-covalent Distances in X-ray Crystal Structures of Proteins with Electronic Structure Analysis. J Chem Inf Model 2019; 59:2199-2211. [DOI: 10.1021/acs.jcim.9b00144] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Helena W. Qi
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
99
|
Kitakawa CK, Maruyama T, Oonari J, Mitsuta Y, Kawakami T, Okumura M, Yamaguchi K, Yamanaka S. Linear Response Functions of Densities and Spin Densities for Systematic Modeling of the QM/MM Approach for Mono- and Poly-Nuclear Transition Metal Systems. Molecules 2019; 24:E821. [PMID: 30823580 PMCID: PMC6412466 DOI: 10.3390/molecules24040821] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2018] [Revised: 01/28/2019] [Accepted: 02/15/2019] [Indexed: 12/02/2022] Open
Abstract
We applied our analysis, based on a linear response function of density and spin density, to two typical transition metal complex systems-the reaction centers of P450, and oxygen evolving center in Photosystem II, both of which contain open-shell transition metal ions. We discuss the relationship between LRF of electron density and spin density and the types of units and interactions of the systems. The computational results are discussed in relation to quantum mechanics (QM) cluster and quantum mechanics/molecular mechanics (QM/MM) modeling that are employed to compute the reaction centers of enzymes.
Collapse
Affiliation(s)
- Colin K Kitakawa
- Graduate School of Science, Osaka University, Osaka 565-0871, Japan.
| | - Tomohiro Maruyama
- Graduate School of Science, Osaka University, Osaka 565-0871, Japan.
| | - Jinta Oonari
- Graduate School of Science, Osaka University, Osaka 565-0871, Japan.
| | - Yuki Mitsuta
- Graduate School of Science, Osaka University, Osaka 565-0871, Japan.
| | - Takashi Kawakami
- Graduate School of Science, Osaka University, Osaka 565-0871, Japan.
| | - Mitsutaka Okumura
- Graduate School of Science, Osaka University, Osaka 565-0871, Japan.
| | | | - Shusuke Yamanaka
- Graduate School of Science, Osaka University, Osaka 565-0871, Japan.
- Quantum information, Quantum Biology division, Institute for Open and Transdisciplinary Research Initiatives, Osaka University, Osaka 560-8531, Japan.
| |
Collapse
|
100
|
Yang Z, Mehmood R, Wang M, Qi HW, Steeves AH, Kulik HJ. Revealing quantum mechanical effects in enzyme catalysis with large-scale electronic structure simulation. REACT CHEM ENG 2019; 4:298-315. [PMID: 31572618 PMCID: PMC6768422 DOI: 10.1039/c8re00213d] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Enzymes have evolved to facilitate challenging reactions at ambient conditions with specificity seldom matched by other catalysts. Computational modeling provides valuable insight into catalytic mechanism, and the large size of enzymes mandates multi-scale, quantum mechanical-molecular mechanical (QM/MM) simulations. Although QM/MM plays an essential role in balancing simulation cost to enable sampling with full QM treatment needed to understand electronic structure in enzyme active sites, the relative importance of these two strategies for understanding enzyme mechanism is not well known. We explore challenges in QM/MM for studying the reactivity and stability of three diverse enzymes: i) Mg2+-dependent catechol O-methyltransferase (COMT), ii) radical enzyme choline trimethylamine lyase (CutC), and iii) DNA methyltransferase (DNMT1), which has structural Zn2+ binding sites. In COMT, strong non-covalent interactions lead to long range coupling of electronic structure properties across the active site, but the more isolated nature of the metallocofactor in DNMT1 leads to faster convergence of some properties. We quantify these effects in COMT by computing covariance matrices of by-residue electronic structure properties during dynamics and along the reaction coordinate. In CutC, we observe spontaneous bond cleavage following initiation events, highlighting the importance of sampling and dynamics. We use electronic structure analysis to quantify the relative importance of CHO and OHO non-covalent interactions in imparting reactivity. These three diverse cases enable us to provide some general recommendations regarding QM/MM simulation of enzymes.
Collapse
Affiliation(s)
- Zhongyue Yang
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139
| | - Rimsha Mehmood
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139
| | - Mengyi Wang
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139
| | - Helena W. Qi
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139
| | - Adam H. Steeves
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139
| | - Heather J. Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139
| |
Collapse
|