1
|
Pederson JP, McDaniel JG. PyDFT-QMMM: A modular, extensible software framework for DFT-based QM/MM molecular dynamics. J Chem Phys 2024; 161:034103. [PMID: 39007371 DOI: 10.1063/5.0219851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2024] [Accepted: 06/24/2024] [Indexed: 07/16/2024] Open
Abstract
PyDFT-QMMM is a Python-based package for performing hybrid quantum mechanics/molecular mechanics (QM/MM) simulations at the density functional level of theory. The program is designed to treat short-range and long-range interactions through user-specified combinations of electrostatic and mechanical embedding procedures within periodic simulation domains, providing necessary interfaces to external quantum chemistry and molecular dynamics software. To enable direct embedding of long-range electrostatics in periodic systems, we have derived and implemented force terms for our previously described QM/MM/PME approach [Pederson and McDaniel, J. Chem. Phys. 156, 174105 (2022)]. Communication with external software packages Psi4 and OpenMM is facilitated through Python application programming interfaces (APIs). The core library contains basic utilities for running QM/MM molecular dynamics simulations, and plug-in entry-points are provided for users to implement custom energy/force calculation and integration routines, within an extensible architecture. The user interacts with PyDFT-QMMM primarily through its Python API, allowing for complex workflow development with Python scripting, for example, interfacing with PLUMED for free energy simulations. We provide benchmarks of forces and energy conservation for the QM/MM/PME and alternative QM/MM electrostatic embedding approaches. We further demonstrate a simple example use case for water solute in a water solvent system, for which radial distribution functions are computed from 100 ps QM/MM simulations; in this example, we highlight how the solvation structure is sensitive to different basis-set choices due to under- or over-polarization of the QM water molecule's electron density.
Collapse
Affiliation(s)
- John P Pederson
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Jesse G McDaniel
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
2
|
Solov’yov AV, Verkhovtsev AV, Mason NJ, Amos RA, Bald I, Baldacchino G, Dromey B, Falk M, Fedor J, Gerhards L, Hausmann M, Hildenbrand G, Hrabovský M, Kadlec S, Kočišek J, Lépine F, Ming S, Nisbet A, Ricketts K, Sala L, Schlathölter T, Wheatley AEH, Solov’yov IA. Condensed Matter Systems Exposed to Radiation: Multiscale Theory, Simulations, and Experiment. Chem Rev 2024; 124:8014-8129. [PMID: 38842266 PMCID: PMC11240271 DOI: 10.1021/acs.chemrev.3c00902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2023] [Revised: 05/02/2024] [Accepted: 05/10/2024] [Indexed: 06/07/2024]
Abstract
This roadmap reviews the new, highly interdisciplinary research field studying the behavior of condensed matter systems exposed to radiation. The Review highlights several recent advances in the field and provides a roadmap for the development of the field over the next decade. Condensed matter systems exposed to radiation can be inorganic, organic, or biological, finite or infinite, composed of different molecular species or materials, exist in different phases, and operate under different thermodynamic conditions. Many of the key phenomena related to the behavior of irradiated systems are very similar and can be understood based on the same fundamental theoretical principles and computational approaches. The multiscale nature of such phenomena requires the quantitative description of the radiation-induced effects occurring at different spatial and temporal scales, ranging from the atomic to the macroscopic, and the interlinks between such descriptions. The multiscale nature of the effects and the similarity of their manifestation in systems of different origins necessarily bring together different disciplines, such as physics, chemistry, biology, materials science, nanoscience, and biomedical research, demonstrating the numerous interlinks and commonalities between them. This research field is highly relevant to many novel and emerging technologies and medical applications.
Collapse
Affiliation(s)
| | | | - Nigel J. Mason
- School
of Physics and Astronomy, University of
Kent, Canterbury CT2 7NH, United
Kingdom
| | - Richard A. Amos
- Department
of Medical Physics and Biomedical Engineering, University College London, London WC1E 6BT, U.K.
| | - Ilko Bald
- Institute
of Chemistry, University of Potsdam, Karl-Liebknecht-Str. 24-25, 14476 Potsdam, Germany
| | - Gérard Baldacchino
- Université
Paris-Saclay, CEA, LIDYL, 91191 Gif-sur-Yvette, France
- CY Cergy Paris Université,
CEA, LIDYL, 91191 Gif-sur-Yvette, France
| | - Brendan Dromey
- Centre
for Light Matter Interactions, School of Mathematics and Physics, Queen’s University Belfast, Belfast BT7 1NN, United Kingdom
| | - Martin Falk
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, 61200 Brno, Czech Republic
- Kirchhoff-Institute
for Physics, Heidelberg University, Im Neuenheimer Feld 227, 69120 Heidelberg, Germany
| | - Juraj Fedor
- J.
Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, 18223 Prague, Czech Republic
| | - Luca Gerhards
- Institute
of Physics, Carl von Ossietzky University, Carl-von-Ossietzky-Str. 9-11, 26129 Oldenburg, Germany
| | - Michael Hausmann
- Kirchhoff-Institute
for Physics, Heidelberg University, Im Neuenheimer Feld 227, 69120 Heidelberg, Germany
| | - Georg Hildenbrand
- Kirchhoff-Institute
for Physics, Heidelberg University, Im Neuenheimer Feld 227, 69120 Heidelberg, Germany
- Faculty
of Engineering, University of Applied Sciences
Aschaffenburg, Würzburger
Str. 45, 63743 Aschaffenburg, Germany
| | | | - Stanislav Kadlec
- Eaton European
Innovation Center, Bořivojova
2380, 25263 Roztoky, Czech Republic
| | - Jaroslav Kočišek
- J.
Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, 18223 Prague, Czech Republic
| | - Franck Lépine
- Université
Claude Bernard Lyon 1, CNRS, Institut Lumière
Matière, F-69622, Villeurbanne, France
| | - Siyi Ming
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield
Road, Cambridge CB2 1EW, United Kingdom
| | - Andrew Nisbet
- Department
of Medical Physics and Biomedical Engineering, University College London, London WC1E 6BT, U.K.
| | - Kate Ricketts
- Department
of Targeted Intervention, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Leo Sala
- J.
Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, 18223 Prague, Czech Republic
| | - Thomas Schlathölter
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh
4, 9747 AG Groningen, The Netherlands
- University
College Groningen, University of Groningen, Hoendiepskade 23/24, 9718 BG Groningen, The Netherlands
| | - Andrew E. H. Wheatley
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield
Road, Cambridge CB2 1EW, United Kingdom
| | - Ilia A. Solov’yov
- Institute
of Physics, Carl von Ossietzky University, Carl-von-Ossietzky-Str. 9-11, 26129 Oldenburg, Germany
| |
Collapse
|
3
|
Jiang H, Ryde U. Putative reaction mechanism of nitrogenase with a half-dissociated S2B ligand. Dalton Trans 2024; 53:11500-11513. [PMID: 38916132 DOI: 10.1039/d4dt00937a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/26/2024]
Abstract
We have studied whether dissociation of the S2B sulfide ligand from one of its two coordinating Fe ions may affect the later parts of the reaction mechanism of nitrogenase. Such dissociation has been shown to be favourable for the E2-E4 states in the reaction mechanism, but previous studies have assumed that S2B either remains bridging or has fully dissociated from the active-site FeMo cluster. We employ combined quantum mechanical and molecular mechanical (QM/MM) calculations with two density-functional theory methods, r2SCAN and TPSSh. To make dissociation of S2B possible, we have added a proton to this group throughout the reaction. We study the reaction starting from the E4 state with N2H2 bound to the cluster. Our results indicate that half-dissociation of S2B is unfavourable in most steps of the reaction mechanism. We observe favourable half-dissociation of S2B only when NH or NH2 is bound to the cluster, bridging Fe2 and Fe6. However, the former state is most likely not involved in the reaction mechanism and the latter state is only an intermittent intermediate of the E7 state. Therefore, half-dissociation of S2B seems to play only a minor role in the later parts of the reaction mechanism of nitrogenase. Our suggested mechanism with a protonated S2B is alternating (the two N atoms of the substrate is protonated in an alternating manner) and the substrate prefers to bind to Fe2, in contrast to the preferred binding to Fe6 observed when S2B is unprotonated and bridging Fe2 and Fe6.
Collapse
Affiliation(s)
- Hao Jiang
- Department of Computational Chemistry, Lund University, Chemical Centre, P. O. Box 124, SE-221 00 Lund, Sweden.
| | - Ulf Ryde
- Department of Computational Chemistry, Lund University, Chemical Centre, P. O. Box 124, SE-221 00 Lund, Sweden.
| |
Collapse
|
4
|
Rovaletti A, Moro G, Cosentino U, Ryde U, Greco C. CO Oxidation Mechanism of Silver-Substituted Mo/Cu CO-Dehydrogenase - Analogies and Differences to the Native Enzyme. Chemphyschem 2024; 25:e202400293. [PMID: 38631392 DOI: 10.1002/cphc.202400293] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2024] [Revised: 04/16/2024] [Accepted: 04/17/2024] [Indexed: 04/19/2024]
Abstract
The aerobic oxidation of carbon monoxide to carbon dioxide is catalysed by the Mo/Cu-containing CO-dehydrogenase enzyme in the soil bacterium Oligotropha carboxidovorans, enabling the organism to grow on the small gas molecule as carbon and energy source. It was shown experimentally that silver can be substituted for copper in the active site of Mo/Cu CODH to yield a functional enzyme. In this study, we employed QM/MM calculations to investigate whether the reaction mechanism of the silver-substituted enzyme is similar to that of the native enzyme. Our results suggest that the Ag-substituted enzyme can oxidize CO and release CO2 following the same reaction steps as the native enzyme, with a computed rate-limiting step of 10.4 kcal/mol, consistent with experimental findings. Surprisingly, lower activation energies for C-O bond formation have been found in the presence of silver. Furthermore, comparison of rate constants for reduction of copper- and silver-containing enzymes suggests a discrepancy in the transition state stabilization upon silver substitution. We also evaluated the effects that differences in the water-active site interaction may exert on the overall energy profile of catalysis. Finally, the formation of a thiocarbonate intermediate along the catalytic pathway was found to be energetically unfavorable for the Ag-substituted enzyme. This finding aligns with the hypothesis proposed for the wild-type form, suggesting that the creation of such species may not be necessary for the enzymatic catalysis of CO oxidation.
Collapse
Affiliation(s)
- Anna Rovaletti
- Department of Earth and Environmental Sciences, Milano-Bicocca University, Piazza della Scienza 1, Milano, 20126, Italy
| | - Giorgio Moro
- Department of Biotechnology and Biosciences, Milano-Bicocca University, Piazza della Scienza 2, Milano, 20126, Italy
| | - Ugo Cosentino
- Department of Earth and Environmental Sciences, Milano-Bicocca University, Piazza della Scienza 1, Milano, 20126, Italy
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00, Lund, Sweden
| | - Claudio Greco
- Department of Earth and Environmental Sciences, Milano-Bicocca University, Piazza della Scienza 1, Milano, 20126, Italy
| |
Collapse
|
5
|
Schmidt B, Hildebrandt A. From GPUs to AI and quantum: three waves of acceleration in bioinformatics. Drug Discov Today 2024; 29:103990. [PMID: 38663581 DOI: 10.1016/j.drudis.2024.103990] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2023] [Revised: 04/05/2024] [Accepted: 04/17/2024] [Indexed: 05/01/2024]
Abstract
The enormous growth in the amount of data generated by the life sciences is continuously shifting the field from model-driven science towards data-driven science. The need for efficient processing has led to the adoption of massively parallel accelerators such as graphics processing units (GPUs). Consequently, the development of bioinformatics methods nowadays often heavily depends on the effective use of these powerful technologies. Furthermore, progress in computational techniques and architectures continues to be highly dynamic, involving novel deep neural network models and artificial intelligence (AI) accelerators, and potentially quantum processing units in the future. These are expected to be disruptive for the life sciences as a whole and for drug discovery in particular. Here, we identify three waves of acceleration and their applications in a bioinformatics context: (i) GPU computing, (ii) AI and (iii) next-generation quantum computers.
Collapse
Affiliation(s)
- Bertil Schmidt
- Institut für Informatik, Johannes Gutenberg University, Mainz, Germany.
| | | |
Collapse
|
6
|
Jafari S, Ryde U, Irani M. QM/MM study of the catalytic reaction of aphid myrosinase. Int J Biol Macromol 2024; 262:130089. [PMID: 38360236 DOI: 10.1016/j.ijbiomac.2024.130089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2023] [Revised: 01/07/2024] [Accepted: 02/08/2024] [Indexed: 02/17/2024]
Abstract
Brevicoryne brassicae, an aphid species, exclusively consumes plants from the Brassicaceae family and employs a sophisticated defense mechanism involving a myrosinase enzyme that breaks down glucosinolates obtained from its host plants. In this work, we employed combined quantum mechanical and molecular mechanical (QM/MM) calculations and molecular dynamics (MD) simulations to study the catalytic reaction of aphid myrosinase. A proper QM region to study the myrosinase reaction should contain the whole substrate, models of Gln-19, His-122, Asp-124, Asn-166, Glu-167, Lys-173, Tyr-180, Val-228, Tyr-309, Tyr-346, Ile-347, Glu-374, Glu-423, Trp-424, and a water molecule. The calculations show that Asp-124 and Glu-423 must be charged, His-122 must be protonated on NE2, and Glu-167 must be protonated on OE2. Our model reproduces the anomeric retaining characteristic of myrosinase and indicates that the deglycosylation reaction is the rate-determining step of the reaction. Based on the calculations, we propose a reaction mechanism for aphid myrosinase-mediated hydrolysis of glucosinolates with an overall barrier of 15.2 kcal/mol. According to the results, removing a proton from Arg-312 or altering it to valine or methionine increases glycosylation barriers but decreases the deglycosylation barrier.
Collapse
Affiliation(s)
- Sonia Jafari
- Department of Chemistry, University of Kurdistan, P.O. Box 66175-416, Sanandaj, Iran
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Mehdi Irani
- Department of Chemistry, University of Kurdistan, P.O. Box 66175-416, Sanandaj, Iran.
| |
Collapse
|
7
|
Jiang H, Ryde U. H 2 formation from the E 2-E 4 states of nitrogenase. Phys Chem Chem Phys 2024; 26:1364-1375. [PMID: 38108422 DOI: 10.1039/d3cp05181a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2023]
Abstract
Nitrogenase is the only enzyme that can cleave the strong triple bond in N2, making nitrogen available for biological lifeforms. The active site is a MoFe7S9C cluster (the FeMo cluster) that binds eight electrons and protons during one catalytic cycle, giving rise to eight intermediate states E0-E7. It is experimentally known that N2 binds to the E4 state and that H2 is a compulsory byproduct of the reaction. However, formation of H2 is also an unproductive side reaction that should be avoided, especially in the early steps of the reaction mechanism (E2 and E3). Here, we study the formation of H2 for various structural interpretations of the E2-E4 states using combined quantum mechanical and molecular mechanical (QM/MM) calculations and four different density-functional theory methods. We find large differences in the predictions of the different methods. B3LYP strongly favours protonation of the central carbide ion and H2 cannot form from such structures. On the other hand, with TPSS, r2SCAN and TPSSh, H2 formation is strongly exothermic for all structures and En and therefore need strict kinetic control to be avoided. For the E2 state, the kinetic barriers for the low-energy structures are high enough to avoid H2 formation. However, for both the E3 and E4 states, all three methods predict that the best structure has two hydride ions bridging the same pair of Fe ions (Fe2 and Fe6) and these two ions can combine to form H2 with an activation barrier of only 29-57 kJ mol-1, corresponding to rates of 7 × 102 to 5 × 107 s-1, i.e. much faster than the turnover rate of the enzyme (1-5 s-1). We have also studied H-atom movements within the FeMo cluster, showing that the various protonation states can quite freely be interconverted (activation barriers of 12-69 kJ mol-1).
Collapse
Affiliation(s)
- Hao Jiang
- Department of Computational Chemistry, Lund University, Chemical Centre, P. O. Box 124, SE-221 00 Lund, Sweden.
| | - Ulf Ryde
- Department of Computational Chemistry, Lund University, Chemical Centre, P. O. Box 124, SE-221 00 Lund, Sweden.
| |
Collapse
|
8
|
Chaudhuri S, Logsdail AJ, Maurer RJ. Stability of Single Gold Atoms on Defective and Doped Diamond Surfaces. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2023; 127:16187-16203. [PMID: 37609382 PMCID: PMC10440818 DOI: 10.1021/acs.jpcc.3c03900] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/09/2023] [Revised: 07/20/2023] [Indexed: 08/24/2023]
Abstract
Polycrystalline boron-doped diamond (BDD) is widely used as a working electrode material in electrochemistry, and its properties, such as its stability, make it an appealing support material for nanostructures in electrocatalytic applications. Recent experiments have shown that electrodeposition can lead to the creation of stable small nanoclusters and even single gold adatoms on the BDD surfaces. We investigate the adsorption energy and kinetic stability of single gold atoms adsorbed onto an atomistic model of BDD surfaces by using density functional theory. The surface model is constructed using hybrid quantum mechanics/molecular mechanics embedding techniques and is based on an oxygen-terminated diamond (110) surface. We use the hybrid quantum mechanics/molecular mechanics method to assess the ability of different density functional approximations to predict the adsorption structure, energy, and barrier for diffusion on pristine and defective surfaces. We find that surface defects (vacancies and surface dopants) strongly anchor adatoms on vacancy sites. We further investigated the thermal stability of gold adatoms, which reveals high barriers associated with lateral diffusion away from the vacancy site. The result provides an explanation for the high stability of experimentally imaged single gold adatoms on BDD and a starting point to investigate the early stages of nucleation during metal surface deposition.
Collapse
Affiliation(s)
- Shayantan Chaudhuri
- Department
of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
- Centre
for Doctoral Training in Diamond Science and Technology, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Andrew J. Logsdail
- Cardiff
Catalysis Institute, School of Chemistry, Cardiff University, Cardiff CF10 3AT, United
Kingdom
| | - Reinhard J. Maurer
- Department
of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
- Department
of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
9
|
Motlová L, Šnajdr I, Kutil Z, Andris E, Ptáček J, Novotná A, Nováková Z, Havlínová B, Tueckmantel W, Dráberová H, Majer P, Schutkowski M, Kozikowski A, Rulíšek L, Bařinka C. Comprehensive Mechanistic View of the Hydrolysis of Oxadiazole-Based Inhibitors by Histone Deacetylase 6 (HDAC6). ACS Chem Biol 2023. [PMID: 37392419 PMCID: PMC10367051 DOI: 10.1021/acschembio.3c00212] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/03/2023]
Abstract
Histone deacetylase (HDAC) inhibitors used in the clinic typically contain a hydroxamate zinc-binding group (ZBG). However, more recent work has shown that the use of alternative ZBGs, and, in particular, the heterocyclic oxadiazoles, can confer higher isoenzyme selectivity and more favorable ADMET profiles. Herein, we report on the synthesis and biochemical, crystallographic, and computational characterization of a series of oxadiazole-based inhibitors selectively targeting the HDAC6 isoform. Surprisingly, but in line with a very recent finding reported in the literature, a crystal structure of the HDAC6/inhibitor complex revealed that hydrolysis of the oxadiazole ring transforms the parent oxadiazole into an acylhydrazide through a sequence of two hydrolytic steps. An identical cleavage pattern was also observed both in vitro using the purified HDAC6 enzyme as well as in cellular systems. By employing advanced quantum and molecular mechanics (QM/MM) and QM calculations, we elucidated the mechanistic details of the two hydrolytic steps to obtain a comprehensive mechanistic view of the double hydrolysis of the oxadiazole ring. This was achieved by fully characterizing the reaction coordinate, including identification of the structures of all intermediates and transition states, together with calculations of their respective activation (free) energies. In addition, we ruled out several (intuitively) competing pathways. The computed data (ΔG‡ ≈ 21 kcal·mol-1 for the rate-determining step of the overall dual hydrolysis) are in very good agreement with the experimentally determined rate constants, which a posteriori supports the proposed reaction mechanism. We also clearly (and quantitatively) explain the role of the -CF3 or -CHF2 substituent on the oxadiazole ring, which is a prerequisite for hydrolysis to occur. Overall, our data provide compelling evidence that the oxadiazole warheads can be efficiently transformed within the active sites of target metallohydrolases to afford reaction products possessing distinct selectivity and inhibition profiles.
Collapse
Affiliation(s)
- Lucia Motlová
- Institute of Biotechnology of the Czech Academy of Sciences, BIOCEV, Prumyslova 595, 252 50 Vestec, Czech Republic
| | - Ivan Šnajdr
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Prague 6, Czech Republic
| | - Zsófia Kutil
- Institute of Biotechnology of the Czech Academy of Sciences, BIOCEV, Prumyslova 595, 252 50 Vestec, Czech Republic
| | - Erik Andris
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Prague 6, Czech Republic
| | - Jakub Ptáček
- Institute of Biotechnology of the Czech Academy of Sciences, BIOCEV, Prumyslova 595, 252 50 Vestec, Czech Republic
| | - Adéla Novotná
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Prague 6, Czech Republic
| | - Zora Nováková
- Institute of Biotechnology of the Czech Academy of Sciences, BIOCEV, Prumyslova 595, 252 50 Vestec, Czech Republic
| | - Barbora Havlínová
- Institute of Biotechnology of the Czech Academy of Sciences, BIOCEV, Prumyslova 595, 252 50 Vestec, Czech Republic
| | - Werner Tueckmantel
- StarWise Therapeutics LLC, University Research Park, Inc., Madison, Wisconsin 53719, United States
| | - Helena Dráberová
- Institute of Biotechnology of the Czech Academy of Sciences, BIOCEV, Prumyslova 595, 252 50 Vestec, Czech Republic
| | - Pavel Majer
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Prague 6, Czech Republic
| | - Mike Schutkowski
- Department of Enzymology, Charles Tanford Protein Center, Institute of Biochemistry and Biotechnology, Martin-Luther-University Halle-Wittenberg, 06120 Halle, Germany
| | - Alan Kozikowski
- StarWise Therapeutics LLC, University Research Park, Inc., Madison, Wisconsin 53719, United States
| | - Lubomír Rulíšek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Prague 6, Czech Republic
| | - Cyril Bařinka
- Institute of Biotechnology of the Czech Academy of Sciences, BIOCEV, Prumyslova 595, 252 50 Vestec, Czech Republic
| |
Collapse
|
10
|
Chow M, Lambros E, Li X, Hammes-Schiffer S. Nuclear-Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics. J Chem Theory Comput 2023. [PMID: 37329317 DOI: 10.1021/acs.jctc.3c00361] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear-electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. These developments and initial examples provide the foundation for future studies of nuclear-electronic quantum dynamics in complex chemical and biological environments.
Collapse
Affiliation(s)
- Mathew Chow
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Eleftherios Lambros
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | | |
Collapse
|
11
|
Shirazi J, Jafari S, Ryde U, Irani M. Catalytic Reaction Mechanism of Glyoxalase II: A Quantum Mechanics/Molecular Mechanics Study. J Phys Chem B 2023; 127:4480-4495. [PMID: 37191640 DOI: 10.1021/acs.jpcb.3c01495] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
Methylglyoxal (MG) is a reactive and toxic compound produced in carbohydrate, lipid, and amino acid metabolism. The glyoxalase system is the main detoxifying route for MG and consists of two enzymes, glyoxalase I (GlxI) and glyoxalase II (GlxII). GlxI catalyzes the formation of S-d-lactoylglutathione from hemithioacetal, and GlxII converts this intermediate to d-lactate. A relationship between the glyoxalase system and some diseases like diabetes has been shown, and inhibiting enzymes of this system may be an effective means of controlling certain diseases. A detailed understanding of the reaction mechanism of an enzyme is essential to the rational design of competitive inhibitors. In this work, we use quantum mechanics/molecular mechanics (QM/MM) calculations and energy refinement utilizing the big-QM and QM/MM thermodynamic cycle perturbation methods to propose a mechanism for the GlxII reaction that starts with a nucleophilic attack of the bridging OH- group on the substrate. The coordination of the substrate to the Zn ions places its electrophilic center close to the hydroxide group, enabling the reaction to proceed. Our estimated reaction energies are in excellent agreement with experimental data, thus demonstrating the reliability of our approach and the proposed mechanism. Additionally, we examined alternative protonation states of Asp-29, Asp-58, Asp-134, and the bridging hydroxide ion in the catalytic process. However, these give less favorable reactions, a poorer reproduction of the crystal structure geometry of the active site, and higher root-mean-squared deviations of the active site residues in molecular dynamics simulations.
Collapse
Affiliation(s)
- Javad Shirazi
- Department of Chemistry, University of Kurdistan, P.O. Box 66175-416, 66177-15177 Sanandaj, Iran
| | - Sonia Jafari
- Department of Chemistry, University of Kurdistan, P.O. Box 66175-416, 66177-15177 Sanandaj, Iran
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Mehdi Irani
- Department of Chemistry, University of Kurdistan, P.O. Box 66175-416, 66177-15177 Sanandaj, Iran
| |
Collapse
|
12
|
Torbjörnsson M, Hagemann MM, Ryde U, Hedegård ED. Histidine oxidation in lytic polysaccharide monooxygenase. J Biol Inorg Chem 2023; 28:317-328. [PMID: 36828975 PMCID: PMC10036459 DOI: 10.1007/s00775-023-01993-4] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2022] [Accepted: 01/10/2023] [Indexed: 02/26/2023]
Abstract
The lytic polysaccharide monooxygenases (LPMOs) comprise a super-family of copper enzymes that boost the depolymerisation of polysaccharides by oxidatively disrupting the glycosidic bonds connecting the sugar units. Industrial use of LPMOs for cellulose depolymerisation has already begun but is still far from reaching its full potential. One issue is that the LPMOs self-oxidise and thereby deactivate. The mechanism of this self-oxidation is unknown, but histidine residues coordinating to the copper atom are the most susceptible. An unusual methyl modification of the NE2 atom in one of the coordinating histidine residues has been proposed to have a protective role. Furthermore, substrate binding is also known to reduce oxidative damage. We here for the first time investigate the mechanism of histidine oxidation with combined quantum and molecular mechanical (QM/MM) calculations, with outset in intermediates previously shown to form from a reaction with peroxide and a reduced LPMO. We show that an intermediate with a [Cu-O]+ moiety is sufficiently potent to oxidise the nearest C-H bond on both histidine residues, but methylation of the NE2 atom of His-1 increases the reaction barrier of this reaction. The substrate further increases the activation barrier. We also investigate a [Cu-OH]2+ intermediate with a deprotonated tyrosine radical. This intermediate was previously proposed to have a protective role, and we also find it to have higher barriers than the corresponding a [Cu-O]+ intermediate.
Collapse
Affiliation(s)
- Magne Torbjörnsson
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Marlisa M Hagemann
- Department of Physics, Chemistry, and Pharmacy, University of Southern Denmark, Campusvej 55, 5230, Odense M, Denmark
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden.
| | - Erik Donovan Hedegård
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden.
- Department of Physics, Chemistry, and Pharmacy, University of Southern Denmark, Campusvej 55, 5230, Odense M, Denmark.
| |
Collapse
|
13
|
Bramley GA, Beynon OT, Stishenko PV, Logsdail AJ. The application of QM/MM simulations in heterogeneous catalysis. Phys Chem Chem Phys 2023; 25:6562-6585. [PMID: 36810655 DOI: 10.1039/d2cp04537k] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
Abstract
The QM/MM simulation method is provenly efficient for the simulation of biological systems, where an interplay of extensive environment and delicate local interactions drives a process of interest through a funnel on a complex energy landscape. Recent advances in quantum chemistry and force-field methods present opportunities for the adoption of QM/MM to simulate heterogeneous catalytic processes, and their related systems, where similar intricacies exist on the energy landscape. Herein, the fundamental theoretical considerations for performing QM/MM simulations, and the practical considerations for setting up QM/MM simulations of catalytic systems, are introduced; then, areas of heterogeneous catalysis are explored where QM/MM methods have been most fruitfully applied. The discussion includes simulations performed for adsorption processes in solvent at metallic interfaces, reaction mechanisms within zeolitic systems, nanoparticles, and defect chemistry within ionic solids. We conclude with a perspective on the current state of the field and areas where future opportunities for development and application exist.
Collapse
Affiliation(s)
- Gabriel Adrian Bramley
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, CF10 3AT, UK.
| | - Owain Tomos Beynon
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, CF10 3AT, UK.
| | | | - Andrew James Logsdail
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, CF10 3AT, UK.
| |
Collapse
|
14
|
Izsák R, Riplinger C, Blunt NS, de Souza B, Holzmann N, Crawford O, Camps J, Neese F, Schopf P. Quantum computing in pharma: A multilayer embedding approach for near future applications. J Comput Chem 2023; 44:406-421. [PMID: 35789492 DOI: 10.1002/jcc.26958] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2022] [Revised: 06/09/2022] [Accepted: 06/10/2022] [Indexed: 01/03/2023]
Abstract
Quantum computers are special purpose machines that are expected to be particularly useful in simulating strongly correlated chemical systems. The quantum computer excels at treating a moderate number of orbitals within an active space in a fully quantum mechanical manner. We present a quantum phase estimation calculation on F2 in a (2,2) active space on Rigetti's Aspen-11 QPU. While this is a promising start, it also underlines the need for carefully selecting the orbital spaces treated by the quantum computer. In this work, a scheme for selecting such an active space automatically is described and simulated results obtained using both the quantum phase estimation (QPE) and variational quantum eigensolver (VQE) algorithms are presented and combined with a subtractive method to enable accurate description of the environment. The active occupied space is selected from orbitals localized on the chemically relevant fragment of the molecule, while the corresponding virtual space is chosen based on the magnitude of interactions with the occupied space calculated from perturbation theory. This protocol is then applied to two chemical systems of pharmaceutical relevance: the enzyme [Fe] hydrogenase and the photosenzitizer temoporfin. While the sizes of the active spaces currently amenable to a quantum computational treatment are not enough to demonstrate quantum advantage, the procedure outlined here is applicable to any active space size, including those that are outside the reach of classical computation.
Collapse
Affiliation(s)
| | | | | | | | - Nicole Holzmann
- Riverlane Research Ltd, Cambridge, UK.,Astex Pharmaceuticals, Cambridge, UK
| | | | | | - Frank Neese
- Max-Planck Institut für Kohlenforschung, Mülheim an der Ruhr, Germany
| | | |
Collapse
|
15
|
Pérez-Barcia Á, Cárdenas G, Nogueira JJ, Mandado M. Effect of the QM Size, Basis Set, and Polarization on QM/MM Interaction Energy Decomposition Analysis. J Chem Inf Model 2023; 63:882-897. [PMID: 36661314 PMCID: PMC9930123 DOI: 10.1021/acs.jcim.2c01184] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
Herein, an Energy Decomposition Analysis (EDA) scheme extended to the framework of QM/MM calculations in the context of electrostatic embeddings (QM/MM-EDA) including atomic charges and dipoles is applied to assess the effect of the QM region size on the convergence of the different interaction energy components, namely, electrostatic, Pauli, and polarization, for cationic, anionic, and neutral systems interacting with a strong polar environment (water). Significant improvements are found when the bulk solvent environment is described by a MM potential in the EDA scheme as compared to pure QM calculations that neglect bulk solvation. The predominant electrostatic interaction requires sizable QM regions. The results reported here show that it is necessary to include a surprisingly large number of water molecules in the QM region to obtain converged values for this energy term, contrary to most cluster models often employed in the literature. Both the improvement of the QM wave function by means of a larger basis set and the introduction of polarization into the MM region through a polarizable force field do not translate to a faster convergence with the QM region size, but they lead to better results for the different interaction energy components. The results obtained in this work provide insight into the effect of each energy component on the convergence of the solute-solvent interaction energy with the QM region size. This information can be used to improve the MM FFs and embedding schemes employed in QM/MM calculations of solvated systems.
Collapse
Affiliation(s)
- Álvaro Pérez-Barcia
- Department
of Physical Chemistry, University of Vigo, Lagoas-Marcosende s\n, ES-36310-Vigo, Galicia, Spain
| | - Gustavo Cárdenas
- Department
of Chemistry, Universidad Autónoma
de Madrid, 28049, Madrid, Spain
| | - Juan J. Nogueira
- Department
of Chemistry, Universidad Autónoma
de Madrid, 28049, Madrid, Spain,Institute
for Advanced Research in Chemistry (IAdChem), Universidad Autónoma de Madrid, 28049Madrid, Spain,E-mail:
| | - Marcos Mandado
- Department
of Physical Chemistry, University of Vigo, Lagoas-Marcosende s\n, ES-36310-Vigo, Galicia, Spain,E-mail:
| |
Collapse
|
16
|
Sharma H, Raju B, Narendra G, Motiwale M, Sharma B, Verma H, Silakari O. QM/MM Studies on Enzyme Catalysis and Insight into Designing of New Inhibitors by ONIOM Approach: Recent Update. ChemistrySelect 2023. [DOI: 10.1002/slct.202203319] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Himani Sharma
- Molecular Modeling Lab (MML) Department of Pharmaceutical Sciences and Drug Research Punjabi University Patiala Punjab 147002 India
| | - Baddipadige Raju
- Molecular Modeling Lab (MML) Department of Pharmaceutical Sciences and Drug Research Punjabi University Patiala Punjab 147002 India
| | - Gera Narendra
- Molecular Modeling Lab (MML) Department of Pharmaceutical Sciences and Drug Research Punjabi University Patiala Punjab 147002 India
| | - Mohit Motiwale
- Molecular Modeling Lab (MML) Department of Pharmaceutical Sciences and Drug Research Punjabi University Patiala Punjab 147002 India
| | - Bhavna Sharma
- Molecular Modeling Lab (MML) Department of Pharmaceutical Sciences and Drug Research Punjabi University Patiala Punjab 147002 India
| | - Himanshu Verma
- Molecular Modeling Lab (MML) Department of Pharmaceutical Sciences and Drug Research Punjabi University Patiala Punjab 147002 India
| | - Om Silakari
- Molecular Modeling Lab (MML) Department of Pharmaceutical Sciences and Drug Research Punjabi University Patiala Punjab 147002 India
| |
Collapse
|
17
|
Kar RK. Benefits of hybrid QM/MM over traditional classical mechanics in pharmaceutical systems. Drug Discov Today 2023; 28:103374. [PMID: 36174967 DOI: 10.1016/j.drudis.2022.103374] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2022] [Revised: 06/27/2022] [Accepted: 09/22/2022] [Indexed: 02/02/2023]
Abstract
Hybrid quantum mechanics/molecular mechanics (QM/MM) is one of the most reliable approaches for accurately modeling and studying the complex pharmaceutical discovery system. Classical mechanics has significantly accelerated the drug discovery process in the past decade. However, the current challenge is the large pool of false positives, which require extensive validation. Hybrid QM/MM is an effective solution for accurately studying ligand binding, structural mechanisms, free energy evaluation, and spectroscopic characterization. This article highlights the methodological details relevant to cost-effective hybrid QM/MM methods. This approach, combined with traditional pharmacoinformatics methods, could be a reliable strategy to balance the cost and accuracy of the calculations.
Collapse
Affiliation(s)
- Rajiv K Kar
- Jyoti and Bhupat Mehta School of Health Sciences and Technology, Indian Institute of Technology Guwahati, Guwahati, Assam 781039, India.
| |
Collapse
|
18
|
Quantum Mechanical Calculations of Redox Potentials of the Metal Clusters in Nitrogenase. MOLECULES (BASEL, SWITZERLAND) 2022; 28:molecules28010065. [PMID: 36615260 PMCID: PMC9822455 DOI: 10.3390/molecules28010065] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/18/2022] [Revised: 12/14/2022] [Accepted: 12/19/2022] [Indexed: 12/24/2022]
Abstract
We have calculated redox potentials of the two metal clusters in Mo-nitrogenase with quantum mechanical (QM) calculations. We employ an approach calibrated for iron-sulfur clusters with 1-4 Fe ions, involving QM-cluster calculations in continuum solvent and large QM systems (400-500 atoms), based on structures from combined QM and molecular mechanics (QM/MM) geometry optimisations. Calculations on the P-cluster show that we can reproduce the experimental redox potentials within 0.33 V. This is similar to the accuracy obtained for the smaller clusters, although two of the redox reactions involve also proton transfer. The calculated P1+/PN redox potential is nearly the same independently of whether P1+ is protonated or deprotonated, explaining why redox titrations do not show any pH dependence. For the FeMo cluster, the calculations clearly show that the formal oxidation state of the cluster in the resting E0 state is MoIIIFe3IIFe4III , in agreement with previous experimental studies and QM calculations. Moreover, the redox potentials of the first five E0-E4 states are nearly constant, as is expected if the electrons are delivered by the same site (the P-cluster). However, the redox potentials are insensitive to the formal oxidation states of the Fe ion (i.e., whether the added protons bind to sulfide or Fe ions). Finally, we show that the later (E4-E8) states of the reaction mechanism have redox potential that are more positive (i.e., more exothermic) than that of the E0/E1 couple.
Collapse
|
19
|
Blunt NS, Camps J, Crawford O, Izsák R, Leontica S, Mirani A, Moylett AE, Scivier SA, Sünderhauf C, Schopf P, Taylor JM, Holzmann N. Perspective on the Current State-of-the-Art of Quantum Computing for Drug Discovery Applications. J Chem Theory Comput 2022; 18:7001-7023. [PMID: 36355616 DOI: 10.1021/acs.jctc.2c00574] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Computational chemistry is an essential tool in the pharmaceutical industry. Quantum computing is a fast evolving technology that promises to completely shift the computational capabilities in many areas of chemical research by bringing into reach currently impossible calculations. This perspective illustrates the near-future applicability of quantum computation of molecules to pharmaceutical problems. We briefly summarize and compare the scaling properties of state-of-the-art quantum algorithms and provide novel estimates of the quantum computational cost of simulating progressively larger embedding regions of a pharmaceutically relevant covalent protein-drug complex involving the drug Ibrutinib. Carrying out these calculations requires an error-corrected quantum architecture that we describe. Our estimates showcase that recent developments on quantum phase estimation algorithms have dramatically reduced the quantum resources needed to run fully quantum calculations in active spaces of around 50 orbitals and electrons, from estimated over 1000 years using the Trotterization approach to just a few days with sparse qubitization, painting a picture of fast and exciting progress in this nascent field.
Collapse
Affiliation(s)
- Nick S Blunt
- Riverlane, St. Andrews House, 59 St. Andrews Street, Cambridge CB2 3BZ, United Kingdom
| | - Joan Camps
- Riverlane, St. Andrews House, 59 St. Andrews Street, Cambridge CB2 3BZ, United Kingdom
| | - Ophelia Crawford
- Riverlane, St. Andrews House, 59 St. Andrews Street, Cambridge CB2 3BZ, United Kingdom
| | - Róbert Izsák
- Riverlane, St. Andrews House, 59 St. Andrews Street, Cambridge CB2 3BZ, United Kingdom
| | - Sebastian Leontica
- Riverlane, St. Andrews House, 59 St. Andrews Street, Cambridge CB2 3BZ, United Kingdom
| | - Arjun Mirani
- Riverlane, St. Andrews House, 59 St. Andrews Street, Cambridge CB2 3BZ, United Kingdom
| | - Alexandra E Moylett
- Riverlane, St. Andrews House, 59 St. Andrews Street, Cambridge CB2 3BZ, United Kingdom
| | - Sam A Scivier
- Riverlane, St. Andrews House, 59 St. Andrews Street, Cambridge CB2 3BZ, United Kingdom
| | - Christoph Sünderhauf
- Riverlane, St. Andrews House, 59 St. Andrews Street, Cambridge CB2 3BZ, United Kingdom
| | - Patrick Schopf
- Astex Pharmaceuticals, 436 Cambridge Science Park, Cambridge CB4 0QA, United Kingdom
| | - Jacob M Taylor
- Riverlane, St. Andrews House, 59 St. Andrews Street, Cambridge CB2 3BZ, United Kingdom
| | - Nicole Holzmann
- Riverlane, St. Andrews House, 59 St. Andrews Street, Cambridge CB2 3BZ, United Kingdom.,Astex Pharmaceuticals, 436 Cambridge Science Park, Cambridge CB4 0QA, United Kingdom
| |
Collapse
|
20
|
Jiang H, Svensson OKG, Ryde U. QM/MM Study of Partial Dissociation of S2B for the E 2 Intermediate of Nitrogenase. Inorg Chem 2022; 61:18067-18076. [DOI: 10.1021/acs.inorgchem.2c02488] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Hao Jiang
- Department of Theoretical Chemistry, Lund University, Chemical Centre, SE-221 00Lund, Sweden
| | - Oskar K. G. Svensson
- Department of Theoretical Chemistry, Lund University, Chemical Centre, SE-221 00Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, SE-221 00Lund, Sweden
| |
Collapse
|
21
|
Falbo E, Fusè M, Lazzari F, Mancini G, Barone V. Integration of Quantum Chemistry, Statistical Mechanics, and Artificial Intelligence for Computational Spectroscopy: The UV-Vis Spectrum of TEMPO Radical in Different Solvents. J Chem Theory Comput 2022; 18:6203-6216. [PMID: 36166322 PMCID: PMC9558374 DOI: 10.1021/acs.jctc.2c00654] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
The ongoing integration of quantum chemistry, statistical
mechanics,
and artificial intelligence is paving the route toward more effective
and accurate strategies for the investigation of the spectroscopic
properties of medium-to-large size chromophores in condensed phases.
In this context we are developing a novel workflow aimed at improving
the generality, reliability, and ease of use of the available computational
tools. In this paper we report our latest developments with specific
reference to unsupervised atomistic simulations employing non periodic
boundary conditions (NPBC) followed by clustering of the trajectories
employing optimized feature spaces. Next accurate variational computations
are performed for a representative point of each cluster, whereas
intracluster fluctuations are taken into account by a cheap yet reliable
perturbative approach. A number of methodological improvements have
been introduced including, e.g., more realistic reaction field effects
at the outer boundary of the simulation sphere, automatic definition
of the feature space by continuous perception of solute–solvent
interactions, full account of polarization and charge transfer in
the first solvation shell, and inclusion of vibronic contributions.
After its validation, this new approach has been applied to the challenging
case of solvatochromic effects on the UV–vis spectra of a prototypical
nitroxide radical (TEMPO) in different solvents. The reliability,
effectiveness, and robustness of the new platform is demonstrated
by the remarkable agreement with experiment of the results obtained
through an unsupervised approach characterized by a strongly reduced
computational cost as compared to that of conventional quantum mechanics
and molecular mechanics models without any accuracy reduction.
Collapse
Affiliation(s)
- Emanuele Falbo
- Scuola Normale Superiore di Pisa, piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Marco Fusè
- Scuola Normale Superiore di Pisa, piazza dei Cavalieri 7, 56126 Pisa, Italy.,Dipartimento di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123 Brescia, Italy
| | - Federico Lazzari
- Scuola Normale Superiore di Pisa, piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Giordano Mancini
- Scuola Normale Superiore di Pisa, piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore di Pisa, piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
22
|
Lawniczak JJ, Zhang X, Christianson M, Bailey B, Bremer S, Barcia S, Mukhopadhyay S, Klosin J, Miller TF. Solution-Phase Conformational/Vibrational Anharmonicity in Comonomer Incorporation Polyolefin Catalysis. J Phys Chem A 2022; 126:6858-6869. [PMID: 36137217 DOI: 10.1021/acs.jpca.2c04038] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The prediction of comonomer incorporation statistics in polyolefin catalysis necessitates an accurate calculation of free energies corresponding to monomer binding and insertion, often requiring sub-kcal/mol resolution to resolve experimental free energies. Batch reactor experiments are used to probe incorporation statistics of ethene and larger α-olefins for three constrained geometry complexes which are employed as model systems. Herein, over 6 ns of quantum mechanics/molecular mechanics (QM/MM) molecular dynamics is performed in combination with the zero-temperature string method to characterize the solution-phase insertion barrier and to analyze the contributions from conformational and vibrational anharmonicity arising both in vacuum and in solution. Conformational sampling in the solution-phase results in 0-2 kcal/mol corrections to the insertion barrier which are on the same scale necessary to resolve experimental free energies. Anharmonic contributions from conformational sampling in the solution phase are crucial energy contributions missing from static density functional theory calculations and implicit solvation models, and the accurate calculation of these contributions is a key step toward the quantitative prediction of comonomer incorporation statistics.
Collapse
Affiliation(s)
- James J Lawniczak
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Xinglong Zhang
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | | | - Brad Bailey
- The Dow Chemical Company, Midland, Michigan 48674, United States
| | - Sean Bremer
- The Dow Chemical Company, Midland, Michigan 48674, United States
| | - Sarah Barcia
- Kelly Services, Inc., Troy, Michigan 48084, United States
| | | | - Jerzy Klosin
- The Dow Chemical Company, Midland, Michigan 48674, United States
| | - Thomas F Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
23
|
Barone V, Carnimeo I, Mancini G, Pagliai M. Development, Validation, and Pilot Application of a Generalized Fluctuating Charge Model for Computational Spectroscopy in Solution. ACS OMEGA 2022; 7:13382-13394. [PMID: 35474835 PMCID: PMC9026056 DOI: 10.1021/acsomega.2c01110] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/07/2022] [Accepted: 03/23/2022] [Indexed: 06/14/2023]
Abstract
A general approach enforcing nonperiodic boundary conditions for the computation of spectroscopic properties in solution has been improved including an effective description of charge-transfer contributions and coordination number adjustment for explicit solvent molecules. Both contributions are obtained from a continuous description of intermolecular hydrogen bonds, which has been employed also for an effective clustering of molecular dynamics trajectories. Fine tuning of the model has been performed for several water clusters, and then its efficiency and reliability have been demonstrated by computing the absorption spectra of different creatinine tautomers in aqueous solution.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Ivan Carnimeo
- Scuola
Internazionale Superiore di Studi Avanzati, via Bonomea 265, 34136 Trieste, Italy
| | - Giordano Mancini
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Marco Pagliai
- Dipartimento
di Chimica “Ugo Schiff”, Università
degli Studi di Firenze, Via della Lastruccia 3, 50019 Sesto Fiorentino, Italy
| |
Collapse
|
24
|
Pederson JP, McDaniel J. DFT-based QM/MM with Particle-Mesh Ewald for Direct, Long-Range Electrostatic Embedding. J Chem Phys 2022; 156:174105. [DOI: 10.1063/5.0087386] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a DFT-based, QM/MM implementation with long-range electrostatic embedding achieved by direct real-space integration of the particle mesh Ewald (PME) computed electrostatic potential. The key transformation is the interpolation of the electrostatic potential from the PME grid to the DFT quadrature grid, from which integrals are easily evaluated utilizing standard DFT machinery. We provide benchmarks of the numerical accuracy with choice of grid size and real-space corrections, and demonstrate that good convergence is achieved while introducing nominal computational overhead. Furthermore, the approach requires only small modification to existing software packages, as is demonstrated with our implementation in the OpenMM and Psi4 software. After presenting convergence benchmarks, we evaluate the importance of long-range electrostatic embedding in three solute/solvent systems modeled with QM/MM. Water and BMIM/BF4 ionic liquid were considered as ``simple' and ``complex' solvents respectively, with water and p-phenylenediamine (PPD) solute molecules treated at QM level of theory. While electrostatic embedding with standard real-space truncation may introduce negligible error for simple systems such as water solute in water solvent, errors become more significant when QM/MM is applied to complex solvents such as ionic liquids. An extreme example is the electrostatic embedding energy for oxidized PPD in BMIM/BF4 for which real-space truncation produces severe error even at 2-3 nm cutoff distances. This latter example illustrates that utilization of QM/MM to compute redox potentials within concentrated electrolytes/ionic media requires carefully chosen long-range electrostatic embedding algorithms, with our presented algorithm providing a general and robust approach.
Collapse
Affiliation(s)
| | - Jesse McDaniel
- Chemistry, Georgia Institute of Technology, United States of America
| |
Collapse
|
25
|
Jafari S, Tavares Santos YA, Bergmann J, Irani M, Ryde U. Benchmark Study of Redox Potential Calculations for Iron-Sulfur Clusters in Proteins. Inorg Chem 2022; 61:5991-6007. [PMID: 35403427 PMCID: PMC9044450 DOI: 10.1021/acs.inorgchem.1c03422] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Redox potentials
have been calculated for 12 different iron–sulfur
sites of 6 different types with 1–4 iron ions. Structures were
optimized with combined quantum mechanical and molecular mechanical
(QM/MM) methods, and the redox potentials were calculated using the
QM/MM energies, single-point QM methods in a continuum solvent or
by QM/MM thermodynamic cycle perturbations. We show that the best
results are obtained with a large QM system (∼300 atoms, but
a smaller QM system, ∼150 atoms, can be used for the QM/MM
geometry optimization) and a large value of the dielectric constant
(80). For absolute redox potentials, the B3LYP density functional
method gives better results than TPSS, and the results are improved
with a larger basis set. However, for relative redox potentials, the
opposite is true. The results are insensitive to the force field (charges
of the surroundings) used for the QM/MM calculations or whether the
protein and solvent outside the QM system are relaxed or kept fixed
at the crystal structure. With the best approach for relative potentials,
mean absolute and maximum deviations of 0.17 and 0.44 V, respectively,
are obtained after removing a systematic error of −0.55 V.
Such an approach can be used to identify the correct oxidation states
involved in a certain redox reaction. We
have studied redox potentials of 12 iron−sulfur
sites of 6 types with 1−4 iron ions. Structures were optimized
with combined quantum mechanical and molecular mechanical (QM/MM)
methods, and the redox potentials were calculated with QM/MM, QM calculations
in a continuum solvent or by QM/MM thermodynamic cycle perturbations.
The best results are obtained with the second approach using ∼300
atoms in the QM model and a large dielectric constant.
Collapse
Affiliation(s)
- Sonia Jafari
- Department of Chemistry, University of Kurdistan, 66175-416 Sanandaj, Iran.,Department of Theoretical Chemistry, Chemical Centre, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Yakini A Tavares Santos
- Department of Theoretical Chemistry, Chemical Centre, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Justin Bergmann
- Department of Theoretical Chemistry, Chemical Centre, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Mehdi Irani
- Department of Chemistry, University of Kurdistan, 66175-416 Sanandaj, Iran
| | - Ulf Ryde
- Department of Theoretical Chemistry, Chemical Centre, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
26
|
Brandt F, Jacob CR. Systematic QM Region Construction in QM/MM Calculations Based on Uncertainty Quantification. J Chem Theory Comput 2022; 18:2584-2596. [PMID: 35271768 DOI: 10.1021/acs.jctc.1c01093] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
While QM/MM studies of enzymatic reactions are widely used in computational chemistry, the results of such studies are subject to numerous sources of uncertainty, and the effect of different choices by the simulation scientist that are required when setting up QM/MM calculations is often unclear. In particular, the selection of the QM region is crucial for obtaining accurate and reliable results. Simply including amino acids by their distance to the active site is mostly not sufficient as necessary residues are missing or unimportant residues are included without evidence. Here, we take a first step toward quantifying uncertainties in QM/MM calculations by assessing the sensitivity of QM/MM reaction energies with respect to variations of the MM point charges. We show that such a point charge variation analysis (PCVA) can be employed to judge the accuracy of QM/MM reaction energies obtained with a selected QM region and devise a protocol to systematically construct QM regions that minimize this uncertainty. We apply such a PCVA to the example of catechol O-methyltransferase and demonstrate that it provides a simple and reliable approach for the construction of the QM region. Our PCVA-based scheme is computationally efficient and requires only calculations for a system with a minimal QM region. Our work highlights the promise of applying methods of uncertainty quantification in computational chemistry.
Collapse
Affiliation(s)
- Felix Brandt
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstr. 17, 38106 Braunschweig, Germany
| | - Christoph R Jacob
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstr. 17, 38106 Braunschweig, Germany
| |
Collapse
|
27
|
Mancini G, Fusè M, Lipparini F, Nottoli M, Scalmani G, Barone V. Molecular Dynamics Simulations Enforcing Nonperiodic Boundary Conditions: New Developments and Application to the Solvent Shifts of Nitroxide Magnetic Parameters. J Chem Theory Comput 2022; 18:2479-2493. [PMID: 35257572 PMCID: PMC9009096 DOI: 10.1021/acs.jctc.2c00046] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
Abstract
![]()
Multiscale methods
combining quantum mechanics and molecular mechanics
(QM/MM) have become the most suitable and effective strategies for
the investigation of the spectroscopic properties of medium-to-large
size chromophores in condensed phases. In this context, we are developing
a novel workflow aimed at improving the generality, reliability, and
ease of use of the available computational tools. In this paper, we
report our latest developments with specific reference to a general
protocol based on atomistic simulations, carried out under nonperiodic
boundary conditions (NPBC). In particular, we add to our in house
MD engine a new efficient treatment of mean field electrostatic contributions
to energy and forces, together with the capability of performing the
simulations either in the canonical (NVT) or in the
isothermal–isobaric (NPT) ensemble. Next,
we provide convincing evidence that the NBPC approach enhanced by
specific tweaks for rigid body propagation, allows for the simulation
of solute–solvent systems with a minimum number of degrees
of freedom and large integration time step. After its validation,
this new approach is applied to the challenging case of solvatochromic
effects on the electron paramagnetic resonance (EPR) spectrum of a
prototypical nitroxide radical. To this end, we propose and validate
also an automated protocol to extract and weight simulation snapshots,
making use of a continuous description of the strength of solute–solvent
hydrogen bridges. While further developments are being worked on,
the effectiveness of our approach, even in its present form, is demonstrated
by the accuracy of the results obtained through an unsupervised approach
characterized by a strongly reduced computational cost as compared
to that of conventional QM/MM models, without any appreciable deterioration
of accuracy.
Collapse
Affiliation(s)
- Giordano Mancini
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy.,Istituto Nazionale di Fisica Nucleare (INFN) sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
| | - Marco Fusè
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Universitaá di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Universitaá di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Giovanni Scalmani
- Gaussian, Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, United States
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy.,Istituto Nazionale di Fisica Nucleare (INFN) sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
| |
Collapse
|
28
|
Jiang H, Ryde U. Thermodynamically Favourable States in the Reaction of Nitrogenase without Dissociation of any Sulfide Ligand. Chemistry 2022; 28:e202103933. [PMID: 35006641 PMCID: PMC9305431 DOI: 10.1002/chem.202103933] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2021] [Indexed: 12/16/2022]
Abstract
We have used combined quantum mechanical and molecular mechanical (QM/MM) calculations to study the reaction mechanism of nitrogenase, assuming that none of the sulfide ligands dissociates. To avoid the problem that there is no consensus regarding the structure and protonation of the E4 state, we start from a state where N2 is bound to the cluster and is protonated to N2H2, after dissociation of H2. We show that the reaction follows an alternating mechanism with HNNH (possibly protonated to HNNH2) and H2NNH2 as intermediates and the two NH3 products dissociate at the E7 and E8 levels. For all intermediates, coordination to Fe6 is preferred, but for the E4 and E8 intermediates, binding to Fe2 is competitive. For the E4, E5 and E7 intermediates we find that the substrate may abstract a proton from the hydroxy group of the homocitrate ligand of the FeMo cluster, thereby forming HNNH2, H2NNH2 and NH3 intermediates. This may explain why homocitrate is a mandatory component of nitrogenase. All steps in the suggested reaction mechanism are thermodynamically favourable compared to protonation of the nearby His‐195 group and in all cases, protonation of the NE2 atom of the latter group is preferred.
Collapse
Affiliation(s)
- Hao Jiang
- Department of Theoretical Chemistry, Lund University Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| |
Collapse
|
29
|
Rovaletti A, Moro G, Cosentino U, Ryde U, Greco C. Can water act as a nucleophile in CO oxidation catalysed by Mo/Cu CO-dehydrogenase? Answers from theory. Chemphyschem 2022; 23:e202200053. [PMID: 35170169 PMCID: PMC9310835 DOI: 10.1002/cphc.202200053] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2022] [Revised: 02/14/2022] [Indexed: 11/14/2022]
Abstract
The aerobic CO dehydrogenase from Oligotropha carboxidovorans is an environmentally crucial bacterial enzyme for maintenance of subtoxic concentration of CO in the lower atmosphere, as it allows for the oxidation of CO to CO2 which takes place at its Mo−Cu heterobimetallic active site. Despite extensive experimental and theoretical efforts, significant uncertainties still concern the reaction mechanism for the CO oxidation. In this work, we used the hybrid quantum mechanical/molecular mechanical approach to evaluate whether a water molecule present in the active site might act as a nucleophile upon formation of the new C−O bond, a hypothesis recently suggested in the literature. Our study shows that activation of H2O can be favoured by the presence of the Mo=Oeq group. However, overall our results suggest that mechanisms other than the nucleophilic attack by Mo=Oeq to the activated carbon of the CO substrate are not likely to constitute reactive channels for the oxidation of CO by the enzyme.
Collapse
Affiliation(s)
- Anna Rovaletti
- University of Milano-Bicocca: Universita degli Studi di Milano-Bicocca, Department of Earth and Environmental Sciences, ITALY
| | - Giorgio Moro
- University of Milano-Bicocca: Universita degli Studi di Milano-Bicocca, Department of Biotechnology and Biosciences, ITALY
| | - Ugo Cosentino
- University of Milano-Bicocca: Universita degli Studi di Milano-Bicocca, Department of Earth and Environmental Sciences, ITALY
| | - Ulf Ryde
- Lund University: Lunds Universitet, Department of Theoretical Chemistry, ITALY
| | - Claudio Greco
- Università degli Studi di Milano-Bicocca: Universita degli Studi di Milano-Bicocca, earth and environmental sciences, Piazza della Scienza 1, 20126, Milan, ITALY
| |
Collapse
|
30
|
González-Espinoza CE, Rumble CA, Borgis D, Wesolowski TA. Quantifying Fluctuations of Average Solvent Environments for Embedding Calculations. J Chem Theory Comput 2022; 18:1072-1088. [PMID: 35044168 DOI: 10.1021/acs.jctc.1c01108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The viability and effectiveness of replacing an ensemble of embedded solute calculations by a single calculation using an average description of the solvent environment are evaluated. This work explores the fluctuations of the average description of the system obtained in two ways: from calculations on an ensemble of geometries and from an average environment constructed from the same ensemble. To this end, classical molecular dynamics simulations of a rigid acetone solute in SPCE water are performed in order to generate an ensemble of solvent environments. From this ensemble of solvent configurations, a number of different approaches for constructing an average solvent environment are employed. We perform a thorough numerical analysis of the fluctuations of the electrostatic potential experienced by the solute, as well as the resulting fluctuations of the solute's electronic density, measured through its dipole moment and fitted atomic point charges. At the same time, we inspect the accuracy of the methods used to construct average environments. Finally, the proposed method for generating the embedding potential from an average environment density is applied to estimate the solvatochromic shift of the first excitation of acetone. In order to account for quantum confinement effects, which may be important in certain cases, the fluctuations in the shift due to the interaction with the solvent are evaluated using frozen-density-embedding theory. Our results demonstrate that, for normally distributed environments, the constructed average environment is a reasonably good representation of a fluctuating molecular solvent environment. We then provide guidance for future comparisons between these theoretical treatments of solute/solvent systems to experimental measurements.
Collapse
Affiliation(s)
| | - Christopher A Rumble
- Université de Genève, Départment de Chimie Physique 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - Daniel Borgis
- Ecole Normale Supérieure, Départment de Chimie, UMR 8640, PSL University, Sorbonne University, CNRS, 75005 Paris, France
| | - Tomasz A Wesolowski
- Université de Genève, Départment de Chimie Physique 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| |
Collapse
|
31
|
Bím D, Navrátil M, Gutten O, Konvalinka J, Kutil Z, Culka M, Navrátil V, Alexandrova AN, Bařinka C, Rulíšek L. Predicting Effects of Site-Directed Mutagenesis on Enzyme Kinetics by QM/MM and QM Calculations: A Case of Glutamate Carboxypeptidase II. J Phys Chem B 2022; 126:132-143. [PMID: 34978450 DOI: 10.1021/acs.jpcb.1c09240] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Quantum and molecular mechanics (QM/MM) and QM-only (cluster model) modeling techniques represent the two workhorses in mechanistic understanding of enzyme catalysis. One of the stringent tests for QM/MM and/or QM approaches is to provide quantitative answers to real-world biochemical questions, such as the effect of single-point mutations on enzyme kinetics. This translates into predicting the relative activation energies to 1-2 kcal·mol-1 accuracy; such predictions can be used for the rational design of novel enzyme variants with desired/improved characteristics. Herein, we employ glutamate carboxypeptidase II (GCPII), a dizinc metallopeptidase, also known as the prostate specific membrane antigen, as a model system. The structure and activity of this major cancer antigen have been thoroughly studied, both experimentally and computationally, which makes it an ideal model system for method development. Its reaction mechanism is quite well understood: the reaction coordinate comprises a "tetrahedral intermediate" and two transition states and experimental activation Gibbs free energy of ∼17.5 kcal·mol-1 can be inferred for the known kcat ≈ 1 s-1. We correlate experimental kinetic data (including the E424H variant, newly characterized in this work) for various GCPII mutants (kcat = 8.6 × 10-5 s-1 to 2.7 s-1) with the energy profiles calculated by QM/MM and QM-only (cluster model) approaches. We show that the near-quantitative agreement between the experimental values and the calculated activation energies (ΔH⧧) can be obtained and recommend the combination of the two protocols: QM/MM optimized structures and cluster model (QM) energetics. The trend in relative activation energies is mostly independent of the QM method (DFT functional) used. Last but not least, a satisfactory correlation between experimental and theoretical data allows us to provide qualitative and fairly simple explanations of the observed kinetic effects which are thus based on a rigorous footing.
Collapse
Affiliation(s)
- Daniel Bím
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic.,Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, United States
| | - Michal Navrátil
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Ondrej Gutten
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Jan Konvalinka
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic.,Department of Biochemistry, Faculty of Science, Charles University, Hlavova 2030, 2120 00 Prague, Czech Republic
| | - Zsófia Kutil
- Institute of Biotechnology of the Czech Academy of Sciences, Průmyslová 595, 252 50 Vestec, Czech Republic
| | - Martin Culka
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Václav Navrátil
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Anastassia N Alexandrova
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569, United States
| | - Cyril Bařinka
- Institute of Biotechnology of the Czech Academy of Sciences, Průmyslová 595, 252 50 Vestec, Czech Republic
| | - Lubomír Rulíšek
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| |
Collapse
|
32
|
Rossano‐Tapia M, Brown A. Quantum mechanical/molecular mechanical studies of photophysical properties of fluorescent proteins. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1557] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
| | - Alex Brown
- Department of Chemistry University of Alberta Edmonton Alberta Canada
| |
Collapse
|
33
|
Soto-Ospina A, Araque Marín P, Bedoya GDJ, Villegas Lanau A. Structural Predictive Model of Presenilin-2 Protein and Analysis of Structural Effects of Familial Alzheimer's Disease Mutations. Biochem Res Int 2021; 2021:9542038. [PMID: 34881055 PMCID: PMC8648483 DOI: 10.1155/2021/9542038] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2021] [Accepted: 10/21/2021] [Indexed: 12/12/2022] Open
Abstract
Alzheimer's disease manifests itself in brain tissue by neuronal death, due to aggregation of β-amyloid, produced by senile plaques, and hyperphosphorylation of the tau protein, which produces neurofibrillary tangles. One of the genetic markers of the disease is the gene that translates the presenilin-2 protein, which has mutations that favor the appearance of the disease and has no reported crystallographic structure. In view of this, protein modeling is performed using prediction and structural refinement tools followed by an energetic and stereochemical characterization for its validation. For the simulation, four reported mutations are chosen, which are Met239Ile, Met239Val, Ser130Leu, and Thr122Arg, all associated with various functional responses. From a theoretical analysis, a preliminary bioinformatic study is made to find the phosphorylation patterns in the protein and the hydropathic index according to the polarity and chemical environment. Molecular visualization was carried out with the Chimera 1.14 software, and the theoretical calculation with the hybrid quantum mechanics/molecular mechanics system from the semi-empirical method, with Spartan18 software and an AustinModel1 basis. These relationships allow for studying the system from a structural approach with the determination of small distance changes, potential surfaces, electrostatic maps, and angle changes, which favor the comparison between wild-type and mutant systems. With the results obtained, it is expected to complement experimental data reported in the literature from models that would allow us to understand the effects of the selected mutations.
Collapse
Affiliation(s)
- Alejandro Soto-Ospina
- University of Antioquia, Faculty of Medicine, Group Molecular Genetics, Medellín, Colombia
- University of Antioquia, Faculty of Medicine, Group Neuroscience of Antioquia, Medellín, Colombia
| | - Pedronel Araque Marín
- EIA University, School of Life Sciences, Research and Innovation in Chemistry Formulations Group, Envigado, Colombia
| | | | - Andrés Villegas Lanau
- University of Antioquia, Faculty of Medicine, Group Molecular Genetics, Medellín, Colombia
- University of Antioquia, Faculty of Medicine, Group Neuroscience of Antioquia, Medellín, Colombia
| |
Collapse
|
34
|
Epifanovsky E, Gilbert ATB, Feng X, Lee J, Mao Y, Mardirossian N, Pokhilko P, White AF, Coons MP, Dempwolff AL, Gan Z, Hait D, Horn PR, Jacobson LD, Kaliman I, Kussmann J, Lange AW, Lao KU, Levine DS, Liu J, McKenzie SC, Morrison AF, Nanda KD, Plasser F, Rehn DR, Vidal ML, You ZQ, Zhu Y, Alam B, Albrecht BJ, Aldossary A, Alguire E, Andersen JH, Athavale V, Barton D, Begam K, Behn A, Bellonzi N, Bernard YA, Berquist EJ, Burton HGA, Carreras A, Carter-Fenk K, Chakraborty R, Chien AD, Closser KD, Cofer-Shabica V, Dasgupta S, de Wergifosse M, Deng J, Diedenhofen M, Do H, Ehlert S, Fang PT, Fatehi S, Feng Q, Friedhoff T, Gayvert J, Ge Q, Gidofalvi G, Goldey M, Gomes J, González-Espinoza CE, Gulania S, Gunina AO, Hanson-Heine MWD, Harbach PHP, Hauser A, Herbst MF, Hernández Vera M, Hodecker M, Holden ZC, Houck S, Huang X, Hui K, Huynh BC, Ivanov M, Jász Á, Ji H, Jiang H, Kaduk B, Kähler S, Khistyaev K, Kim J, Kis G, Klunzinger P, Koczor-Benda Z, Koh JH, Kosenkov D, Koulias L, Kowalczyk T, Krauter CM, Kue K, Kunitsa A, Kus T, Ladjánszki I, Landau A, Lawler KV, Lefrancois D, Lehtola S, Li RR, Li YP, Liang J, Liebenthal M, Lin HH, Lin YS, Liu F, Liu KY, Loipersberger M, Luenser A, Manjanath A, Manohar P, Mansoor E, Manzer SF, Mao SP, Marenich AV, Markovich T, Mason S, Maurer SA, McLaughlin PF, Menger MFSJ, Mewes JM, Mewes SA, Morgante P, Mullinax JW, Oosterbaan KJ, Paran G, Paul AC, Paul SK, Pavošević F, Pei Z, Prager S, Proynov EI, Rák Á, Ramos-Cordoba E, Rana B, Rask AE, Rettig A, Richard RM, Rob F, Rossomme E, Scheele T, Scheurer M, Schneider M, Sergueev N, Sharada SM, Skomorowski W, Small DW, Stein CJ, Su YC, Sundstrom EJ, Tao Z, Thirman J, Tornai GJ, Tsuchimochi T, Tubman NM, Veccham SP, Vydrov O, Wenzel J, Witte J, Yamada A, Yao K, Yeganeh S, Yost SR, Zech A, Zhang IY, Zhang X, Zhang Y, Zuev D, Aspuru-Guzik A, Bell AT, Besley NA, Bravaya KB, Brooks BR, Casanova D, Chai JD, Coriani S, Cramer CJ, Cserey G, DePrince AE, DiStasio RA, Dreuw A, Dunietz BD, Furlani TR, Goddard WA, Hammes-Schiffer S, Head-Gordon T, Hehre WJ, Hsu CP, Jagau TC, Jung Y, Klamt A, Kong J, Lambrecht DS, Liang W, Mayhall NJ, McCurdy CW, Neaton JB, Ochsenfeld C, Parkhill JA, Peverati R, Rassolov VA, Shao Y, Slipchenko LV, Stauch T, Steele RP, Subotnik JE, Thom AJW, Tkatchenko A, Truhlar DG, Van Voorhis T, Wesolowski TA, Whaley KB, Woodcock HL, Zimmerman PM, Faraji S, Gill PMW, Head-Gordon M, Herbert JM, Krylov AI. Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package. J Chem Phys 2021; 155:084801. [PMID: 34470363 PMCID: PMC9984241 DOI: 10.1063/5.0055522] [Citation(s) in RCA: 451] [Impact Index Per Article: 150.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Abstract
This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.
Collapse
Affiliation(s)
- Evgeny Epifanovsky
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | | | | | - Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Yuezhi Mao
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Pavel Pokhilko
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Alec F. White
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Marc P. Coons
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Adrian L. Dempwolff
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Zhengting Gan
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Diptarka Hait
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Paul R. Horn
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Leif D. Jacobson
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | | | - Jörg Kussmann
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Adrian W. Lange
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Ka Un Lao
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Daniel S. Levine
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Simon C. McKenzie
- Research School of Chemistry, Australian National University, Canberra, Australia
| | | | - Kaushik D. Nanda
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Dirk R. Rehn
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Marta L. Vidal
- Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kgs Lyngby, Denmark
| | | | - Ying Zhu
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Bushra Alam
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Benjamin J. Albrecht
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | | | - Ethan Alguire
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Josefine H. Andersen
- Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kgs Lyngby, Denmark
| | - Vishikh Athavale
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Dennis Barton
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Khadiza Begam
- Department of Physics, Kent State University, Kent, Ohio 44242, USA
| | - Andrew Behn
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Nicole Bellonzi
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Yves A. Bernard
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Hugh G. A. Burton
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Abel Carreras
- Donostia International Physics Center, 20080 Donostia, Euskadi, Spain
| | - Kevin Carter-Fenk
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | | | - Alan D. Chien
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | | | - Vale Cofer-Shabica
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Saswata Dasgupta
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Marc de Wergifosse
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Jia Deng
- Research School of Chemistry, Australian National University, Canberra, Australia
| | | | - Hainam Do
- School of Chemistry, University of Nottingham, Nottingham, United Kingdom
| | - Sebastian Ehlert
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Beringstr. 4, 53115 Bonn, Germany
| | - Po-Tung Fang
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | | | - Qingguo Feng
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Triet Friedhoff
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - James Gayvert
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - Qinghui Ge
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Gergely Gidofalvi
- Department of Chemistry and Biochemistry, Gonzaga University, Spokane, Washington 99258, USA
| | - Matthew Goldey
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Joe Gomes
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Sahil Gulania
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Anastasia O. Gunina
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Phillip H. P. Harbach
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Andreas Hauser
- Institute of Experimental Physics, Graz University of Technology, Graz, Austria
| | | | - Mario Hernández Vera
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Manuel Hodecker
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Zachary C. Holden
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Shannon Houck
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Xunkun Huang
- Department of Chemistry, Xiamen University, Xiamen 361005, China
| | - Kerwin Hui
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Bang C. Huynh
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Maxim Ivanov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Ádám Jász
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | - Hyunjun Ji
- Graduate School of Energy, Environment, Water and Sustainability (EEWS), Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Hanjie Jiang
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Benjamin Kaduk
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Sven Kähler
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Kirill Khistyaev
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Jaehoon Kim
- Graduate School of Energy, Environment, Water and Sustainability (EEWS), Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Gergely Kis
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | | | - Zsuzsanna Koczor-Benda
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Joong Hoon Koh
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Dimitri Kosenkov
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - Laura Koulias
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | | | - Caroline M. Krauter
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Karl Kue
- Institute of Chemistry, Academia Sinica, 128, Academia Road Section 2, Nangang District, Taipei 11529, Taiwan
| | - Alexander Kunitsa
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - Thomas Kus
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Arie Landau
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Keith V. Lawler
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Daniel Lefrancois
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | | | - Run R. Li
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | - Yi-Pei Li
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Jiashu Liang
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Marcus Liebenthal
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | - Hung-Hsuan Lin
- Institute of Chemistry, Academia Sinica, 128, Academia Road Section 2, Nangang District, Taipei 11529, Taiwan
| | - You-Sheng Lin
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Fenglai Liu
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | | | | | - Arne Luenser
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Aaditya Manjanath
- Institute of Chemistry, Academia Sinica, 128, Academia Road Section 2, Nangang District, Taipei 11529, Taiwan
| | - Prashant Manohar
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Erum Mansoor
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Sam F. Manzer
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Shan-Ping Mao
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | | | - Thomas Markovich
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Stephen Mason
- School of Chemistry, University of Nottingham, Nottingham, United Kingdom
| | - Simon A. Maurer
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Peter F. McLaughlin
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | | | - Jan-Michael Mewes
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Stefanie A. Mewes
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Pierpaolo Morgante
- Department of Chemistry, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | - J. Wayne Mullinax
- Department of Chemistry, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | | | | | - Alexander C. Paul
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Suranjan K. Paul
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Fabijan Pavošević
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA
| | - Zheng Pei
- School of Electrical and Computer Engineering, University of Oklahoma, Norman, Oklahoma 73019, USA
| | - Stefan Prager
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Emil I. Proynov
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Ádám Rák
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | - Eloy Ramos-Cordoba
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Bhaskar Rana
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Alan E. Rask
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Adam Rettig
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Ryan M. Richard
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Fazle Rob
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Elliot Rossomme
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Tarek Scheele
- Institute for Physical and Theoretical Chemistry, University of Bremen, Bremen, Germany
| | - Maximilian Scheurer
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Matthias Schneider
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Nickolai Sergueev
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Shaama M. Sharada
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Wojciech Skomorowski
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - David W. Small
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Christopher J. Stein
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Yu-Chuan Su
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Eric J. Sundstrom
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Zhen Tao
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA
| | - Jonathan Thirman
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Gábor J. Tornai
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | - Takashi Tsuchimochi
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Norm M. Tubman
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Oleg Vydrov
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Jan Wenzel
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Jon Witte
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Atsushi Yamada
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Kun Yao
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Sina Yeganeh
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Shane R. Yost
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Alexander Zech
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Geneva 4, Switzerland
| | - Igor Ying Zhang
- Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Xing Zhang
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Yu Zhang
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Dmitry Zuev
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Alán Aspuru-Guzik
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Alexis T. Bell
- Department of Chemical Engineering, University of California, Berkeley, California 94720, USA
| | - Nicholas A. Besley
- School of Chemistry, University of Nottingham, Nottingham, United Kingdom
| | - Ksenia B. Bravaya
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - Bernard R. Brooks
- Laboratory of Computational Biophysics, National Institute of Health, Bethesda, Maryland 20892, USA
| | - David Casanova
- Donostia International Physics Center, 20080 Donostia, Euskadi, Spain
| | | | - Sonia Coriani
- Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kgs Lyngby, Denmark
| | | | | | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | - Robert A. DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Barry D. Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Thomas R. Furlani
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - William A. Goddard
- Materials and Process Simulation Center, California Institute of Technology, Pasadena, California 91125, USA
| | | | - Teresa Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | | | | | - Yousung Jung
- Graduate School of Energy, Environment, Water and Sustainability (EEWS), Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Andreas Klamt
- COSMOlogic GmbH & Co. KG, Imbacher Weg 46, D-51379 Leverkusen, Germany
| | - Jing Kong
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Daniel S. Lambrecht
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | | | | | - C. William McCurdy
- Department of Chemistry, University of California, Davis, California 95616, USA
| | - Jeffrey B. Neaton
- Department of Physics, University of California, Berkeley, California 94720, USA
| | - Christian Ochsenfeld
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - John A. Parkhill
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Roberto Peverati
- Department of Chemistry, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | - Vitaly A. Rassolov
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, South Carolina 29208, USA
| | | | | | | | - Ryan P. Steele
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| | - Joseph E. Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Alex J. W. Thom
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Donald G. Truhlar
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Troy Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Tomasz A. Wesolowski
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Geneva 4, Switzerland
| | - K. Birgitta Whaley
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - H. Lee Woodcock
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, USA
| | - Paul M. Zimmerman
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Shirin Faraji
- Zernike Institute for Advanced Materials, University of Groningen, 9774AG Groningen, The Netherlands
| | | | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA,Author to whom correspondence should be addressed:
| |
Collapse
|
35
|
Kelpšas V, Caldararu O, Blakeley MP, Coquelle N, Wierenga RK, Ryde U, von Wachenfeldt C, Oksanen E. Neutron structures of Leishmania mexicana triosephosphate isomerase in complex with reaction-intermediate mimics shed light on the proton-shuttling steps. IUCRJ 2021; 8:633-643. [PMID: 34258011 PMCID: PMC8256706 DOI: 10.1107/s2052252521004619] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/22/2020] [Accepted: 04/30/2021] [Indexed: 06/13/2023]
Abstract
Triosephosphate isomerase (TIM) is a key enzyme in glycolysis that catalyses the interconversion of glyceraldehyde 3-phosphate and dihydroxy-acetone phosphate. This simple reaction involves the shuttling of protons mediated by protolysable side chains. The catalytic power of TIM is thought to stem from its ability to facilitate the deprotonation of a carbon next to a carbonyl group to generate an enediolate intermediate. The enediolate intermediate is believed to be mimicked by the inhibitor 2-phosphoglycolate (PGA) and the subsequent enediol intermediate by phosphoglycolohydroxamate (PGH). Here, neutron structures of Leishmania mexicana TIM have been determined with both inhibitors, and joint neutron/X-ray refinement followed by quantum refinement has been performed. The structures show that in the PGA complex the postulated general base Glu167 is protonated, while in the PGH complex it remains deprotonated. The deuteron is clearly localized on Glu167 in the PGA-TIM structure, suggesting an asymmetric hydrogen bond instead of a low-barrier hydrogen bond. The full picture of the active-site protonation states allowed an investigation of the reaction mechanism using density-functional theory calculations.
Collapse
Affiliation(s)
- Vinardas Kelpšas
- Department of Biology, Lund University, Sölvegatan 35, 223 62 Lund, Sweden
| | - Octav Caldararu
- Department of Chemistry, Lund University, 221 00 Lund, Sweden
| | - Matthew P. Blakeley
- Large-Scale Structures Group, Institut Laue–Langevin, 71 Avenue des Martyrs, 38042 Grenoble, France
| | - Nicolas Coquelle
- Large-Scale Structures Group, Institut Laue–Langevin, 71 Avenue des Martyrs, 38042 Grenoble, France
| | - Rikkert K. Wierenga
- Faculty of Biochemistry and Molecular Medicine, University of Oulu, Pentti Kaiteran katu 1, 90570 Oulu, Finland
| | - Ulf Ryde
- Department of Chemistry, Lund University, 221 00 Lund, Sweden
| | | | - Esko Oksanen
- European Spallation Source Consortium ESS ERIC, Odarslövsvägen 113, 224 84 Lund, Sweden
| |
Collapse
|
36
|
Soto-Ospina A, Araque Marín P, Bedoya G, Sepulveda-Falla D, Villegas Lanau A. Protein Predictive Modeling and Simulation of Mutations of Presenilin-1 Familial Alzheimer's Disease on the Orthosteric Site. Front Mol Biosci 2021; 8:649990. [PMID: 34150846 PMCID: PMC8206637 DOI: 10.3389/fmolb.2021.649990] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2021] [Accepted: 04/22/2021] [Indexed: 11/13/2022] Open
Abstract
Alzheimer's disease pathology is characterized by β-amyloid plaques and neurofibrillary tangles. Amyloid precursor protein is processed by β and γ secretase, resulting in the production of β-amyloid peptides with a length ranging from 38 to 43 amino acids. Presenilin 1 (PS1) is the catalytic unit of γ-secretase, and more than 200 PS1 pathogenic mutations have been identified as causative for Alzheimer's disease. A complete monocrystal structure of PS1 has not been determined so far due to the presence of two flexible domains. We have developed a complete structural model of PS1 using a computational approach with structure prediction software. Missing fragments Met1-Glut72 and Ser290-Glu375 were modeled and validated by their energetic and stereochemical characteristics. Then, with the complete structure of PS1, we defined that these fragments do not have a direct effect in the structure of the pore. Next, we used our hypothetical model for the analysis of the functional effects of PS1 mutations Ala246GLu, Leu248Pro, Leu248Arg, Leu250Val, Tyr256Ser, Ala260Val, and Val261Phe, localized in the catalytic pore. For this, we used a quantum mechanics/molecular mechanics (QM/MM) hybrid method, evaluating modifications in the topology, potential surface density, and electrostatic potential map of mutated PS1 proteins. We found that each mutation exerts changes resulting in structural modifications of the active site and in the shape of the pore. We suggest this as a valid approach for functional studies of PS1 in view of the possible impact in substrate processing and for the design of targeted therapeutic strategies.
Collapse
Affiliation(s)
- Alejandro Soto-Ospina
- Faculty of Medicine, Group Molecular Genetics, University of Antioquia, Medellín, Colombia
- Faculty of Medicine, Group Neuroscience of Antioquia, University of Antioquia, Medellín, Colombia
| | - Pedronel Araque Marín
- School of Life Sciences, Research and Innovation in Chemistry Formulations Group, EIA University, Envigado, Colombia
| | - Gabriel Bedoya
- Faculty of Medicine, Group Molecular Genetics, University of Antioquia, Medellín, Colombia
| | - Diego Sepulveda-Falla
- Faculty of Medicine, Group Neuroscience of Antioquia, University of Antioquia, Medellín, Colombia
- Molecular Neuropathology of Alzheimer’s Disease, Institute of Neuropathology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Andrés Villegas Lanau
- Faculty of Medicine, Group Molecular Genetics, University of Antioquia, Medellín, Colombia
- Faculty of Medicine, Group Neuroscience of Antioquia, University of Antioquia, Medellín, Colombia
| |
Collapse
|
37
|
Elsässer B, Goettig P. Mechanisms of Proteolytic Enzymes and Their Inhibition in QM/MM Studies. Int J Mol Sci 2021; 22:3232. [PMID: 33810118 PMCID: PMC8004986 DOI: 10.3390/ijms22063232] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2021] [Revised: 02/24/2021] [Accepted: 02/26/2021] [Indexed: 12/11/2022] Open
Abstract
Experimental evidence for enzymatic mechanisms is often scarce, and in many cases inadvertently biased by the employed methods. Thus, apparently contradictory model mechanisms can result in decade long discussions about the correct interpretation of data and the true theory behind it. However, often such opposing views turn out to be special cases of a more comprehensive and superior concept. Molecular dynamics (MD) and the more advanced molecular mechanical and quantum mechanical approach (QM/MM) provide a relatively consistent framework to treat enzymatic mechanisms, in particular, the activity of proteolytic enzymes. In line with this, computational chemistry based on experimental structures came up with studies on all major protease classes in recent years; examples of aspartic, metallo-, cysteine, serine, and threonine protease mechanisms are well founded on corresponding standards. In addition, experimental evidence from enzyme kinetics, structural research, and various other methods supports the described calculated mechanisms. One step beyond is the application of this information to the design of new and powerful inhibitors of disease-related enzymes, such as the HIV protease. In this overview, a few examples demonstrate the high potential of the QM/MM approach for sophisticated pharmaceutical compound design and supporting functions in the analysis of biomolecular structures.
Collapse
Affiliation(s)
| | - Peter Goettig
- Structural Biology Group, Department of Biosciences, University of Salzburg, Billrothstrasse 11, 5020 Salzburg, Austria;
| |
Collapse
|
38
|
Jafari S, Ryde U, Irani M. QM/MM Study of the Catalytic Reaction of Myrosinase; Importance of Assigning Proper Protonation States of Active-Site Residues. J Chem Theory Comput 2021; 17:1822-1841. [PMID: 33543623 PMCID: PMC8023669 DOI: 10.1021/acs.jctc.0c01121] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
![]()
Myrosinase from Sinapis alba hydrolyzes glycosidic
bonds of β-d-S-glucosides. The enzyme
shows an enhanced activity in the presence of l-ascorbic
acid. In this work, we employed combined quantum mechanical and molecular
mechanical (QM/MM) calculations and molecular dynamics simulations
to study the catalytic reaction of wild-type myrosinase and its E464A,
Q187A, and Q187E mutants. Test calculations show that a proper QM
region to study the myrosinase reaction must contain the whole substrate,
models of Gln-187, Glu-409, Gln-39, His-141, Asn-186, Tyr-330, Glu-464,
Arg-259, and a water molecule. Furthermore, to make the deglycosylation
step possible, Arg-259 must be charged, Glu-464 must be protonated
on OE2, and His-141 must be protonated on the NE2 atom. The results
indicate that assigning proper protonation states of the residues
is more important than the size of the model QM system. Our model
reproduces the anomeric retaining characteristic of myrosinase and
also reproduces the experimental fact that ascorbate increases the
rate of the reaction. A water molecule in the active site, positioned
by Gln-187, helps the aglycon moiety of the substrate to stabilize
the buildup of negative charge during the glycosylation reaction and
this in turn makes the moiety a better leaving group. The water molecule
also lowers the glycosylation barrier by ∼9 kcal/mol. The results
indicate that the Q187E and E464A mutants but not the Q187A mutant
can perform the glycosylation step. However, the energy profiles for
the deglycosylation step of the mutants are not similar to that of
the wild-type enzyme. The Glu-464 residue lowers the barriers of the
glycosylation and deglycosylation steps. The ascorbate ion can act
as a general base in the reaction of the wild-type enzyme only if
the Glu-464 and His-141 residues are properly protonated.
Collapse
Affiliation(s)
- Sonia Jafari
- Department of Chemistry, University of Kurdistan, 66175-416 Sanandaj, Iran
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Mehdi Irani
- Department of Chemistry, University of Kurdistan, 66175-416 Sanandaj, Iran
| |
Collapse
|
39
|
Rovaletti A, Greco C, Ryde U. QM/MM study of the binding of H 2 to MoCu CO dehydrogenase: development and applications of improved H 2 van der Waals parameters. J Mol Model 2021; 27:68. [PMID: 33538901 PMCID: PMC7862525 DOI: 10.1007/s00894-020-04655-3] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2020] [Accepted: 12/20/2020] [Indexed: 11/28/2022]
Abstract
The MoCu CO dehydrogenase enzyme not only transforms CO into CO2 but it can also oxidise H2. Even if its hydrogenase activity has been known for decades, a debate is ongoing on the most plausible mode for the binding of H2 to the enzyme active site and the hydrogen oxidation mechanism. In the present work, we provide a new perspective on the MoCu-CODH hydrogenase activity by improving the in silico description of the enzyme. Energy refinement—by means of the BigQM approach—was performed on the intermediates involved in the dihydrogen oxidation catalysis reported in our previously published work (Rovaletti, et al. “Theoretical Insights into the Aerobic Hydrogenase Activity of Molybdenum–Copper CO Dehydrogenase.” Inorganics 7 (2019) 135). A suboptimal description of the H2–HN(backbone) interaction was observed when the van der Waals parameters described in previous literature for H2 were employed. Therefore, a new set of van der Waals parameters is developed here in order to better describe the hydrogen–backbone interaction. They give rise to improved binding modes of H2 in the active site of MoCu CO dehydrogenase. Implications of the resulting outcomes for a better understanding of hydrogen oxidation catalysis mechanisms are proposed and discussed.
Collapse
Affiliation(s)
- Anna Rovaletti
- Department of Earth and Environmental Sciences, Milano-Bicocca University, Piazza della Scienza 1, 20126, Milan, Italy
| | - Claudio Greco
- Department of Earth and Environmental Sciences, Milano-Bicocca University, Piazza della Scienza 1, 20126, Milan, Italy.
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00, Lund, Sweden.
| |
Collapse
|
40
|
Abstract
Computational methods for modeling biochemical processes implemented in GAMESS package are reviewed; in particular, quantum mechanics combined with molecular mechanics (QM/MM), semi-empirical, and fragmentation approaches. A detailed summary of capabilities is provided for the QM/MM implementation in QuanPol program and the fragment molecular orbital (FMO) method. Molecular modeling and visualization packages useful for biochemical simulations with GAMESS are described. GAMESS capabilities with corresponding references are tabulated for reader's convenience.
Collapse
|
41
|
Nakliang P, Yoon S, Choi S. Emerging computational approaches for the study of regio- and stereoselectivity in organic synthesis. Org Chem Front 2021. [DOI: 10.1039/d1qo00531f] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Computational chemistry has become important in organic synthesis as it provides a detailed understanding of molecular structures and properties and detailed reaction mechanisms.
Collapse
Affiliation(s)
- Pratanphorn Nakliang
- Global AI Drug Discovery Center, College of Pharmacy and Graduate School of Pharmaceutical Sciences, Ewha Womans University, Seoul 03760, Republic of Korea
| | - Sanghee Yoon
- Global AI Drug Discovery Center, College of Pharmacy and Graduate School of Pharmaceutical Sciences, Ewha Womans University, Seoul 03760, Republic of Korea
| | - Sun Choi
- Global AI Drug Discovery Center, College of Pharmacy and Graduate School of Pharmaceutical Sciences, Ewha Womans University, Seoul 03760, Republic of Korea
- State Key Laboratory of Chemical Oncogenomics, Peking University Shenzhen Graduate School, Xili, Nanshan District, Shenzhen, 518055, China
| |
Collapse
|
42
|
Decherchi S, Cavalli A. Thermodynamics and Kinetics of Drug-Target Binding by Molecular Simulation. Chem Rev 2020; 120:12788-12833. [PMID: 33006893 PMCID: PMC8011912 DOI: 10.1021/acs.chemrev.0c00534] [Citation(s) in RCA: 95] [Impact Index Per Article: 23.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2020] [Indexed: 12/19/2022]
Abstract
Computational studies play an increasingly important role in chemistry and biophysics, mainly thanks to improvements in hardware and algorithms. In drug discovery and development, computational studies can reduce the costs and risks of bringing a new medicine to market. Computational simulations are mainly used to optimize promising new compounds by estimating their binding affinity to proteins. This is challenging due to the complexity of the simulated system. To assess the present and future value of simulation for drug discovery, we review key applications of advanced methods for sampling complex free-energy landscapes at near nonergodicity conditions and for estimating the rate coefficients of very slow processes of pharmacological interest. We outline the statistical mechanics and computational background behind this research, including methods such as steered molecular dynamics and metadynamics. We review recent applications to pharmacology and drug discovery and discuss possible guidelines for the practitioner. Recent trends in machine learning are also briefly discussed. Thanks to the rapid development of methods for characterizing and quantifying rare events, simulation's role in drug discovery is likely to expand, making it a valuable complement to experimental and clinical approaches.
Collapse
Affiliation(s)
- Sergio Decherchi
- Computational
and Chemical Biology, Fondazione Istituto
Italiano di Tecnologia, 16163 Genoa, Italy
| | - Andrea Cavalli
- Computational
and Chemical Biology, Fondazione Istituto
Italiano di Tecnologia, 16163 Genoa, Italy
- Department
of Pharmacy and Biotechnology, University
of Bologna, 40126 Bologna, Italy
| |
Collapse
|
43
|
Huix-Rotllant M, Ferré N. Analytic Energy, Gradient, and Hessian of Electrostatic Embedding QM/MM Based on Electrostatic Potential-Fitted Atomic Charges Scaling Linearly with the MM Subsystem Size. J Chem Theory Comput 2020; 17:538-548. [DOI: 10.1021/acs.jctc.0c01075] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
44
|
Cao L, Ryde U. Putative reaction mechanism of nitrogenase after dissociation of a sulfide ligand. J Catal 2020. [DOI: 10.1016/j.jcat.2020.08.028] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
45
|
Cao L, Ryde U. Quantum refinement with multiple conformations: application to the P-cluster in nitrogenase. Acta Crystallogr D Struct Biol 2020; 76:1145-1156. [PMID: 33135685 PMCID: PMC7604908 DOI: 10.1107/s2059798320012917] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2020] [Accepted: 09/21/2020] [Indexed: 11/17/2022] Open
Abstract
X-ray crystallography is the main source of atomistic information on the structure of proteins. Normal crystal structures are obtained as a compromise between the X-ray scattering data and a set of empirical restraints that ensure chemically reasonable bond lengths and angles. However, such restraints are not always available or accurate for nonstandard parts of the structure, for example substrates, inhibitors and metal sites. The method of quantum refinement, in which these empirical restraints are replaced by quantum-mechanical (QM) calculations, has previously been suggested for small but interesting parts of the protein. Here, this approach is extended to allow for multiple conformations in the QM region by performing separate QM calculations for each conformation. This approach is shown to work properly and leads to improved structures in terms of electron-density maps and real-space difference density Z-scores. It is also shown that the quality of the structures can be gauged using QM strain energies. The approach, called ComQumX-2QM, is applied to the P-cluster in two different crystal structures of the enzyme nitrogenase, i.e. an Fe8S7Cys6 cluster, used for electron transfer. One structure is at a very high resolution (1.0 Å) and shows a mixture of two different oxidation states, the fully reduced PN state (Fe82+, 20%) and the doubly oxidized P2+ state (80%). In the original crystal structure the coordinates differed for only two iron ions, but here it is shown that the two states also show differences in other atoms of up to 0.7 Å. The second structure is at a more modest resolution, 2.1 Å, and was originally suggested to show only the one-electron oxidized state, P1+. Here, it is shown that it is rather a 50/50% mixture of the P1+ and P2+ states and that many of the Fe-Fe and Fe-S distances in the original structure were quite inaccurate (by up to 0.8 Å). This shows that the new ComQumX-2QM approach can be used to sort out what is actually seen in crystal structures with dual conformations and to give locally improved coordinates.
Collapse
Affiliation(s)
- Lili Cao
- Department of Theoretical Chemistry, Lund University, PO Box 124, 221 00 Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, PO Box 124, 221 00 Lund, Sweden
| |
Collapse
|
46
|
Dong G, Lin LR, Xu LY, Li EM. Reaction mechanism of lysyl oxidase-like 2 (LOXL2) studied by computational methods. J Inorg Biochem 2020; 211:111204. [PMID: 32801097 DOI: 10.1016/j.jinorgbio.2020.111204] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2020] [Revised: 07/14/2020] [Accepted: 07/14/2020] [Indexed: 02/05/2023]
Abstract
Lysyl oxidase-like 2 (LOXL2) is a copper-dependent amine oxidase that catalyzes the oxidative deamination of the ε-amino group of lysines/hydroxylysines on substrate proteins (collagen and elastin) to form aldehyde groups. The generated aldehyde groups are of significance in crosslinking with the adjacent aldehyde or ε-amino group on proteins in extracellular matrix. In this paper, we have studied the reaction mechanism of LOXL2 by means of quantum mechanics (QM) and combined QM and molecular mechanics (QM/MM) methods. This study is divided into two parts, i.e. the biosynthesis of lysine tyrosylquinone (LTQ) cofactor and oxidative deamination of ε-amino group of lysine by LTQ. For the former part, the reaction is driven by a large exothermicity of about 284 kJ/mol. Dopaquinone radical (DPQr) is suggested to be an intermediate state in this reaction. In addition, His652 residue is predicted to serve as proton acceptor. The rate-determining step for the biosynthesis of LTQ is found to be hydrogen-atom abstraction from the benzene ring on substrate by Cu2+-hydroxide, which is a proton-coupled electron transfer (PCET) process with an energy barrier of 84 kJ/mol. For the latter part, the reaction is exothermic by about 145 kJ/mol, and the copper ion is proposed to play a role of redox catalyst in the last step to generate the product of aldehyde. However, the copper ion might not be indispensable for the latter part, which is consistent with the previous study.
Collapse
Affiliation(s)
- Geng Dong
- Department of Biochemistry and Molecular Biology, Shantou University Medical College, Shantou 515041, PR China; Medical Informatics Research Center, Shantou University Medical College, Shantou 515041, PR China.
| | - Li-Rui Lin
- Department of Biochemistry and Molecular Biology, Shantou University Medical College, Shantou 515041, PR China; Medical Informatics Research Center, Shantou University Medical College, Shantou 515041, PR China
| | - Li-Yan Xu
- Key Laboratory of Molecular Biology in High Cancer Incidence Coastal Area of Guangdong Higher Education Institutes, Shantou University Medical College, Shantou 515041, PR China; Cancer Research Center, Shantou University Medical College, Shantou 515041, PR China
| | - En-Min Li
- Department of Biochemistry and Molecular Biology, Shantou University Medical College, Shantou 515041, PR China; Key Laboratory of Molecular Biology in High Cancer Incidence Coastal Area of Guangdong Higher Education Institutes, Shantou University Medical College, Shantou 515041, PR China.
| |
Collapse
|
47
|
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: 96] [Impact Index Per Article: 24.0] [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
|
48
|
Bím D, Chalupský J, Culka M, Solomon EI, Rulíšek L, Srnec M. Proton-Electron Transfer to the Active Site Is Essential for the Reaction Mechanism of Soluble Δ 9-Desaturase. J Am Chem Soc 2020; 142:10412-10423. [PMID: 32406236 DOI: 10.1021/jacs.0c01786] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
A full understanding of the catalytic action of non-heme iron (NHFe) and non-heme diiron (NHFe2) enzymes is still beyond the grasp of contemporary computational and experimental techniques. Many of these enzymes exhibit fascinating chemo-, regio-, and stereoselectivity, in spite of employing highly reactive intermediates which are necessary for activations of most stable chemical bonds. Herein, we study in detail one intriguing representative of the NHFe2 family of enzymes: soluble Δ9 desaturase (Δ9D), which desaturates rather than performing the thermodynamically favorable hydroxylation of substrate. Its catalytic mechanism has been explored in great detail by using QM(DFT)/MM and multireference wave function methods. Starting from the spectroscopically observed 1,2-μ-peroxo diferric P intermediate, the proton-electron uptake by P is the favored mechanism for catalytic activation, since it allows a significant reduction of the barrier of the initial (and rate-determining) H-atom abstraction from the stearoyl substrate as compared to the "proton-only activated" pathway. Also, we ruled out that a Q-like intermediate (high-valent diamond-core bis-μ-oxo-[FeIV]2 unit) is involved in the reaction mechanism. Our mechanistic picture is consistent with the experimental data available for Δ9D and satisfies fairly stringent conditions required by Nature: the chemo-, stereo-, and regioselectivity of the desaturation of stearic acid. Finally, the mechanisms evaluated are placed into a broader context of NHFe2 chemistry, provided by an amino acid sequence analysis through the families of the NHFe2 enzymes. Our study thus represents an important contribution toward understanding the catalytic action of the NHFe2 enzymes and may inspire further work in NHFe(2) biomimetic chemistry.
Collapse
Affiliation(s)
- Daniel Bím
- J. Heyrovský Institute of Physical Chemistry, The Czech Academy of Sciences, Dolejškova 3, Prague 8 182 23, Czech Republic.,Institute of Organic Chemistry and Biochemistry, The Czech Academy of Sciences, Flemingovo nám. 2, Prague 6 166 10, Czech Republic
| | - Jakub Chalupský
- Institute of Organic Chemistry and Biochemistry, The Czech Academy of Sciences, Flemingovo nám. 2, Prague 6 166 10, Czech Republic
| | - Martin Culka
- Institute of Organic Chemistry and Biochemistry, The Czech Academy of Sciences, Flemingovo nám. 2, Prague 6 166 10, Czech Republic
| | - Edward I Solomon
- Department of Chemistry, Stanford University, 333 Campus Drive, Stanford, California 94305-5080, United States
| | - Lubomír Rulíšek
- Institute of Organic Chemistry and Biochemistry, The Czech Academy of Sciences, Flemingovo nám. 2, Prague 6 166 10, Czech Republic
| | - Martin Srnec
- J. Heyrovský Institute of Physical Chemistry, The Czech Academy of Sciences, Dolejškova 3, Prague 8 182 23, Czech Republic
| |
Collapse
|
49
|
Cao L, Ryde U. N 2H 2 binding to the nitrogenase FeMo cluster studied by QM/MM methods. J Biol Inorg Chem 2020; 25:521-540. [PMID: 32266560 PMCID: PMC7186253 DOI: 10.1007/s00775-020-01780-5] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2019] [Accepted: 03/25/2020] [Indexed: 12/15/2022]
Abstract
We have made a systematic combined quantum mechanical and molecular mechanical (QM/MM) investigation of possible structures of the N2 bound state of nitrogenase. We assume that N2 is immediately protonated to a N2H2 state, thereby avoiding the problem of determining the position of the protons in the cluster. We have systematically studied both end-on and side-on structures, as well as both HNNH and NNH2 states. Our results indicate that the binding of N2H2 is determined more by interactions and steric clashes with the surrounding protein than by the intrinsic preferences of the ligand and the cluster. The best binding mode with both the TPSS and B3LYP density-functional theory methods has trans-HNNH terminally bound to Fe2. It is stabilised by stacking of the substrate with His-195 and Ser-278. However, several other structures come rather close in energy (within 3-35 kJ/mol) at least in some calculations: The corresponding cis-HNNH structure terminally bound to Fe2 is second best with B3LYP. A structure with HNNH2 terminally bound to Fe6 is second most stable with TPSS (where the third proton is transferred to the substrate from the homocitrate ligand). Structures with trans-HNNH, bound to Fe4 or Fe6, or cis-HNNH bound to Fe6 are also rather stable. Finally, with the TPSS functional, a structure with cis-HNNH side-on binding to the Fe3-Fe4-Fe5-Fe7 face of the cluster is also rather low in energy, but all side-on structures are strongly disfavoured by the B3LYP method.
Collapse
Affiliation(s)
- Lili Cao
- Department of Theoretical Chemistry, Chemical Centre, Lund University, P. O. Box 124, 221 00, Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Chemical Centre, Lund University, P. O. Box 124, 221 00, Lund, Sweden.
| |
Collapse
|
50
|
Abstract
![]()
Nitrogenase
is the only enzyme that can cleave the strong triple
bond in N2. The active site contains a complicated MoFe7S9C cluster. It is believed that it needs to accept
four protons and electrons, forming the E4 state, before
it can bind N2. However, there is no consensus on the atomic
structure of the E4 state. Experimental studies indicate
that it should contain two hydride ions bridging two pairs of Fe ions,
and it has been suggested that both hydride ions as well as the two
protons bind on the same face of the cluster. On the other hand, density
functional theory (DFT) studies have indicated that it is energetically
more favorable with either three hydride ions or with a triply protonated
carbide ion, depending on the DFT functional. We have performed a
systematic combined quantum mechanical and molecular mechanical (QM/MM)
study of possible E4 states with two bridging hydride ions.
Our calculations suggest that the most favorable structure has hydride
ions bridging the Fe2/6 and Fe3/7 ion pairs. In fact, such structures
are 14 kJ/mol more stable than structures with three hydride ions,
showing that pure DFT functionals give energetically most favorable
structures in agreement with experiments. An important reason for
this finding is that we have identified a new type of broken-symmetry
state that involves only two Fe ions with minority spin, in contrast
to the previously studied states with three Fe ions with minority
spin. The energetically best structures have the two hydride ions
on different faces of the FeMo cluster, whereas better agreement with
ENDOR data is obtained if they are on the same face; such structures
are only 6–22 kJ/mol less stable.
Collapse
Affiliation(s)
- Lili Cao
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|