1
|
Wittmann L, Gordiy I, Friede M, Helmich-Paris B, Grimme S, Hansen A, Bursch M. Extension of the D3 and D4 London dispersion corrections to the full actinides series. Phys Chem Chem Phys 2024. [PMID: 39092890 DOI: 10.1039/d4cp01514b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/04/2024]
Abstract
Efficient dispersion corrections are an indispensable component of modern density functional theory, semi-empirical quantum mechanical, and even force field methods. In this work, we extend the well established D3 and D4 London dispersion corrections to the full actinides series, francium, and radium. To keep consistency with the existing versions, the original parameterization strategy of the D4 model was only slightly modified. This includes improved reference Hirshfeld atomic partial charges at the ωB97M-V/ma-def-TZVP level to fit the required electronegativity equilibration charge (EEQ) model. In this context, we developed a new actinide data set called AcQM, which covers the most common molecular actinide compound space. Furthermore, the efficient calculation of dynamic polarizabilities that are needed to construct CAB6 dispersion coefficients was implemented into the ORCA program package. The extended models are assessed for the computation of dissociation curves of actinide atoms and ions, geometry optimizations of crystal structure cutouts, gas-phase structures of small uranium compounds, and an example extracted from a small actinide complex protein assembly. We found that the novel parameterizations perform on par with the computationally more demanding density-dependent VV10 dispersion correction. With the presented extension, the excellent cost-accuracy ratio of the D3 and D4 models can now be utilized in various fields of computational actinide chemistry and, e.g., in efficient composite DFT methods such as r2SCAN-3c. They are implemented in our freely available standalone codes (dftd4, s-dftd3) and the D4 version will be also available in the upcoming ORCA 6.0 program package.
Collapse
Affiliation(s)
- Lukas Wittmann
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Igor Gordiy
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Marvin Friede
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Benjamin Helmich-Paris
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Markus Bursch
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
- FACCTs GmbH, 50677, Köln, Germany
| |
Collapse
|
2
|
de Carvalho Matias EG, Bezerra KS, Costa AHL, Clemente Junior WS, Oliveira JIN, Ribeiro Junior LA, Galvão DS, Fulco UL. Quantum biochemical analysis of the TtgR regulator and effectors. Sci Rep 2024; 14:8519. [PMID: 38609407 PMCID: PMC11015042 DOI: 10.1038/s41598-024-58441-9] [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: 01/25/2024] [Accepted: 03/29/2024] [Indexed: 04/14/2024] Open
Abstract
The recent expansion of multidrug-resistant (MDR) pathogens poses significant challenges in treating healthcare-associated infections. Although antibacterial resistance occurs by numerous mechanisms, active efflux of the drugs is a critical concern. A single species of efflux pump can produce a simultaneous resistance to several drugs. One of the best-studied efflux pumps is the TtgABC: a tripartite resistance-nodulation-division (RND) efflux pump implicated in the intrinsic antibiotic resistance in Pseudomonas putida DOT-T1E. The expression of the TtgABC gene is down-regulated by the HTH-type transcriptional repressor TtgR. In this context, by employing quantum chemistry methods based on the Density Functional Theory (DFT) within the Molecular Fragmentation with Conjugate Caps (MFCC) approach, we investigate the coupling profiles of the transcriptional regulator TtgR in complex with quercetin (QUE), a natural polyphenolic flavonoid, tetracycline (TAC), and chloramphenicol (CLM), two broad-spectrum antimicrobial agents. Our quantum biochemical computational results show the: [i] convergence radius, [ii] total binding energy, [iii] relevance (energetically) of the ligands regions, and [iv] most relevant amino acids residues of the TtgR-QUE/TAC/CLM complexes, pointing out distinctions and similarities among them. These findings improve the understanding of the binding mechanism of effectors and facilitate the development of new chemicals targeting TtgR, helping in the battle against the rise of resistance to antimicrobial drugs. These advances are crucial in the ongoing fight against rising antimicrobial drug resistance, providing hope for a future where healthcare-associated infections can be more beneficially treated.
Collapse
Affiliation(s)
- E G de Carvalho Matias
- Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, Natal, RN, 59072-970, Brazil
| | - K S Bezerra
- Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, Natal, RN, 59072-970, Brazil
| | - A H Lima Costa
- Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, Natal, RN, 59072-970, Brazil
| | - W S Clemente Junior
- Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, Natal, RN, 59072-970, Brazil
| | - J I N Oliveira
- Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, Natal, RN, 59072-970, Brazil
| | - L A Ribeiro Junior
- Institute of Physics, University of Brasília, Brasília, 70919-970, Brazil.
| | - D S Galvão
- Applied Physics Department, University of Campinas, Campinas, São Paulo, Brazil
| | - U L Fulco
- Departamento de Biofísica e Farmacologia, Universidade Federal do Rio Grande do Norte, Natal, RN, 59072-970, Brazil
| |
Collapse
|
3
|
Prasad VK, Cheng F, Fekl U, Jacobsen HA. Applications of noisy quantum computing and quantum error mitigation to "adamantaneland": a benchmarking study for quantum chemistry. Phys Chem Chem Phys 2024; 26:4071-4082. [PMID: 38225897 DOI: 10.1039/d3cp03523a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2024]
Abstract
The field of quantum computing has the potential to transform quantum chemistry. The variational quantum eigensolver (VQE) algorithm has allowed quantum computing to be applied to chemical problems in the noisy intermediate-scale quantum (NISQ) era. Applications of VQE have generally focused on predicting absolute energies instead of chemical properties that are relative energy differences and that are most interesting to chemists studying a chemical problem. We address this shortcoming by constructing a molecular benchmark data set in this work containing isomers of C10H16 and carbocationic rearrangements of C10H15+, calculated at a high-level of theory. Using the data set, we compared noiseless VQE simulations to conventionally performed density functional and wavefunction theory-based methods to understand the quality of results. We also investigated the effectiveness of a quantum state tomography-based error mitigation technique in applications of VQE under noise (simulated and real). Our findings reveal that the use of quantum error mitigation is crucial in the NISQ era and advantageous to yield almost noiseless quality results.
Collapse
Affiliation(s)
- Viki Kumar Prasad
- The Edward S. Rogers Sr. Department of Electrical and Computer Engineering, University of Toronto, 10 Kings College Road, Toronto, Ontario, Canada, M5S 3G4. arno,
- Department of Chemical and Physical Sciences, University of Toronto Mississauga, 3359 Mississauga Road, Mississauga, Ontario, Canada, L5L 1C6.
| | - Freeman Cheng
- Department of Computer Science, University of Toronto, 40 St. George Street, Toronto, Ontario, Canada, M5S 2E4
| | - Ulrich Fekl
- Department of Chemical and Physical Sciences, University of Toronto Mississauga, 3359 Mississauga Road, Mississauga, Ontario, Canada, L5L 1C6.
| | - Hans-Arno Jacobsen
- The Edward S. Rogers Sr. Department of Electrical and Computer Engineering, University of Toronto, 10 Kings College Road, Toronto, Ontario, Canada, M5S 3G4. arno,
- Department of Computer Science, University of Toronto, 40 St. George Street, Toronto, Ontario, Canada, M5S 2E4
| |
Collapse
|
4
|
Azad I, Khan T, Ahmad N, Khan AR, Akhter Y. Updates on drug designing approach through computational strategies: a review. Future Sci OA 2023; 9:FSO862. [PMID: 37180609 PMCID: PMC10167725 DOI: 10.2144/fsoa-2022-0085] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2022] [Accepted: 04/12/2023] [Indexed: 05/16/2023] Open
Abstract
The drug discovery and development (DDD) process in pursuit of novel drug candidates is a challenging procedure requiring lots of time and resources. Therefore, computer-aided drug design (CADD) methodologies are used extensively to promote proficiency in drug development in a systematic and time-effective manner. The point in reference is SARS-CoV-2 which has emerged as a global pandemic. In the absence of any confirmed drug moiety to treat the infection, the science fraternity adopted hit and trial methods to come up with a lead drug compound. This article is an overview of the virtual methodologies, which assist in finding novel hits and help in the progression of drug development in a short period with a specific medicinal solution.
Collapse
Affiliation(s)
- Iqbal Azad
- Department of Chemistry, Integral University, Dasauli, P.O. Bas-ha, Kursi Road, Lucknow, 226026, UP, India
| | - Tahmeena Khan
- Department of Chemistry, Integral University, Dasauli, P.O. Bas-ha, Kursi Road, Lucknow, 226026, UP, India
| | - Naseem Ahmad
- Department of Chemistry, Integral University, Dasauli, P.O. Bas-ha, Kursi Road, Lucknow, 226026, UP, India
| | - Abdul Rahman Khan
- Department of Chemistry, Integral University, Dasauli, P.O. Bas-ha, Kursi Road, Lucknow, 226026, UP, India
| | - Yusuf Akhter
- Department of Biotechnology, Babasaheb Bhimrao Ambedkar University, Vidya Vihar, Raebareli Road, Lucknow, UP, 2260025, India
| |
Collapse
|
5
|
Kubař T, Elstner M, Cui Q. Hybrid Quantum Mechanical/Molecular Mechanical Methods For Studying Energy Transduction in Biomolecular Machines. Annu Rev Biophys 2023; 52:525-551. [PMID: 36791746 PMCID: PMC10810093 DOI: 10.1146/annurev-biophys-111622-091140] [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] [Indexed: 02/17/2023]
Abstract
Hybrid quantum mechanical/molecular mechanical (QM/MM) methods have become indispensable tools for the study of biomolecules. In this article, we briefly review the basic methodological details of QM/MM approaches and discuss their applications to various energy transduction problems in biomolecular machines, such as long-range proton transports, fast electron transfers, and mechanochemical coupling. We highlight the particular importance for these applications of balancing computational efficiency and accuracy. Using several recent examples, we illustrate the value and limitations of QM/MM methodologies for both ground and excited states, as well as strategies for calibrating them in specific applications. We conclude with brief comments on several areas that can benefit from further efforts to make QM/MM analyses more quantitative and applicable to increasingly complex biological problems.
Collapse
Affiliation(s)
- T Kubař
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany;
| | - M Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany;
- Institute of Biological Interfaces (IBG-2), Karlsruhe Institute of Technology, Karlsruhe, Germany;
| | - Q Cui
- Department of Chemistry, Boston University, Boston, Massachusetts, USA;
- Department of Physics, Boston University, Boston, Massachusetts, USA
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts, USA
| |
Collapse
|
6
|
Rennekamp B, Karfusehr C, Kurth M, Ünal A, Monego D, Riedmiller K, Gryn'ova G, Hudson DM, Gräter F. Collagen breaks at weak sacrificial bonds taming its mechanoradicals. Nat Commun 2023; 14:2075. [PMID: 37045839 PMCID: PMC10097693 DOI: 10.1038/s41467-023-37726-z] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2022] [Accepted: 03/29/2023] [Indexed: 04/14/2023] Open
Abstract
Collagen is a force-bearing, hierarchical structural protein important to all connective tissue. In tendon collagen, high load even below macroscopic failure level creates mechanoradicals by homolytic bond scission, similar to polymers. The location and type of initial rupture sites critically decide on both the mechanical and chemical impact of these micro-ruptures on the tissue, but are yet to be explored. We here use scale-bridging simulations supported by gel electrophoresis and mass spectrometry to determine breakage points in collagen. We find collagen crosslinks, as opposed to the backbone, to harbor the weakest bonds, with one particular bond in trivalent crosslinks as the most dominant rupture site. We identify this bond as sacrificial, rupturing prior to other bonds while maintaining the material's integrity. Also, collagen's weak bonds funnel ruptures such that the potentially harmful mechanoradicals are readily stabilized. Our results suggest this unique failure mode of collagen to be tailored towards combatting an early onset of macroscopic failure and material ageing.
Collapse
Affiliation(s)
- Benedikt Rennekamp
- Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118, Heidelberg, Germany
- Interdisciplinary Center for Scientific Computing, Heidelberg University, INF 205, 69120, Heidelberg, Germany
- Max Planck School Matter to Life, Jahnstrasse 29, 69120, Heidelberg, Germany
| | - Christoph Karfusehr
- Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118, Heidelberg, Germany
- Max Planck School Matter to Life, Jahnstrasse 29, 69120, Heidelberg, Germany
- Physics Department and ZNN, Technical University Munich, Coulombwall 4a, 85748, Garching, Germany
| | - Markus Kurth
- Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118, Heidelberg, Germany
- Interdisciplinary Center for Scientific Computing, Heidelberg University, INF 205, 69120, Heidelberg, Germany
| | - Aysecan Ünal
- Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118, Heidelberg, Germany
- Interdisciplinary Center for Scientific Computing, Heidelberg University, INF 205, 69120, Heidelberg, Germany
- Max Planck School Matter to Life, Jahnstrasse 29, 69120, Heidelberg, Germany
| | - Debora Monego
- Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118, Heidelberg, Germany
| | - Kai Riedmiller
- Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118, Heidelberg, Germany
| | - Ganna Gryn'ova
- Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118, Heidelberg, Germany
- Interdisciplinary Center for Scientific Computing, Heidelberg University, INF 205, 69120, Heidelberg, Germany
| | - David M Hudson
- Department of Orthopaedics and Sports Medicine, University of Washington, Seattle, WA, 98195, USA
| | - Frauke Gräter
- Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118, Heidelberg, Germany.
- Interdisciplinary Center for Scientific Computing, Heidelberg University, INF 205, 69120, Heidelberg, Germany.
- Max Planck School Matter to Life, Jahnstrasse 29, 69120, Heidelberg, Germany.
| |
Collapse
|
7
|
Vuong VQ, Cevallos C, Hourahine B, Aradi B, Jakowski J, Irle S, Camacho C. Accelerating the density-functional tight-binding method using graphical processing units. J Chem Phys 2023; 158:084802. [PMID: 36859078 DOI: 10.1063/5.0130797] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023] Open
Abstract
Acceleration of the density-functional tight-binding (DFTB) method on single and multiple graphical processing units (GPUs) was accomplished using the MAGMA linear algebra library. Two major computational bottlenecks of DFTB ground-state calculations were addressed in our implementation: the Hamiltonian matrix diagonalization and the density matrix construction. The code was implemented and benchmarked on two different computer systems: (1) the SUMMIT IBM Power9 supercomputer at the Oak Ridge National Laboratory Leadership Computing Facility with 1-6 NVIDIA Volta V100 GPUs per computer node and (2) an in-house Intel Xeon computer with 1-2 NVIDIA Tesla P100 GPUs. The performance and parallel scalability were measured for three molecular models of 1-, 2-, and 3-dimensional chemical systems, represented by carbon nanotubes, covalent organic frameworks, and water clusters.
Collapse
Affiliation(s)
- Van-Quan Vuong
- Bredesen Center for Interdisciplinary Research and Graduate Education, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - Caterina Cevallos
- School of Chemistry, University of Costa Rica, San José 11501-2060, Costa Rica
| | - Ben Hourahine
- SUPA, Department of Physics, The John Anderson Building, 107 Rottenrow East, Glasgow G4 0NG, United Kingdom
| | - Bálint Aradi
- Bremen Center for Computational Materials Science, Universität Bremen, Bremen, Germany
| | - Jacek Jakowski
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Stephan Irle
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Cristopher Camacho
- School of Chemistry, University of Costa Rica, San José 11501-2060, Costa Rica
| |
Collapse
|
8
|
Lasisi KH, Abass OK, Zhang K, Ajibade TF, Ajibade FO, Ojediran JO, Okonofua ES, Adewumi JR, Ibikunle PD. Recent advances on graphyne and its family members as membrane materials for water purification and desalination. Front Chem 2023; 11:1125625. [PMID: 36742031 PMCID: PMC9895114 DOI: 10.3389/fchem.2023.1125625] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2022] [Accepted: 01/04/2023] [Indexed: 01/21/2023] Open
Abstract
Graphyne and its family members (GFMs) are allotropes of carbon (a class of 2D materials) having unique properties in form of structures, pores and atom hybridizations. Owing to their unique properties, GFMs have been widely utilized in various practical and theoretical applications. In the past decade, GFMs have received considerable attention in the area of water purification and desalination, especially in theoretical and computational aspects. More recently, GFMs have shown greater prospects in achieving optimal separation performance than the experimentally derived commercial polyamide membranes. In this review, recent theoretical and computational advances made in the GFMs research as it relates to water purification and desalination are summarized. Brief details on the properties of GFMs and the commonly used computational methods were described. More specifically, we systematically reviewed the various computational approaches employed with emphasis on the predicted permeability and selectivity of the GFM membranes. Finally, the current challenges limiting their large-scale practical applications coupled with the possible research directions for overcoming the challenges are proposed.
Collapse
Affiliation(s)
- Kayode Hassan Lasisi
- Key Laboratory of Urban Pollutant Conversion, Institute of Urban Environment, Chinese Academy of Sciences, Xiamen, China
| | - Olusegun K. Abass
- Department of Civil Engineering, and ReNEWACT Laboratory, Landmark University, Omu-Aran, Kwara State, Nigeria,*Correspondence: Olusegun K. Abass, ,
| | - Kaisong Zhang
- Key Laboratory of Marine Environment and Ecology, Ministry of Education, Ocean University of China, Qingdao, China
| | - Temitope Fausat Ajibade
- Department of Civil and Environmental Engineering, Federal University of Technology, Akure, Nigeria
| | | | - John O. Ojediran
- Department of Agricultural and Biosystems Engineering, Landmark University, Omu-Aran, Kwara State, Nigeria
| | | | - James Rotimi Adewumi
- Department of Civil and Environmental Engineering, Federal University of Technology, Akure, Nigeria
| | - Peter D. Ibikunle
- Department of Civil Engineering, and ReNEWACT Laboratory, Landmark University, Omu-Aran, Kwara State, Nigeria
| |
Collapse
|
9
|
Conformational analysis and water dynamics: a molecular dynamics study on the survival of a β-lactoglobulin peptide in the archaeological record. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111602] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
10
|
Pan X, Van R, Epifanovsky E, Liu J, Pu J, Nam K, Shao Y. Accelerating Ab Initio Quantum Mechanical and Molecular Mechanical (QM/MM) Molecular Dynamics Simulations with Multiple Time Step Integration and a Recalibrated Semiempirical QM/MM Hamiltonian. J Phys Chem B 2022; 126:10.1021/acs.jpcb.2c02262. [PMID: 35653199 PMCID: PMC9715852 DOI: 10.1021/acs.jpcb.2c02262] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular dynamics (MD) simulations employing ab initio quantum mechanical and molecular mechanical (ai-QM/MM) potentials are considered to be the state of the art, but the high computational cost associated with the ai-QM calculations remains a theoretical challenge for their routine application. Here, we present a modified protocol of the multiple time step (MTS) method for accelerating ai-QM/MM MD simulations of condensed-phase reactions. Within a previous MTS protocol [Nam J. Chem. Theory Comput. 2014, 10, 4175], reference forces are evaluated using a low-level (semiempirical QM/MM) Hamiltonian and employed at inner time steps to propagate the nuclear motions. Correction forces, which arise from the force differences between high-level (ai-QM/MM) and low-level Hamiltonians, are applied at outer time steps, where the MTS algorithm allows the time-reversible integration of the correction forces. To increase the outer step size, which is bound by the highest-frequency component in the correction forces, the semiempirical QM Hamiltonian is recalibrated in this work to minimize the magnitude of the correction forces. The remaining high-frequency modes, which are mainly bond stretches involving hydrogen atoms, are then removed from the correction forces. When combined with a Langevin or SIN(R) thermostat, the modified MTS-QM/MM scheme remains robust with an up to 8 (with Langevin) or 10 fs (with SIN(R)) outer time step (with 1 fs inner time steps) for the chorismate mutase system. This leads to an over 5-fold speedup over standard ai-QM/MM simulations, without sacrificing the accuracy in the predicted free energy profile of the reaction.
Collapse
Affiliation(s)
- Xiaoliang Pan
- Department of Chemistry and Biochemistry, University of Oklahoma, 101 Stephenson Parkway, Norman, Oklahoma 73019-5251, United States
| | - Richard Van
- Department of Chemistry and Biochemistry, University of Oklahoma, 101 Stephenson Parkway, Norman, Oklahoma 73019-5251, United States
| | - Evgeny Epifanovsky
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, United States
| | - Jian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Jingzhi Pu
- Department of Chemistry and Chemical Biology, Indiana University-Purdue University Indianapolis, 402 N Blackford St., LD326, Indianapolis, Indiana 46202, United States
| | - Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, Texas 76019, United States
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, 101 Stephenson Parkway, Norman, Oklahoma 73019-5251, United States
| |
Collapse
|
11
|
Prasad VK, Otero-de-la-Roza A, DiLabio GA. Small-Basis Set Density-Functional Theory Methods Corrected with Atom-Centered Potentials. J Chem Theory Comput 2022; 18:2913-2930. [PMID: 35412817 DOI: 10.1021/acs.jctc.2c00036] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Density functional theory (DFT) is currently the most popular method for modeling noncovalent interactions and thermochemistry. The accurate calculation of noncovalent interaction energies, reaction energies, and barrier heights requires choosing an appropriate functional and, typically, a relatively large basis set. Deficiencies of the density-functional approximation and the use of a limited basis set are the leading sources of error in the calculation of noncovalent and thermochemical properties in molecular systems. In this article, we present three new DFT methods based on the BLYP, M06-2X, and CAM-B3LYP functionals in combination with the 6-31G* basis set and corrected with atom-centered potentials (ACPs). ACPs are one-electron potentials that have the same form as effective-core potentials, except they do not replace any electrons. The ACPs developed in this work are used to generate energy corrections to the underlying DFT/basis-set method such that the errors in predicted chemical properties are minimized while maintaining the low computational cost of the parent methods. ACPs were developed for the elements H, B, C, N, O, F, Si, P, S, and Cl. The ACP parameters were determined using an extensive training set of 118655 data points, mostly of complete basis set coupled-cluster level quality. The target molecular properties for the ACP-corrected methods include noncovalent interaction energies, molecular conformational energies, reaction energies, barrier heights, and bond separation energies. The ACPs were tested first on the training set and then on a validation set of 42567 additional data points. We show that the ACP-corrected methods can predict the target molecular properties with accuracy close to complete basis set wavefunction theory methods, but at a computational cost of double-ζ DFT methods. This makes the new BLYP/6-31G*-ACP, M06-2X/6-31G*-ACP, and CAM-B3LYP/6-31G*-ACP methods uniquely suited to the calculation of noncovalent, thermochemical, and kinetic properties in large molecular systems.
Collapse
Affiliation(s)
- Viki Kumar Prasad
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada
| | - Alberto Otero-de-la-Roza
- Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, MALTA Consolider Team, Oviedo E-33006, Spain
| | - Gino A DiLabio
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada
| |
Collapse
|
12
|
Prasad VK, Otero-de-la-Roza A, DiLabio GA. Fast and Accurate Quantum Mechanical Modeling of Large Molecular Systems Using Small Basis Set Hartree-Fock Methods Corrected with Atom-Centered Potentials. J Chem Theory Comput 2022; 18:2208-2232. [PMID: 35313106 DOI: 10.1021/acs.jctc.1c01128] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
There has been significant interest in developing fast and accurate quantum mechanical methods for modeling large molecular systems. In this work, by utilizing a machine learning regression technique, we have developed new low-cost quantum mechanical approaches to model large molecular systems. The developed approaches rely on using one-electron Gaussian-type functions called atom-centered potentials (ACPs) to correct for the basis set incompleteness and the lack of correlation effects in the underlying minimal or small basis set Hartree-Fock (HF) methods. In particular, ACPs are proposed for ten elements common in organic and bioorganic chemistry (H, B, C, N, O, F, Si, P, S, and Cl) and four different base methods: two minimal basis sets (MINIs and MINIX) plus a double-ζ basis set (6-31G*) in combination with dispersion-corrected HF (HF-D3/MINIs, HF-D3/MINIX, HF-D3/6-31G*) and the HF-3c method. The new ACPs are trained on a very large set (73 832 data points) of noncovalent properties (interaction and conformational energies) and validated additionally on a set of 32 048 data points. All reference data are of complete basis set coupled-cluster quality, mostly CCSD(T)/CBS. The proposed ACP-corrected methods are shown to give errors in the tenths of a kcal/mol range for noncovalent interaction energies and up to 2 kcal/mol for molecular conformational energies. More importantly, the average errors are similar in the training and validation sets, confirming the robustness and applicability of these methods outside the boundaries of the training set. In addition, the performance of the new ACP-corrected methods is similar to complete basis set density functional theory (DFT) but at a cost that is orders of magnitude lower, and the proposed ACPs can be used in any computational chemistry program that supports effective-core potentials without modification. It is also shown that ACPs improve the description of covalent and noncovalent bond geometries of the underlying methods and that the improvement brought about by the application of the ACPs is directly related to the number of atoms to which they are applied, allowing the treatment of systems containing some atoms for which ACPs are not available. Overall, the ACP-corrected methods proposed in this work constitute an alternative accurate, economical, and reliable quantum mechanical approach to describe the geometries, interaction energies, and conformational energies of systems with hundreds to thousands of atoms.
Collapse
Affiliation(s)
- Viki Kumar Prasad
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
| | - Alberto Otero-de-la-Roza
- MALTA Consolider Team, Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, E-33006 Oviedo, Spain
| | - Gino A DiLabio
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
| |
Collapse
|
13
|
Deng J, Cui Q. Electronic Polarization Is Essential for the Stabilization and Dynamics of Buried Ion Pairs in Staphylococcal Nuclease Mutants. J Am Chem Soc 2022; 144:4594-4610. [PMID: 35239338 PMCID: PMC9616648 DOI: 10.1021/jacs.2c00312] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Buried charged residues play important roles in the modulation of protein stabilities and conformational dynamics and make crucial contributions to protein functions. Considering the generally nonpolar nature of protein interior, a key question concerns the contribution of electronic polarization to the stabilization and properties of buried charges. We answer this question by conducting free energy simulations using the latest polarizable CHARMM force field based on Drude oscillators for a series of Staphylococcal nuclease mutants that involve a buried Glu-Lys pair in different titration states and orientations. While a nonpolarizable model suggests that the ionized form of the buried Glu-Lys pair is more than 40 kcal/mol less stable than the charge-neutral form, the two titration states are comparable in stability when electronic polarization is included explicitly, a result better reconcilable with available experimental data. Analysis of free energy components suggests that additional stabilization of the ionized Glu-Lys pair has contributions from both the enhanced salt-bridge strength and stronger interaction between the ion-pair and surrounding protein residues and penetrated water. Despite the stronger direct interaction between Glu and Lys, the ion-pair exhibits considerably larger and faster structural fluctuations when polarization is included, due to compensation of interactions in the cavity. Collectively, observations from this work provide compelling evidence that electronic polarization is essential to the stability, hydration, dynamics, and therefore function of buried charges in proteins. Therefore, our study advocates for the explicit consideration of electronic polarization for mechanistic and engineering studies that implicate buried charged residues, such as enzymes and ion transporters.
Collapse
Affiliation(s)
- Jiahua Deng
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Qiang Cui
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States.,Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States.,Department of Biomedical Engineering, Boston University, 44 Cummington Mall, Boston, Massachusetts 02215, United States
| |
Collapse
|
14
|
Meuwly M. Atomistic Simulations for Reactions and Vibrational Spectroscopy in the Era of Machine Learning─ Quo Vadis?. J Phys Chem B 2022; 126:2155-2167. [PMID: 35286087 DOI: 10.1021/acs.jpcb.2c00212] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Atomistic simulations using accurate energy functions can provide molecular-level insight into functional motions of molecules in the gas and in the condensed phase. This Perspective delineates the present status of the field from the efforts of others and some of our own work and discusses open questions and future prospects. The combination of physics-based long-range representations using multipolar charge distributions and kernel representations for the bonded interactions is shown to provide realistic models for the exploration of the infrared spectroscopy of molecules in solution. For reactions, empirical models connecting dedicated energy functions for the reactant and product states allow statistically meaningful sampling of conformational space whereas machine-learned energy functions are superior in accuracy. The future combination of physics-based models with machine-learning techniques and integration into all-purpose molecular simulation software provides a unique opportunity to bring such dynamics simulations closer to reality.
Collapse
Affiliation(s)
- Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland
| |
Collapse
|
15
|
Prasetyo N, Hidayat Y. Lability of the first solvation shell of silver cations in liquid ammonia: A quantum mechanical charge field molecular dynamics simulation study. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.118517] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
16
|
Demapan D, Kussmann J, Ochsenfeld C, Cui Q. Factors That Determine the Variation of Equilibrium and Kinetic Properties of QM/MM Enzyme Simulations: QM Region, Conformation, and Boundary Condition. J Chem Theory Comput 2022; 18:2530-2542. [PMID: 35226489 PMCID: PMC9652774 DOI: 10.1021/acs.jctc.1c00714] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
To analyze the impact of various technical details on the results of quantum mechanical (QM)/molecular mechanical (MM) enzyme simulations, including the QM region size, catechol-O-methyltransferase (COMT) is studied as a model system using an approximate QM/MM method (DFTB3/CHARMM). The results show that key equilibrium and kinetic properties for methyl transfer in COMT exhibit limited variations with respect to the size of the QM region, which ranges from ∼100 to ∼500 atoms in this study. With extensive sampling, local and global structural characteristics of the enzyme are largely conserved across the studied QM regions, while the nature of the transition state (e.g., secondary kinetic isotope effect) and reaction exergonicity are largely maintained. Deviations in the free energy profile with different QM region sizes are similar in magnitude to those observed with changes in other simulation protocols, such as different initial enzyme conformations and boundary conditions. Electronic structural properties, such as the covariance matrix of residual charge fluctuations, appear to exhibit rather long-range correlations, especially when the peptide backbone is included in the QM region; this observation holds when a range-separated DFT approach is used as the QM region, suggesting that delocalization error is unlikely the origin. Overall, the analyses suggest that multiple simulation details determine the results of QM/MM enzyme simulations with comparable contributions.
Collapse
Affiliation(s)
- Darren Demapan
- Department of Chemistry, University of Munich (LMU), Butenandtstr. 7 (C), D-81377 Munich, Germany.,Department of Chemistry, University of Wisconsin, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Jörg Kussmann
- Department of Chemistry, University of Munich (LMU), Butenandtstr. 7 (C), D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Department of Chemistry, University of Munich (LMU), Butenandtstr. 7 (C), D-81377 Munich, Germany
| | - Qiang Cui
- Departments of Chemistry, Physics and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
17
|
Maag D, Mast T, Elstner M, Cui Q, Kubař T. O to bR transition in bacteriorhodopsin occurs through a proton hole mechanism. Proc Natl Acad Sci U S A 2021; 118:e2024803118. [PMID: 34561302 PMCID: PMC8488608 DOI: 10.1073/pnas.2024803118] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/04/2021] [Indexed: 12/27/2022] Open
Abstract
Extensive classical and quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations are used to establish the structural features of the O state in bacteriorhodopsin (bR) and its conversion back to the bR ground state. The computed free energy surface is consistent with available experimental data for the kinetics and thermodynamics of the O to bR transition. The simulation results highlight the importance of the proton release group (PRG, consisting of Glu194/204) and the conserved arginine 82 in modulating the hydration level of the protein cavity. In particular, in the O state, deprotonation of the PRG and downward rotation of Arg82 lead to elevated hydration level and a continuous water network that connects the PRG to the protonated Asp85. Proton exchange through this water network is shown by ∼0.1-μs semiempirical QM/MM free energy simulations to occur through the generation and propagation of a proton hole, which is relayed by Asp212 and stabilized by Arg82. This mechanism provides an explanation for the observation that the D85S mutant of bacteriorhodopsin pumps chloride ions. The electrostatics-hydration coupling mechanism and the involvement of all titration states of water are likely applicable to many biomolecules involved in bioenergetic transduction.
Collapse
Affiliation(s)
- Denis Maag
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
| | - Thilo Mast
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
- Institute for Biological Interfaces, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
| | - Qiang Cui
- Department of Chemistry, Boston University, Boston, MA 02215
- Department of Physics, Boston University, Boston, MA 02215
- Department of Biomedical Engineering, Boston University, Boston, MA 02215
| | - Tomáš Kubař
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany;
| |
Collapse
|
18
|
Rezaeisadat M, Salehi N, Bordbar AK. Inclusion of Levodopa into β-Cyclodextrin: A Comprehensive Computational Study. ACS OMEGA 2021; 6:23814-23825. [PMID: 34568661 PMCID: PMC8459354 DOI: 10.1021/acsomega.1c02637] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/20/2021] [Accepted: 07/05/2021] [Indexed: 06/13/2023]
Abstract
This study focused on the inclusion of levodopa (LVDP) into β-cyclodextrin (BCD) using various computational methods such as quantum mechanics (QM), molecular dynamics/steered molecular dynamics (MD/SMD), and QM/molecular mechanics/Poison-Boltzmann surface area (QM/MM/PBSA). The QM results assigned the most significant charge-transfer atoms and the higher stability of LVDP in the aqueous phase. The MD results indicate the formation of a 1:1 complex with a reasonable estimation of the effective radius of the complex, the significant contribution of hydrogen bonding in the binding energy, and the enhancement of the water solubility of LVDP. By accounting for the water hydrogen bonds and their dipolar effects, QM/MM calculations lead to the more accurate IR spectrum and binding energy of the BCD-LVDP complex. By considering carboxylic and amine functional groups' more precise arrangement, QM/MM assigns stronger hydrogen bonds between LVDP and BCD. While all the methods provide a reasonable estimation of the binding energy, the most accurate value (-4.14 kcal/mol) is obtained from QM/MM/PBSA.
Collapse
Affiliation(s)
| | - Nafiseh Salehi
- Department of Chemistry, University
of Isfahan, Isfahan 8174673441, Iran
| | | |
Collapse
|
19
|
Jing Z, Liu C, Ren P. Advanced Electrostatic Model for Monovalent Ions Based on Ab Initio Energy Decomposition. J Chem Inf Model 2021; 61:2806-2817. [PMID: 34096706 PMCID: PMC8323402 DOI: 10.1021/acs.jcim.1c00426] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Ions play important roles in the structures and functions of biomolecules. In biomolecular simulations, ions either directly interact with biomolecules or provide an ionic environment that influences electrostatic interactions of solutes. The AMOEBA+ water model has demonstrated significant advancement of the classical force field for describing molecular interactions due to its improvements on the functional forms to account for essential physics. This work expands the applicability of the AMOEBA+ model toward alkali metal (Li, Na, K, Rb, and Cs) and halogen (F, Cl, Br, and I) ions. Various quantum chemical data on ion-ion and ion-water interactions, experimental ion hydration free energies, and lattice energies of salt crystals are used in the parametrization. The final parameters are verified with other properties outside of the parametrization data, including lattice energies of additional salt crystals and ionic activity coefficients in solution. The new model captures a wide range of ion properties from the gas phase to solution phase and crystals. More importantly, AMOEBA+ provides energy components that are consistent with ab initio energy decomposition. Thus, we expect AMOEBA+ to be more general, transferable, and valuable for the interpretation of intermolecular forces in efficient classical simulations.
Collapse
Affiliation(s)
- Zhifeng Jing
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712
| |
Collapse
|
20
|
Abstract
Machine learning (ML) techniques applied to chemical reactions have a long history. The present contribution discusses applications ranging from small molecule reaction dynamics to computational platforms for reaction planning. ML-based techniques can be particularly relevant for problems involving both computation and experiments. For one, Bayesian inference is a powerful approach to develop models consistent with knowledge from experiments. Second, ML-based methods can also be used to handle problems that are formally intractable using conventional approaches, such as exhaustive characterization of state-to-state information in reactive collisions. Finally, the explicit simulation of reactive networks as they occur in combustion has become possible using machine-learned neural network potentials. This review provides an overview of the questions that can and have been addressed using machine learning techniques, and an outlook discusses challenges in this diverse and stimulating field. It is concluded that ML applied to chemistry problems as practiced and conceived today has the potential to transform the way with which the field approaches problems involving chemical reactions, in both research and academic teaching.
Collapse
Affiliation(s)
- Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland.,Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
21
|
Friederich P, Häse F, Proppe J, Aspuru-Guzik A. Machine-learned potentials for next-generation matter simulations. NATURE MATERIALS 2021; 20:750-761. [PMID: 34045696 DOI: 10.1038/s41563-020-0777-6] [Citation(s) in RCA: 112] [Impact Index Per Article: 37.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/21/2019] [Accepted: 07/17/2020] [Indexed: 05/18/2023]
Abstract
The choice of simulation methods in computational materials science is driven by a fundamental trade-off: bridging large time- and length-scales with highly accurate simulations at an affordable computational cost. Venturing the investigation of complex phenomena on large scales requires fast yet accurate computational methods. We review the emerging field of machine-learned potentials, which promises to reach the accuracy of quantum mechanical computations at a substantially reduced computational cost. This Review will summarize the basic principles of the underlying machine learning methods, the data acquisition process and active learning procedures. We highlight multiple recent applications of machine-learned potentials in various fields, ranging from organic chemistry and biomolecules to inorganic crystal structure predictions and surface science. We furthermore discuss the developments required to promote a broader use of ML potentials, and the possibility of using them to help solve open questions in materials science and facilitate fully computational materials design.
Collapse
Affiliation(s)
- Pascal Friederich
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario, Canada
- Department of Computer Science, University of Toronto, Toronto, Ontario, Canada
- Institute of Nanotechnology, Karlsruhe Institute of Technology, Eggenstein-Leopoldshafen, Germany
- Institute of Theoretical Informatics, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | - Florian Häse
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario, Canada
- Department of Computer Science, University of Toronto, Toronto, Ontario, Canada
- Vector Institute for Artificial Intelligence, Toronto, Ontario, Canada
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA, USA
| | - Jonny Proppe
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario, Canada
- Department of Computer Science, University of Toronto, Toronto, Ontario, Canada
- Institute of Physical Chemistry, Georg-August University, Göttingen, Germany
| | - Alán Aspuru-Guzik
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario, Canada.
- Department of Computer Science, University of Toronto, Toronto, Ontario, Canada.
- Vector Institute for Artificial Intelligence, Toronto, Ontario, Canada.
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA, USA.
- Lebovic Fellow, Canadian Institute for Advanced Research (CIFAR), Toronto, Ontario, Canada.
| |
Collapse
|
22
|
Yagi K, Ito S, Sugita Y. Exploring the Minimum-Energy Pathways and Free-Energy Profiles of Enzymatic Reactions with QM/MM Calculations. J Phys Chem B 2021; 125:4701-4713. [PMID: 33914537 PMCID: PMC10986901 DOI: 10.1021/acs.jpcb.1c01862] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Understanding molecular mechanisms of enzymatic reactions is of vital importance in biochemistry and biophysics. Here, we introduce new functions of hybrid quantum mechanical/molecular mechanical (QM/MM) calculations in the GENESIS program to compute the minimum-energy pathways (MEPs) and free-energy profiles of enzymatic reactions. For this purpose, an interface in GENESIS is developed to utilize a highly parallel electronic structure program, QSimulate-QM (https://qsimulate.com), calling it as a shared library from GENESIS. Second, algorithms to search the MEP are implemented, combining the string method (E et al. J. Chem. Phys. 2007, 126, 164103) with the energy minimization of the buffer MM region. The method implemented in GENESIS is applied to an enzyme, triosephosphate isomerase, which converts dihyroxyacetone phosphate to glyceraldehyde 3-phosphate in four proton-transfer processes. QM/MM-molecular dynamics simulations show performances of greater than 1 ns/day with the density functional tight binding (DFTB), and 10-30 ps/day with the hybrid density functional theory, B3LYP-D3. These performances allow us to compute not only MEP but also the potential of mean force (PMF) of the enzymatic reactions using the QM/MM calculations. The barrier height obtained as 13 kcal mol-1 with B3LYP-D3 in the QM/MM calculation is in agreement with the experimental results. The impact of conformational sampling in PMF calculations and the level of electronic structure calculations (DFTB vs B3LYP-D3) suggests reliable computational protocols for enzymatic reactions without high computational costs.
Collapse
Affiliation(s)
- Kiyoshi Yagi
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Shingo Ito
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Yuji Sugita
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, 7-1-26 minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, 1-6-5 minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
23
|
Abstract
QM/MM simulations have become an indispensable tool in many chemical and biochemical investigations. Considering the tremendous degree of success, including recognition by a 2013 Nobel Prize in Chemistry, are there still "burning challenges" in QM/MM methods, especially for biomolecular systems? In this short Perspective, we discuss several issues that we believe greatly impact the robustness and quantitative applicability of QM/MM simulations to many, if not all, biomolecules. We highlight these issues with observations and relevant advances from recent studies in our group and others in the field. Despite such limited scope, we hope the discussions are of general interest and will stimulate additional developments that help push the field forward in meaningful directions.
Collapse
Affiliation(s)
- Qiang Cui
- Departments of Chemistry, Physics, and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Tanmoy Pal
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Luke Xie
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
24
|
Kwapien K, Gavara L, Docquier J, Berthomieu D, Hernandez J, Gresh N. Intermolecular interactions of the extended recognition site of
VIM
‐2
metallo‐β‐lactamase
with 1,2,4‐triazole‐3‐thione inhibitors. Validations of a polarizable molecular mechanics potential by ab initio
QC. J Comput Chem 2020; 42:86-106. [DOI: 10.1002/jcc.26437] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2020] [Revised: 09/30/2020] [Accepted: 10/01/2020] [Indexed: 12/12/2022]
Affiliation(s)
- Karolina Kwapien
- Laboratoire de Chimie et Biochimie Pharmacologiques et Toxicologiques Université de Paris UMR 8601 Paris France
- Laboratoire de Chimie Théorique Paris France
- Institut Charles Gerhardt, UMR 5253, CNRS, Université de Montpellier, ENSCM Montpellier France
| | - Laurent Gavara
- Institut des Biomolécules Max Mousseron, UMR 5247 CNRS, Université de Montpellier, ENSCM, Faculté de Pharmacie Montpellier France
| | | | - Dorothée Berthomieu
- Institut Charles Gerhardt, UMR 5253, CNRS, Université de Montpellier, ENSCM Montpellier France
| | - Jean‐François Hernandez
- Institut des Biomolécules Max Mousseron, UMR 5247 CNRS, Université de Montpellier, ENSCM, Faculté de Pharmacie Montpellier France
| | - Nohad Gresh
- Laboratoire de Chimie Théorique Paris France
| |
Collapse
|
25
|
Lai R, Cui Q. Differences in the Nature of the Phosphoryl Transfer Transition State in Protein Phosphatase 1 and Alkaline Phosphatase: Insights from QM Cluster Models. J Phys Chem B 2020; 124:9371-9384. [PMID: 33030898 PMCID: PMC7647665 DOI: 10.1021/acs.jpcb.0c07863] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Quantum mechanical (QM) cluster models are used to probe effects on the catalytic properties of protein phosphatase 1 (PP1) and alkaline phosphatase (AP) due to metal ions and active site residues. The calculations suggest that the phosphoryl transfer transition states in PP1 are synchronous in nature with a significant degree of P-Olg cleavage, while those in AP are tighter with a modest degree of P-Olg cleavage and a range of P-Onuc formation. Similar to observations made in our recent work, a significant degree of cross talk between the forming and breaking P-O bonds complicates the interpretation of the Brønsted relation, especially in regard to AP for which the computed βlg/βEQ,lg value does not correlate with the degree of P-Olg cleavage regardless of the metal ions in the active site. By comparison, the correlation between βlg/βEQ,lg and the P-Olg bond order is more applicable to PP1, which generally exhibits less variation in the transition state than AP. Results for computational models with swapped metal ions between PP1 and AP suggest that the metal ions modulate both the nature of the transition state and the degrees of sensitivity of the transition state to the leaving group. In the reactant state, the degree of the scissile bond polarization is also different in the two enzymes, although this difference appears to be largely determined by the active site residues rather than the metal ions. Therefore, both the identity of the metal ion and the positioning of polar or charged residues in the active site contribute to the distinct catalytic characteristics of these enzymes. Several discrepancies observed between the QM cluster results and the available experimental data highlight the need for further QM/MM method developments for the quantitative analysis of metalloenzymes that contain open-shell transition metal ions.
Collapse
Affiliation(s)
- Rui Lai
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Qiang Cui
- Departments of Chemistry, Physics, and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
26
|
Lynch C, Rao S, Sansom MSP. Water in Nanopores and Biological Channels: A Molecular Simulation Perspective. Chem Rev 2020; 120:10298-10335. [PMID: 32841020 PMCID: PMC7517714 DOI: 10.1021/acs.chemrev.9b00830] [Citation(s) in RCA: 88] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2019] [Indexed: 12/18/2022]
Abstract
This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.
Collapse
Affiliation(s)
- Charlotte
I. Lynch
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| | - Shanlin Rao
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| | - Mark S. P. Sansom
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, U.K.
| |
Collapse
|
27
|
Ono J, Imai M, Nishimura Y, Nakai H. Hydroxide Ion Carrier for Proton Pumps in Bacteriorhodopsin: Primary Proton Transfer. J Phys Chem B 2020; 124:8524-8539. [DOI: 10.1021/acs.jpcb.0c05507] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Junichi Ono
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts & Batteries (ESICB), Kyoto University, 1-30 Goryo-Ohara, Nishikyo-ku, Kyoto 615-8245, Japan
| | - Minori Imai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts & Batteries (ESICB), Kyoto University, 1-30 Goryo-Ohara, Nishikyo-ku, Kyoto 615-8245, Japan
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| |
Collapse
|
28
|
Cuevas-Zuviría B, Pacios LF. Analytical Model of Electron Density and Its Machine Learning Inference. J Chem Inf Model 2020; 60:3831-3842. [DOI: 10.1021/acs.jcim.0c00197] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Bruno Cuevas-Zuviría
- Centro de Biotecnologı́a y Genómica de Plantas (CBGP, UPM-INIA), Instituto Nacional de Investigación y Tecnologı́a Agraria y Alimentaria (INIA), Universidad Politécnica de Madrid (UPM), Campus de Montegancedo-UPM, Pozuelo de Alarcón, 28223 Madrid, Spain
| | - Luis F. Pacios
- Centro de Biotecnologı́a y Genómica de Plantas (CBGP, UPM-INIA), Instituto Nacional de Investigación y Tecnologı́a Agraria y Alimentaria (INIA), Universidad Politécnica de Madrid (UPM), Campus de Montegancedo-UPM, Pozuelo de Alarcón, 28223 Madrid, Spain
- Departamento de Biotecnologı́a-Biologı́a Vegetal, Escuela Técnica Superior de Ingenierı́a Agraria, Alimentaria y de Biosistemas (ETSIAAB), Universidad Politécnica de Madrid (UPM), 28040 Madrid, Spain
| |
Collapse
|
29
|
Ito S, Cui Q. Multi-level free energy simulation with a staged transformation approach. J Chem Phys 2020; 153:044115. [PMID: 32752685 DOI: 10.1063/5.0012494] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Combining multiple levels of theory in free energy simulations to balance computational accuracy and efficiency is a promising approach for studying processes in the condensed phase. While the basic idea has been proposed and explored for quite some time, it remains challenging to achieve convergence for such multi-level free energy simulations as it requires a favorable distribution overlap between different levels of theory. Previous efforts focused on improving the distribution overlap by either altering the low-level of theory for the specific system of interest or ignoring certain degrees of freedom. Here, we propose an alternative strategy that first identifies the degrees of freedom that lead to gaps in the distributions of different levels of theory and then treats them separately with either constraints or restraints or by introducing an intermediate model that better connects the low and high levels of theory. As a result, the conversion from the low level to the high level model is done in a staged fashion that ensures a favorable distribution overlap along the way. Free energy components associated with different steps are mostly evaluated explicitly, and thus, the final result can be meaningfully compared to the rigorous free energy difference between the two levels of theory with limited and well-defined approximations. The additional free energy component calculations involve simulations at the low level of theory and therefore do not incur high computational costs. The approach is illustrated with two simple but non-trivial solution examples, and factors that dictate the reliability of the result are discussed.
Collapse
Affiliation(s)
- Shingo Ito
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| | - Qiang Cui
- Departments of Chemistry, Physics and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| |
Collapse
|
30
|
Kraus J, Gupta R, Lu M, Gronenborn AM, Akke M, Polenova T. Accurate Backbone 13 C and 15 N Chemical Shift Tensors in Galectin-3 Determined by MAS NMR and QM/MM: Details of Structure and Environment Matter. Chemphyschem 2020; 21:1436-1443. [PMID: 32363727 PMCID: PMC8080305 DOI: 10.1002/cphc.202000249] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2020] [Revised: 04/27/2020] [Indexed: 01/07/2023]
Abstract
Chemical shift tensors obtained from solid-state NMR spectroscopy are very sensitive reporters of structure and dynamics in proteins. While accurate 13 C and 15 N chemical shift tensors are accessible by magic angle spinning (MAS) NMR, their quantum mechanical calculations remain challenging, particularly for 15 N atoms. Here we compare experimentally determined backbone 13 Cα and 15 NH chemical shift tensors by MAS NMR with hybrid quantum mechanics/molecular mechanics/molecular dynamics (MD-QM/MM) calculations for the carbohydrate-binding domain of galectin-3. Excellent agreement between experimental and computed 15 NH chemical shift anisotropy values was obtained using the Amber ff15ipq force field when solvent dynamics was taken into account in the calculation. Our results establish important benchmark conditions for improving the accuracy of chemical shift calculations in proteins and may aid in the validation of protein structure models derived by MAS NMR.
Collapse
Affiliation(s)
- Jodi Kraus
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
- Pittsburgh Center for HIV Protein Interactions, University of Pittsburgh School of Medicine, 1051 Biomedical Science Tower 3, 3501 Fifth Ave., Pittsburgh, PA 15261, United States
| | - Rupal Gupta
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
- Department of Chemistry, The College of Staten Island, 2800 Victory Blvd, Staten Island, NY 10314
| | - Manman Lu
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
- Pittsburgh Center for HIV Protein Interactions, University of Pittsburgh School of Medicine, 1051 Biomedical Science Tower 3, 3501 Fifth Ave., Pittsburgh, PA 15261, United States
- Department of Structural Biology, University of Pittsburgh School of Medicine, 3501 Fifth Ave., Pittsburgh, PA 15261, United States
| | - Angela M. Gronenborn
- Pittsburgh Center for HIV Protein Interactions, University of Pittsburgh School of Medicine, 1051 Biomedical Science Tower 3, 3501 Fifth Ave., Pittsburgh, PA 15261, United States
- Department of Structural Biology, University of Pittsburgh School of Medicine, 3501 Fifth Ave., Pittsburgh, PA 15261, United States
| | - Mikael Akke
- Division of Biophysical Chemistry, Center for Molecular Protein Science, Department of Chemistry, Lund University, P.O. Box 124, SE-22100 Lund, Sweden
| | - Tatyana Polenova
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
- Pittsburgh Center for HIV Protein Interactions, University of Pittsburgh School of Medicine, 1051 Biomedical Science Tower 3, 3501 Fifth Ave., Pittsburgh, PA 15261, United States
| |
Collapse
|
31
|
Exploring the Mechanism of Catalysis with the Unified Reaction Valley Approach (URVA)—A Review. Catalysts 2020. [DOI: 10.3390/catal10060691] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
The unified reaction valley approach (URVA) differs from mainstream mechanistic studies, as it describes a chemical reaction via the reaction path and the surrounding reaction valley on the potential energy surface from the van der Waals region to the transition state and far out into the exit channel, where the products are located. The key feature of URVA is the focus on the curving of the reaction path. Moving along the reaction path, any electronic structure change of the reacting molecules is registered by a change in their normal vibrational modes and their coupling with the path, which recovers the curvature of the reaction path. This leads to a unique curvature profile for each chemical reaction with curvature minima reflecting minimal change and curvature maxima, the location of important chemical events such as bond breaking/forming, charge polarization and transfer, rehybridization, etc. A unique decomposition of the path curvature into internal coordinate components provides comprehensive insights into the origins of the chemical changes taking place. After presenting the theoretical background of URVA, we discuss its application to four diverse catalytic processes: (i) the Rh catalyzed methanol carbonylation—the Monsanto process; (ii) the Sharpless epoxidation of allylic alcohols—transition to heterogenous catalysis; (iii) Au(I) assisted [3,3]-sigmatropic rearrangement of allyl acetate; and (iv) the Bacillus subtilis chorismate mutase catalyzed Claisen rearrangement—and show how URVA leads to a new protocol for fine-tuning of existing catalysts and the design of new efficient and eco-friendly catalysts. At the end of this article the pURVA software is introduced. The overall goal of this article is to introduce to the chemical community a new protocol for fine-tuning existing catalytic reactions while aiding in the design of modern and environmentally friendly catalysts.
Collapse
|
32
|
Unke OT, Koner D, Patra S, Käser S, Meuwly M. High-dimensional potential energy surfaces for molecular simulations: from empiricism to machine learning. MACHINE LEARNING-SCIENCE AND TECHNOLOGY 2020. [DOI: 10.1088/2632-2153/ab5922] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
|
33
|
Spicher S, Abdullin D, Grimme S, Schiemann O. Modeling of spin–spin distance distributions for nitroxide labeled biomacromolecules. Phys Chem Chem Phys 2020; 22:24282-24290. [DOI: 10.1039/d0cp04920d] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Combining CREST and MD simulations based on GFN-FF for the automated computation of distance distributions for nitroxide labeled (metallo-) proteins.
Collapse
Affiliation(s)
- Sebastian Spicher
- Mulliken Center for Theoretical Chemistry
- Institute of Physical and Theoretical Chemistry
- University of Bonn
- 53115 Bonn
- Germany
| | - Dinar Abdullin
- Institute of Physical and Theoretical Chemistry
- University of Bonn
- 53115 Bonn
- Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry
- Institute of Physical and Theoretical Chemistry
- University of Bonn
- 53115 Bonn
- Germany
| | - Olav Schiemann
- Institute of Physical and Theoretical Chemistry
- University of Bonn
- 53115 Bonn
- Germany
| |
Collapse
|
34
|
Quantum mechanics/molecular mechanics multiscale modeling of biomolecules. ADVANCES IN PHYSICAL ORGANIC CHEMISTRY 2020. [DOI: 10.1016/bs.apoc.2020.08.002] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
|
35
|
Bezerra KS, Vianna JF, Lima Neto JX, Ivan Nobre Oliveira J, Albuquerque EL, Fulco UL. Interaction energies between two antiandrogenic and one androgenic agonist receptor in the presence of a T877A mutation in prostate cancer: a quantum chemistry analysis. NEW J CHEM 2020. [DOI: 10.1039/c9nj06276a] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
We provide a detailed picture of the interactions between the androgen receptor T877A mutation and three different compounds used in the treatment of prostate cancer.
Collapse
Affiliation(s)
- K. S. Bezerra
- Departamento de Biofísica e Farmacologia
- Universidade Federal do Rio Grande do Norte
- Natal-RN
- Brazil
| | - J. F. Vianna
- Departamento de Biofísica e Farmacologia
- Universidade Federal do Rio Grande do Norte
- Natal-RN
- Brazil
| | - J. X. Lima Neto
- Departamento de Biofísica e Farmacologia
- Universidade Federal do Rio Grande do Norte
- Natal-RN
- Brazil
| | - J. Ivan Nobre Oliveira
- Departamento de Biofísica e Farmacologia
- Universidade Federal do Rio Grande do Norte
- Natal-RN
- Brazil
| | - E. L. Albuquerque
- Departamento de Biofísica e Farmacologia
- Universidade Federal do Rio Grande do Norte
- Natal-RN
- Brazil
| | - U. L. Fulco
- Departamento de Biofísica e Farmacologia
- Universidade Federal do Rio Grande do Norte
- Natal-RN
- Brazil
| |
Collapse
|
36
|
Rennekamp B, Kutzki F, Obarska-Kosinska A, Zapp C, Gräter F. Hybrid Kinetic Monte Carlo/Molecular Dynamics Simulations of Bond Scissions in Proteins. J Chem Theory Comput 2019; 16:553-563. [DOI: 10.1021/acs.jctc.9b00786] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Benedikt Rennekamp
- Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
- Institute for Theoretical Physics, Heidelberg University, Philosophenweg 16, 69120 Heidelberg, Germany
| | - Fabian Kutzki
- Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| | - Agnieszka Obarska-Kosinska
- Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
- Hamburg Unit c/o DESY, European Molecular Biology Laboratory, Notkestrasse 85, 22607 Hamburg, Germany
| | - Christopher Zapp
- Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
- Institute for Theoretical Physics, Heidelberg University, Philosophenweg 16, 69120 Heidelberg, Germany
| | - Frauke Gräter
- Heidelberg Institute for Theoretical Studies, Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
- Interdisciplinary Center for Scientific Computing, Heidelberg University, INF 205, 69120 Heidelberg, Germany
| |
Collapse
|
37
|
Lu X, Duchimaza-Heredia J, Cui Q. Analysis of Density Functional Tight Binding with Natural Bonding Orbitals. J Phys Chem A 2019; 123:7439-7453. [PMID: 31373822 DOI: 10.1021/acs.jpca.9b05072] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
The description of chemical bonding by the density functional tight binding (DFTB) model is analyzed using natural bonding orbitals (NBOs) and compared to results from density functional theory (B3LYP/aug-cc-pVTZ) calculations. Several molecular systems have been chosen to represent fairly diverse bonding scenarios that include standard covalent bonds, hypervalent interactions, multicenter bonds, metal-ligand interactions (with and without the pseudo-Jahn-Teller effect), and through-space donor-acceptor interactions. Overall, the results suggest that DFTB3/3OB provides physically sound descriptions for the different bonding scenarios analyzed here, as reflected by the general agreement between DFTB3 and B3LYP NBO properties, such as the nature of the NBOs, the magnitudes of natural charges and bond orders, and the dominant donor-acceptor interactions. The degree of ligand-to-metal charge transfer and the ionic nature of pentavalent phosphate are overestimated, likely reflecting the minimal-basis nature of DFTB3/3OB. Moreover, certain orbital interactions, such as geminal interactions, are observed to be grossly overestimated by DFTB3 for hypervalent phosphate and several transition metal compounds that involve copper and nickel. The study indicates that results from NBO analysis can be instructive for identifying electronic structure descriptions at the approximate quantum-mechanical level that require improvement and thus for guiding the systematic improvement of these methods.
Collapse
Affiliation(s)
- Xiya Lu
- Department of Chemistry and Theoretical Chemistry Institute , University of Wisconsin-Madison , 1101 University Avenue , Madison , Wisconsin 53706 , United States
| | | | | |
Collapse
|
38
|
Zhou S, Wang L. Unraveling the structural and chemical features of biological short hydrogen bonds. Chem Sci 2019; 10:7734-7745. [PMID: 31588321 PMCID: PMC6764281 DOI: 10.1039/c9sc01496a] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2019] [Accepted: 06/30/2019] [Indexed: 02/06/2023] Open
Abstract
Short hydrogen bonds are ubiquitous in biological macromolecules and exhibit distinctive proton potential energy surfaces and proton sharing properties.
The three-dimensional architecture of biomolecules often creates specialized structural elements, notably short hydrogen bonds that have donor–acceptor separations below 2.7 Å. In this work, we statistically analyze 1663 high-resolution biomolecular structures from the Protein Data Bank and demonstrate that short hydrogen bonds are prevalent in proteins, protein–ligand complexes and nucleic acids. From these biological macromolecules, we characterize the preferred location, connectivity and amino acid composition in short hydrogen bonds and hydrogen bond networks, and assess their possible functional importance. Using electronic structure calculations, we further uncover how the interplay of the structural and chemical features determines the proton potential energy surfaces and proton sharing conditions in biological short hydrogen bonds.
Collapse
Affiliation(s)
- Shengmin Zhou
- Department of Chemistry and Chemical Biology , Institute for Quantitative Biomedicine , Rutgers University , Piscataway , NJ 08854 , USA .
| | - Lu Wang
- Department of Chemistry and Chemical Biology , Institute for Quantitative Biomedicine , Rutgers University , Piscataway , NJ 08854 , USA .
| |
Collapse
|
39
|
Watanabe HC, Cui Q. Quantitative Analysis of QM/MM Boundary Artifacts and Correction in Adaptive QM/MM Simulations. J Chem Theory Comput 2019; 15:3917-3928. [DOI: 10.1021/acs.jctc.9b00180] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Hiroshi C. Watanabe
- Quantum Computing Center, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama 223-8522, Japan
- Japan Science and Technology Agency, PRESTO, 4-18, Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Qiang Cui
- Departments of Chemistry, Physics, and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
40
|
Corradi V, Sejdiu BI, Mesa-Galloso H, Abdizadeh H, Noskov SY, Marrink SJ, Tieleman DP. Emerging Diversity in Lipid-Protein Interactions. Chem Rev 2019; 119:5775-5848. [PMID: 30758191 PMCID: PMC6509647 DOI: 10.1021/acs.chemrev.8b00451] [Citation(s) in RCA: 245] [Impact Index Per Article: 49.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2018] [Indexed: 02/07/2023]
Abstract
Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.
Collapse
Affiliation(s)
- Valentina Corradi
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Besian I. Sejdiu
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Haydee Mesa-Galloso
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Haleh Abdizadeh
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Sergei Yu. Noskov
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Siewert J. Marrink
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - D. Peter Tieleman
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| |
Collapse
|
41
|
Abstract
This Review illustrates the evaluation of permeability of lipid membranes from molecular dynamics (MD) simulation primarily using water and oxygen as examples. Membrane entrance, translocation, and exit of these simple permeants (one hydrophilic and one hydrophobic) can be simulated by conventional MD, and permeabilities can be evaluated directly by Fick's First Law, transition rates, and a global Bayesian analysis of the inhomogeneous solubility-diffusion model. The assorted results, many of which are applicable to simulations of nonbiological membranes, highlight the limitations of the homogeneous solubility diffusion model; support the utility of inhomogeneous solubility diffusion and compartmental models; underscore the need for comparison with experiment for both simple solvent systems (such as water/hexadecane) and well-characterized membranes; and demonstrate the need for microsecond simulations for even simple permeants like water and oxygen. Undulations, subdiffusion, fractional viscosity dependence, periodic boundary conditions, and recent developments in the field are also discussed. Last, while enhanced sampling methods and increasingly sophisticated treatments of diffusion add substantially to the repertoire of simulation-based approaches, they do not address directly the critical need for force fields with polarizability and multipoles, and constant pH methods.
Collapse
Affiliation(s)
- Richard M Venable
- Laboratory of Computational Biology, National Lung, Heart, and Blood Institute , National Institutes of Health , Bethesda , Maryland 20892 , United States
| | - Andreas Krämer
- Laboratory of Computational Biology, National Lung, Heart, and Blood Institute , National Institutes of Health , Bethesda , Maryland 20892 , United States
| | - Richard W Pastor
- Laboratory of Computational Biology, National Lung, Heart, and Blood Institute , National Institutes of Health , Bethesda , Maryland 20892 , United States
| |
Collapse
|
42
|
Thakur A, Hevel JM, Acevedo O. Examining Product Specificity in Protein Arginine Methyltransferase 7 (PRMT7) Using Quantum and Molecular Mechanical Simulations. J Chem Inf Model 2019; 59:2913-2923. [PMID: 31033288 DOI: 10.1021/acs.jcim.9b00137] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Protein arginine methyltransferase 7 (PRMT7) catalyzes the formation of monomethylarginine (MMA) but is incapable of performing a dimethylation. Given that PRMT7 performs vital functions in mammalian cells and has been implicated in a variety of diseases, including breast cancer and age-related obesity, elucidating the origin of its strict monomethylation activity is of considerable interest. Three active site residues, Glu172, Phe71, and Gln329, have been reported as particularly important for product specificity and enzymatic activity. To better understand their roles, mixed quantum and molecular mechanical (QM/MM) calculations coupled to molecular dynamics and free energy perturbation theory were carried out for the WT, F71I, and Q329S trypanosomal PRMT7 (TbPRMT7) enzymes bound with S-adenosyl- L-methionine (AdoMet) and an arginine substrate in an unmethylated or methylated form. The Q329S mutation, which experimentally abolished enzymatic activity, was appropriately computed to give an outsized Δ G‡ of 30.1 kcal/mol for MMA formation compared to 16.9 kcal/mol for WT. The F71I mutation, which has been experimentally shown to convert the enzyme from a type III PRMT into a mixed type I/II capable of forming dimethylated arginine products, yielded a reasonable Δ G‡ of 21.9 kcal/mol for the second turnover compared to 28.8 kcal/mol in the WT enzyme. Similar active site orientations for both WT and F71I TbPRMT7 allowed Glu172 and Gln329 to better orient the substrate for SN2 methylation, enhanced the nucleophilicity of the attacking guanidino group by reducing positive charge, and facilitated the binding of the subsequent methylated products.
Collapse
Affiliation(s)
- Abhishek Thakur
- Department of Chemistry , University of Miami , Coral Gables , Florida 33146 , United States
| | - Joan M Hevel
- Department of Chemistry and Biochemistry , Utah State University , Logan , Utah 84322 , United States
| | - Orlando Acevedo
- Department of Chemistry , University of Miami , Coral Gables , Florida 33146 , United States
| |
Collapse
|
43
|
Yagi K, Yamada K, Kobayashi C, Sugita Y. Anharmonic Vibrational Analysis of Biomolecules and Solvated Molecules Using Hybrid QM/MM Computations. J Chem Theory Comput 2019; 15:1924-1938. [PMID: 30730746 PMCID: PMC8864611 DOI: 10.1021/acs.jctc.8b01193] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
![]()
Quantum
mechanics/molecular mechanics (QM/MM) calculations are
applied for anharmonic vibrational analyses of biomolecules and solvated
molecules. The QM/MM method is implemented into a molecular dynamics
(MD) program, GENESIS, by interfacing with external electronic structure
programs. Following the geometry optimization and the harmonic normal-mode
analysis based on a partial Hessian, the anharmonic potential energy
surface (PES) is generated from QM/MM energies and gradients calculated
at grid points. The PES is used for vibrational self-consistent field
(VSCF) and post-VSCF calculations to compute the vibrational spectrum.
The method is first applied to a phosphate ion in solution. With both
the ion and neighboring water molecules taken as a QM region, IR spectra
of representative hydration structures are calculated by the second-order
vibrational quasi-degenerate perturbation theory (VQDPT2) at the level
of B3LYP/cc-pVTZ and TIP3P force field. A weight-average of IR spectra
over the structures reproduces the experimental spectrum with a mean
absolute deviation of 16 cm–1. Then, the method
is applied to an enzyme, P450 nitric oxide reductase (P450nor), with
the NO molecule bound to a ferric (FeIII) heme. Starting
from snapshot structures obtained from MD simulations of P450nor in
solution, QM/MM calculations have been carried out at the level of
B3LYP-D3/def2-SVP(D). The spin state of FeIII(NO) is likely
a closed-shell singlet state based on a ratio of N–O and Fe–NO
stretching frequencies (νN–O and νFe–NO) calculated for closed- and open-shell singlet
states. The calculated νN–O and νFe–NO overestimate the experimental ones by 120 and
75 cm–1, respectively. The electronic structure
and solvation of FeIII(NO) affect the structure around
the heme of P450nor leading to an increase in νN–O and νFe–NO.
Collapse
Affiliation(s)
- Kiyoshi Yagi
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Kenta Yamada
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Chigusa Kobayashi
- Computational Biophysics Research Team, RIKEN Center for Computational Science, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Computational Biophysics Research Team, RIKEN Center for Computational Science, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics Research, 1-6-5 Minatojima-Minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
44
|
Komeiji Y, Okiyama Y, Mochizuki Y, Fukuzawa K. Interaction between a Single-Stranded DNA and a Binding Protein Viewed by the Fragment Molecular Orbital Method. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2018. [DOI: 10.1246/bcsj.20180150] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Yuto Komeiji
- Biomedical Research Institute, National Institute of Industrial Science and Technology, AIST Tsukuba Central 6, Tsukuba, Ibaraki 305-8566, Japan
| | - Yoshio Okiyama
- Division of Medicinal Safety Science, National Institute of Health Sciences, 3-25-26 Tonomachi, Kawasaki-ku, Kanagawa 210-9501, Japan
| | - Yuji Mochizuki
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
- Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Kaori Fukuzawa
- Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
- School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
| |
Collapse
|
45
|
Freindorf M, Tao Y, Sethio D, Cremer D, Kraka E. New mechanistic insights into the Claisen rearrangement of chorismate – a Unified Reaction Valley Approach study. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1530464] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Marek Freindorf
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| | - Yunwen Tao
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| | - Daniel Sethio
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| | - Dieter Cremer
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| |
Collapse
|
46
|
Vujović M, Huynh M, Steiner S, Garcia-Fernandez P, Elstner M, Cui Q, Gruden M. Exploring the applicability of density functional tight binding to transition metal ions. Parameterization for nickel with the spin-polarized DFTB3 model. J Comput Chem 2018; 40:400-413. [PMID: 30299559 DOI: 10.1002/jcc.25614] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2018] [Revised: 08/18/2018] [Accepted: 09/11/2018] [Indexed: 01/04/2023]
Abstract
In this work, we explore the applicability and limitations of the current third order density functional tight binding (DFTB3) formalism for treating transition metal ions using nickel as an example. To be consistent with recent parameterization of DFTB3 for copper, the parametrization for nickel is conducted in a spin-polarized formulation and with orbital-resolved Hubbard parameters and their charge derivatives. The performance of the current parameter set is evaluated based on structural and energetic properties of a set of nickel-containing compounds that involve biologically relevant ligands. Qualitatively similar to findings in previous studies of copper complexes, the DFTB3 results are more reliable for nickel complexes with neutral ligands than for charged ligands; nevertheless, encouraging agreement is noted in comparison to the reference method, B3LYP/aug-cc-pVTZ, especially for structural properties, including cases that exhibit Jahn-Teller distortions; the structures also compare favorably to available X-ray data in the Cambridge Crystallographic Database for a number of nickel-containing compounds. As to limitations, we find it is necessary to use different d shell Hubbard charge derivatives for Ni(I) and Ni(II), due to the distinct electronic configurations for the nickel ion in the respective complexes, and substantial errors are observed for ligand binding energies, especially for charged ligands, d orbital splitting energies and splitting between singlet and triplet spin states for Ni(II) compounds. These observations highlight that future improvement in intra-d correlation and ligand polarization is required to enable the application of the DFTB3 model to complex transition metal ions. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Milena Vujović
- Center for Computational Chemistry and Bioinformatics, Faculty of Chemistry, University of Belgrade, Studentski trg 12-16 11001, Belgrade, Serbia
| | - Mioy Huynh
- Departments of Chemistry, Physics, Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts
| | - Sebastian Steiner
- Institute of Physical Chemistry & Institute of Biological Interfaces (IBG-2), Karlsruhe Institute of Technology, Kaiserstr. 12, Karlsruhe, 76131, Germany
| | - Pablo Garcia-Fernandez
- Departamento de Ciencias de la Tierra y Fısica de la Materia Condensada, Universidad de Cantabria,Cantabria Campus Internacional, Avenida de los Castros s/n 39005, Santander, Spain
| | - Marcus Elstner
- Institute of Physical Chemistry & Institute of Biological Interfaces (IBG-2), Karlsruhe Institute of Technology, Kaiserstr. 12, Karlsruhe, 76131, Germany
| | - Qiang Cui
- Departments of Chemistry, Physics, Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts
| | - Maja Gruden
- Center for Computational Chemistry and Bioinformatics, Faculty of Chemistry, University of Belgrade, Studentski trg 12-16 11001, Belgrade, Serbia
| |
Collapse
|
47
|
Hofer TS, de Visser SP. Editorial: Quantum Mechanical/Molecular Mechanical Approaches for the Investigation of Chemical Systems - Recent Developments and Advanced Applications. Front Chem 2018; 6:357. [PMID: 30271768 PMCID: PMC6146044 DOI: 10.3389/fchem.2018.00357] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2018] [Accepted: 07/30/2018] [Indexed: 12/28/2022] Open
Affiliation(s)
- Thomas S Hofer
- Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Sam P de Visser
- School of Chemical Engineering and Analytical Science, Manchester Institute of Biotechnology, The University of Manchester, Manchester, United Kingdom
| |
Collapse
|
48
|
Meuwly M. Reactive molecular dynamics: From small molecules to proteins. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1386] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Markus Meuwly
- Department of Chemistry University of Basel Basel Switzerland
- Department of Chemistry Brown University Providence Rhode Island
| |
Collapse
|
49
|
Roston D, Lu X, Fang D, Demapan D, Cui Q. Analysis of Phosphoryl-Transfer Enzymes with QM/MM Free Energy Simulations. Methods Enzymol 2018; 607:53-90. [PMID: 30149869 DOI: 10.1016/bs.mie.2018.05.005] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We discuss the application of quantum mechanics/molecular mechanics (QM/MM) free energy simulations to the analysis of phosphoryl transfers catalyzed by two enzymes: alkaline phosphatase and myosin. We focus on the nature of the transition state and the issue of mechanochemical coupling, respectively, in the two enzymes. The results illustrate unique insights that emerged from the QM/MM simulations, especially concerning the interpretation of experimental data regarding the nature of enzymatic transition states and coupling between global structural transition and catalysis in the active site. We also highlight a number of technical issues worthy of attention when applying QM/MM free energy simulations, and comment on a number of technical and mechanistic issues that require further studies.
Collapse
Affiliation(s)
- Daniel Roston
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin, Madison, Madison, WI, United States
| | - Xiya Lu
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin, Madison, Madison, WI, United States
| | - Dong Fang
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin, Madison, Madison, WI, United States
| | - Darren Demapan
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin, Madison, Madison, WI, United States
| | - Qiang Cui
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin, Madison, Madison, WI, United States.
| |
Collapse
|
50
|
Hofer TS, Hünenberger PH. Absolute proton hydration free energy, surface potential of water, and redox potential of the hydrogen electrode from first principles: QM/MM MD free-energy simulations of sodium and potassium hydration. J Chem Phys 2018; 148:222814. [DOI: 10.1063/1.5000799] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Thomas S. Hofer
- Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry, Centre for Chemistry and Biomedicine, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria
| | | |
Collapse
|