1
|
Hu X, Amin KS, Schneider M, Lim C, Salahub D, Baldauf C. System-Specific Parameter Optimization for Nonpolarizable and Polarizable Force Fields. J Chem Theory Comput 2024; 20:1448-1464. [PMID: 38279917 PMCID: PMC10867808 DOI: 10.1021/acs.jctc.3c01141] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2023] [Revised: 12/04/2023] [Accepted: 12/05/2023] [Indexed: 01/29/2024]
Abstract
The accuracy of classical force fields (FFs) has been shown to be limited for the simulation of cation-protein systems despite their importance in understanding the processes of life. Improvements can result from optimizing the parameters of classical FFs or by extending the FF formulation by terms describing charge transfer (CT) and polarization (POL) effects. In this work, we introduce our implementation of the CTPOL model in OpenMM, which extends the classical additive FF formula by adding CT and POL. Furthermore, we present an open-source parametrization tool, called FFAFFURR, that enables the (system-specific) parametrization of OPLS-AA and CTPOL models. The performance of our workflow was evaluated by its ability to reproduce quantum chemistry energies and by molecular dynamics simulations of a zinc-finger protein.
Collapse
Affiliation(s)
- Xiaojuan Hu
- Fritz-Haber-Institut
der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Kazi S. Amin
- Centre
for Molecular Simulation and Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Markus Schneider
- Fritz-Haber-Institut
der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Carmay Lim
- Institute
of Biomedical Sciences, Academia Sinica, Taipei 115, Taiwan
- Department
of Chemistry, National Tsing Hua University, Hsinchu 300, Taiwan
| | - Dennis Salahub
- Centre
for Molecular Simulation and Department of Chemistry, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Carsten Baldauf
- Fritz-Haber-Institut
der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| |
Collapse
|
2
|
Brooks CL, MacKerell AD, Post CB, Nilsson L. Biomolecular dynamics in the 21st century. Biochim Biophys Acta Gen Subj 2024; 1868:130534. [PMID: 38065235 PMCID: PMC10842176 DOI: 10.1016/j.bbagen.2023.130534] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Revised: 11/28/2023] [Accepted: 11/29/2023] [Indexed: 01/03/2024]
Abstract
The relevance of motions in biological macromolecules has been clear since the early structural analyses of proteins by X-ray crystallography. Computer simulations have been applied to provide a deeper understanding of the dynamics of biological macromolecules since 1976, and are now a standard tool in many labs working on the structure and function of biomolecules. In this mini-review we highlight some areas of current interest and active development for simulations, in particular all-atom molecular dynamics simulations.
Collapse
Affiliation(s)
- Charles L Brooks
- University of Michigan, Department of Chemistry, Ann Arbor, MI 48109, USA.
| | | | - Carol B Post
- Purdue University, Department of Medicinal Chemistry and Molecular Pharmacology, West Lafayette, IN 47907-2091, USA.
| | - Lennart Nilsson
- Karolinska Institutet, Department of Biosciences and Nutrition, SE-14183 Huddinge, Sweden.
| |
Collapse
|
3
|
Chen J, Qiu Z, Huang J. Structure and Dynamics of Confined Water Inside Diphenylalanine Peptide Nanotubes. ACS OMEGA 2023; 8:42936-42950. [PMID: 38024738 PMCID: PMC10652825 DOI: 10.1021/acsomega.3c06071] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Revised: 09/22/2023] [Accepted: 10/19/2023] [Indexed: 12/01/2023]
Abstract
Diphenylalanine (FF) peptides exhibit a unique ability to self-assemble into nanotubes with confined water molecules playing pivotal roles in their structure and function. This study investigates the structure and dynamics of diphenylalanine peptide nanotubes (FFPNTs) using all-atom molecular dynamics (MD) and grand canonical Monte Carlo combined with MD (GCMC/MD) simulations with both the CHARMM additive and Drude polarizable force fields. The occupancy and dynamics of confined water molecules were also examined. It was found that less than 2 confined water molecules per FF help stabilize the FFPNTs on the x-y plane. Analyses of the kinetics of confined water molecules revealed distinctive transport behaviors for bound and free water, and their respective diffusion coefficients were compared. Our results validate the importance of polarizable force field models in studying peptide nanotubes and provide insights into our understanding of nanoconfined water.
Collapse
Affiliation(s)
- Jinfeng Chen
- College
of Life Sciences, Zhejiang University, Hangzhou, Zhejiang 310027, China
- Key
Laboratory of Structural Biology of Zhejiang Province, School of Life
Sciences, Westlake University, Hangzhou, Zhejiang 310024, China
- Westlake
AI Therapeutics Lab, Westlake Laboratory
of Life Sciences and Biomedicine, Hangzhou, Zhejiang 310024, China
| | - Zongyang Qiu
- Key
Laboratory of Structural Biology of Zhejiang Province, School of Life
Sciences, Westlake University, Hangzhou, Zhejiang 310024, China
- Westlake
AI Therapeutics Lab, Westlake Laboratory
of Life Sciences and Biomedicine, Hangzhou, Zhejiang 310024, China
| | - Jing Huang
- Key
Laboratory of Structural Biology of Zhejiang Province, School of Life
Sciences, Westlake University, Hangzhou, Zhejiang 310024, China
- Westlake
AI Therapeutics Lab, Westlake Laboratory
of Life Sciences and Biomedicine, Hangzhou, Zhejiang 310024, China
| |
Collapse
|
4
|
Csizi K, Reiher M. Universal
QM
/
MM
approaches for general nanoscale applications. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2023. [DOI: 10.1002/wcms.1656] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
| | - Markus Reiher
- Laboratorium für Physikalische Chemie ETH Zürich Zürich Switzerland
| |
Collapse
|
5
|
Nicoli L, Giovannini T, Cappelli C. Assessing the quality of QM/MM approaches to describe vacuo-to-water solvatochromic shifts. J Chem Phys 2022; 157:214101. [PMID: 36511555 DOI: 10.1063/5.0118664] [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/13/2022] Open
Abstract
The performance of different quantum mechanics/molecular mechanics embedding models to compute vacuo-to-water solvatochromic shifts is investigated. In particular, both nonpolarizable and polarizable approaches are analyzed and computed results are compared to reference experimental data. We show that none of the approaches outperform the others and that errors strongly depend on the nature of the molecular transition to be described. Thus, we prove that the best choice of embedding model highly depends on the molecular system and that the use of a specific approach as a black box can lead to significant errors and, sometimes, totally wrong predictions.
Collapse
Affiliation(s)
- Luca Nicoli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
6
|
Cignoni E, Slama V, Cupellini L, Mennucci B. The atomistic modeling of light-harvesting complexes from the physical models to the computational protocol. J Chem Phys 2022; 156:120901. [DOI: 10.1063/5.0086275] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
The function of light-harvesting complexes is determined by a complex network of dynamic interactions among all the different components: the aggregate of pigments, the protein, and the surrounding environment. Complete and reliable predictions on these types of composite systems can be only achieved with an atomistic description. In the last few decades, there have been important advances in the atomistic modeling of light-harvesting complexes. These advances have involved both the completeness of the physical models and the accuracy and effectiveness of the computational protocols. In this Perspective, we present an overview of the main theoretical and computational breakthroughs attained so far in the field, with particular focus on the important role played by the protein and its dynamics. We then discuss the open problems in their accurate modeling that still need to be addressed. To illustrate an effective computational workflow for the modeling of light harvesting complexes, we take as an example the plant antenna complex CP29 and its H111N mutant.
Collapse
Affiliation(s)
- Edoardo Cignoni
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Vladislav Slama
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
7
|
Ambrosetti M, Skoko S, Giovannini T, Cappelli C. Quantum Mechanics/Fluctuating Charge Protocol to Compute Solvatochromic Shifts. J Chem Theory Comput 2021; 17:7146-7156. [PMID: 34619965 PMCID: PMC8582258 DOI: 10.1021/acs.jctc.1c00763] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
![]()
Despite the potentialities
of the quantum mechanics (QM)/fluctuating
charge (FQ) approach to model the spectral properties of solvated
systems, its extensive use has been hampered by the lack of reliable
parametrizations of solvents other than water. In this paper, we substantially
extend the applicability of QM/FQ to solvating environments of different
polarities and hydrogen-bonding capabilities. The reliability and
robustness of the approach are demonstrated by challenging the model
to simulate solvatochromic shifts of four organic chromophores, which
display large shifts when dissolved in apolar, aprotic or polar, protic
solvents.
Collapse
Affiliation(s)
| | - Sulejman Skoko
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
8
|
Schlick T, Portillo-Ledesma S, Myers CG, Beljak L, Chen J, Dakhel S, Darling D, Ghosh S, Hall J, Jan M, Liang E, Saju S, Vohr M, Wu C, Xu Y, Xue E. Biomolecular Modeling and Simulation: A Prospering Multidisciplinary Field. Annu Rev Biophys 2021; 50:267-301. [PMID: 33606945 PMCID: PMC8105287 DOI: 10.1146/annurev-biophys-091720-102019] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We reassess progress in the field of biomolecular modeling and simulation, following up on our perspective published in 2011. By reviewing metrics for the field's productivity and providing examples of success, we underscore the productive phase of the field, whose short-term expectations were overestimated and long-term effects underestimated. Such successes include prediction of structures and mechanisms; generation of new insights into biomolecular activity; and thriving collaborations between modeling and experimentation, including experiments driven by modeling. We also discuss the impact of field exercises and web games on the field's progress. Overall, we note tremendous success by the biomolecular modeling community in utilization of computer power; improvement in force fields; and development and application of new algorithms, notably machine learning and artificial intelligence. The combined advances are enhancing the accuracy andscope of modeling and simulation, establishing an exemplary discipline where experiment and theory or simulations are full partners.
Collapse
Affiliation(s)
- Tamar Schlick
- Department of Chemistry, New York University, New York, New York 10003, USA;
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA
- New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 200122, China
| | | | - Christopher G Myers
- Department of Chemistry, New York University, New York, New York 10003, USA;
| | - Lauren Beljak
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Justin Chen
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Sami Dakhel
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Daniel Darling
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Sayak Ghosh
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Joseph Hall
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Mikaeel Jan
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Emily Liang
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Sera Saju
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Mackenzie Vohr
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Chris Wu
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Yifan Xu
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Eva Xue
- College of Arts and Science, New York University, New York, New York 10003, USA
| |
Collapse
|
9
|
Giovannini T, Koch H. Energy-Based Molecular Orbital Localization in a Specific Spatial Region. J Chem Theory Comput 2021; 17:139-150. [PMID: 33337150 DOI: 10.1021/acs.jctc.0c00737] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
We present a novel energy-based localization procedure able to localize molecular orbitals into predefined spatial regions. The method is defined in a multiscale framework based on the multilevel Hartree-Fock approach. In particular, the system is partitioned into active and inactive fragments. The localized molecular orbitals are obtained maximizing the repulsion between the two fragments. The method is applied to several cases including both conjugated and non-conjugated systems. Our multiscale approach is compared with reference values for both ground-state properties, such as dipole moments, and local excitation energies. The proposed approach is useful to extend the application range of high-level electron correlation methods. In fact, the reduced number of molecular orbitals can lead to a large reduction in the computational cost of correlated calculations.
Collapse
Affiliation(s)
- Tommaso Giovannini
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - Henrik Koch
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
10
|
Gheorghiu A, Coveney PV, Arabi AA. The influence of base pair tautomerism on single point mutations in aqueous DNA. Interface Focus 2020; 10:20190120. [PMID: 33178413 PMCID: PMC7653342 DOI: 10.1098/rsfs.2019.0120] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/24/2020] [Indexed: 12/11/2022] Open
Abstract
The relationship between base pair hydrogen bond proton transfer and the rate of spontaneous single point mutations at ambient temperatures and pressures in aqueous DNA is investigated. By using an ensemble-based multiscale computational modelling method, statistically robust rates of proton transfer for the A:T and G:C base pairs within a solvated DNA dodecamer are calculated. Several different proton transfer pathways are observed within the same base pair. It is shown that, in G:C, the double proton transfer tautomer is preferred, while the single proton transfer process is favoured in A:T. The reported range of rate coefficients for double proton transfer is consistent with recent experimental data. Notwithstanding the approximately 1000 times more common presence of single proton transfer products from A:T, observationally there is bias towards G:C to A:T mutations in a wide range of living organisms. We infer that the double proton transfer reactions between G:C base pairs have a negligible contribution towards this bias for the following reasons: (i) the maximum half-life of the G*:C* tautomer is in the range of picoseconds, which is significantly smaller than the milliseconds it takes for DNA to unwind during replication, (ii) statistically, the majority of G*:C* tautomers revert back to their canonical forms through a barrierless process, and (iii) the thermodynamic instability of the tautomers with respect to the canonical base pairs. Through similar reasoning, we also deduce that proton transfer in the A:T base pair does not contribute to single point mutations in DNA.
Collapse
Affiliation(s)
- A Gheorghiu
- Centre for Computational Science, University College London, London, UK
| | - P V Coveney
- Centre for Computational Science, University College London, London, UK.,Informatics Institute, University of Amsterdam, Amsterdam, The Netherlands
| | - A A Arabi
- Centre for Computational Science, University College London, London, UK.,College of Medicine and Health Sciences, Biochemistry Department, United Arab Emirates University, PO Box 17666, Al Ain, United Arab Emirates
| |
Collapse
|
11
|
Conti I, Cerullo G, Nenov A, Garavelli M. Ultrafast Spectroscopy of Photoactive Molecular Systems from First Principles: Where We Stand Today and Where We Are Going. J Am Chem Soc 2020; 142:16117-16139. [PMID: 32841559 PMCID: PMC7901644 DOI: 10.1021/jacs.0c04952] [Citation(s) in RCA: 41] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
![]()
Computational spectroscopy is becoming a mandatory tool for the interpretation of the
complex, and often congested, spectral maps delivered by modern non-linear multi-pulse
techniques. The fields of Electronic Structure Methods,
Non-Adiabatic Molecular Dynamics, and Theoretical
Spectroscopy represent the three pillars of the virtual ultrafast
optical spectrometer, able to deliver transient spectra in
silico from first principles. A successful simulation strategy requires a
synergistic approach that balances between the three fields, each one having its very
own challenges and bottlenecks. The aim of this Perspective is to demonstrate that,
despite these challenges, an impressive agreement between theory and experiment is
achievable now regarding the modeling of ultrafast photoinduced processes in complex
molecular architectures. Beyond that, some key recent developments in the three fields
are presented that we believe will have major impacts on spectroscopic simulations in
the very near future. Potential directions of development, pending challenges, and
rising opportunities are illustrated.
Collapse
Affiliation(s)
- Irene Conti
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
| | - Giulio Cerullo
- Dipartimento di Fisica, Politecnico di Milano, IFN-CNR, Piazza Leonardo da Vinci 32, I-20133 Milano, Italy
| | - Artur Nenov
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
| | - Marco Garavelli
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
| |
Collapse
|
12
|
Marefat Khah A, Reinholdt P, Nuernberger P, Kongsted J, Hättig C. Relaxation Dynamics of the Triazene Compound Berenil in DNA-Minor-Groove Confinement after Photoexcitation. J Chem Theory Comput 2020; 16:5203-5211. [DOI: 10.1021/acs.jctc.0c00489] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, 5230 Odense, Denmark
| | - Patrick Nuernberger
- Institut für Physikalische und Theoretische Chemie, Universität Regensburg, 93040 Regensburg, Germany
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, 5230 Odense, Denmark
| | - Christof Hättig
- Quantum Chemistry Group, Ruhr University of Bochum, D-44780 Bochum, Germany
| |
Collapse
|
13
|
Marefat Khah A, Reinholdt P, Olsen JMH, Kongsted J, Hättig C. Avoiding Electron Spill-Out in QM/MM Calculations on Excited States with Simple Pseudopotentials. J Chem Theory Comput 2020; 16:1373-1381. [PMID: 32031806 DOI: 10.1021/acs.jctc.9b01162] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
QM/MM calculations of electronic excitations with diffuse basis sets often have large errors due to spill-out of electrons from the quantum subsystem. The Pauli repulsion of the electrons by the environment has to be included to avoid this. We propose transferable atomic all-electron pseudopotentials that can readily be combined with most MM force fields to avoid electron spill-out. QM/MM excitation energies computed with time-dependent Hartree-Fock and the algebraic diagrammatic construction through second-order are benchmarked against supermolecular calculations to validate these new pseudopotentials. The QM/MM calculations with pseudopotentials give accurate results that are stable with augmentation of the basis set with diffuse functions. We show that the largest contribution to residual deviations from full QM calculations is caused by the missing London dispersion interaction.
Collapse
Affiliation(s)
| | - Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense, Denmark
| | - Jógvan Magnus Haugaard Olsen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense, Denmark
| | - Christof Hättig
- Quantum Chemistry Group, Ruhr University of Bochum, D-44801 Bochum, Germany
| |
Collapse
|
14
|
Inakollu VS, Geerke DP, Rowley CN, Yu H. Polarisable force fields: what do they add in biomolecular simulations? Curr Opin Struct Biol 2020; 61:182-190. [PMID: 32044671 DOI: 10.1016/j.sbi.2019.12.012] [Citation(s) in RCA: 51] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2019] [Revised: 12/18/2019] [Accepted: 12/19/2019] [Indexed: 12/11/2022]
Abstract
The quality of biomolecular simulations critically depends on the accuracy of the force field used to calculate the potential energy of the molecular configurations. Currently, most simulations employ non-polarisable force fields, which describe electrostatic interactions as the sum of Coulombic interactions between fixed atomic charges. Polarisation of these charge distributions is incorporated only in a mean-field manner. In the past decade, extensive efforts have been devoted to developing simple, efficient, and yet generally applicable polarisable force fields for biomolecular simulations. In this review, we summarise the latest developments in accounting for key biomolecular interactions with polarisable force fields and applications to address challenging biological questions. In the end, we provide an outlook for future development in polarisable force fields.
Collapse
Affiliation(s)
- Vs Sandeep Inakollu
- School of Chemistry and Molecular Bioscience, University of Wollongong, Wollongong NSW 2522, Australia; Molecular Horizons, University of Wollongong, Wollongong NSW 2522 Australia; Illawarra Health and Medical Research Institute, Wollongong NSW 2522, Australia
| | - Daan P Geerke
- AIMMS Division of Molecular and Computational Toxicology, Department of Chemistry and Pharmaceutical Sciences, Vrije Universiteit Amsterdam, De Boelelaan 1108, 1081 HZ Amsterdam, the Netherlands.
| | - Christopher N Rowley
- Department of Chemistry, Memorial University of Newfoundland, St. John's, Newfoundland and Labrador, Canada.
| | - Haibo Yu
- School of Chemistry and Molecular Bioscience, University of Wollongong, Wollongong NSW 2522, Australia; Molecular Horizons, University of Wollongong, Wollongong NSW 2522 Australia; Illawarra Health and Medical Research Institute, Wollongong NSW 2522, Australia.
| |
Collapse
|
15
|
Giovannini T, Puglisi A, Ambrosetti M, Cappelli C. Polarizable QM/MM Approach with Fluctuating Charges and Fluctuating Dipoles: The QM/FQFμ Model. J Chem Theory Comput 2019; 15:2233-2245. [DOI: 10.1021/acs.jctc.8b01149] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
| | | | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
16
|
Dziedzic J, Head-Gordon T, Head-Gordon M, Skylaris CK. Mutually polarizable QM/MM model with in situ optimized localized basis functions. J Chem Phys 2019; 150:074103. [PMID: 30795653 DOI: 10.1063/1.5080384] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
We extend our recently developed quantum-mechanical/molecular mechanics (QM/MM) approach [Dziedzic et al., J. Chem. Phys. 145, 124106 (2016)] to enable in situ optimization of the localized orbitals. The quantum subsystem is described with onetep linear-scaling density functional theory and the classical subsystem - with the AMOEBA polarizable force field. The two subsystems interact via multipolar electrostatics and are fully mutually polarizable. A total energy minimization scheme is employed for the Hamiltonian of the coupled QM/MM system. We demonstrate that, compared to simpler models using fixed basis sets, the additional flexibility offered by in situ optimized basis functions improves the accuracy of the QM/MM interface, but also poses new challenges, making the QM subsystem more prone to overpolarization and unphysical charge transfer due to increased charge penetration. We show how these issues can be efficiently solved by replacing the classical repulsive van der Waals term for QM/MM interactions with an interaction of the electronic density with a fixed, repulsive MM potential that mimics Pauli repulsion, together with a modest increase in the damping of QM/MM polarization. We validate our method, with particular attention paid to the hydrogen bond, in tests on water-ion pairs, the water dimer, first solvation shells of neutral and charged species, and solute-solvent interaction energies. As a proof of principle, we determine suitable repulsive potential parameters for water, K+, and Cl-. The mechanisms we employed to counteract the unphysical overpolarization of the QM subsystem are demonstrated to be adequate, and our approach is robust. We find that the inclusion of explicit polarization in the MM part of QM/MM improves agreement with fully QM calculations. Our model permits the use of minimal size QM regions and, remarkably, yields good energetics across the well-balanced QM/MM interface.
Collapse
Affiliation(s)
- Jacek Dziedzic
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Teresa Head-Gordon
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Pitzer Theory Center and Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
17
|
Xie L, Cheng H, Fang D, Chen ZN, Yang M. Enhanced QM/MM sampling for free energy calculation of chemical reactions: A case study of double proton transfer. J Chem Phys 2019; 150:044111. [PMID: 30709281 DOI: 10.1063/1.5072779] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Free energy calculations for chemical reactions with a steep energy barrier require well defined reaction coordinates (RCs). However, when multiple parallel channels exist along selected RC, the application of conventional enhanced samplings is difficult to generate correct sampling within limited simulation time and thus cannot give correct prediction about the favorable pathways, the relative stability of multiple products or intermediates. Here, we implement the selective integrated tempering sampling (SITS) method with quantum mechanical and molecular mechanical (QM/MM) potential to investigate the chemical reactions in solution. The combined SITS-QM/MM scheme is used to identify possible reaction paths, intermediate and product states, and the free energy profiles for the different reaction paths. Two double proton transfer reactions were studied to validate the implemented method and simulation protocol, from which the independent and correlated proton transfer processes are identified in two representative systems, respectively. This protocol can be generalized to various kinds of chemical reactions for both academic studies and industry applications, such as in exploration and optimization of potential reactions in DNA encoded compound library and halogen or deuterium substitution of the hit discovery and lead optimization stages of drug design via providing a better understanding of the reaction mechanism along the designed chemical reaction pathways.
Collapse
Affiliation(s)
- Liangxu Xie
- Institute of Bioinformatics and Medical Engineering, School of Electrical and Information Engineering, Jiangsu University of Technology, Changzhou 21300, China
| | - Huimin Cheng
- XtalPi Inc. (Shenzhen Jingtai Technology Co., Ltd.), Times Science & Tech Mansion E. 20F, 7028 Shennan Ave, Futian District, Shenzhen, China
| | - Dong Fang
- XtalPi Inc. (Shenzhen Jingtai Technology Co., Ltd.), Times Science & Tech Mansion E. 20F, 7028 Shennan Ave, Futian District, Shenzhen, China
| | - Zhe-Ning Chen
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, 155 Yangqiao Road West, Fuzhou 350002, China
| | - Mingjun Yang
- XtalPi Inc. (Shenzhen Jingtai Technology Co., Ltd.), Times Science & Tech Mansion E. 20F, 7028 Shennan Ave, Futian District, Shenzhen, China
| |
Collapse
|
18
|
Sahoo SK, Nair NN. Interfacing the Core-Shell or the Drude Polarizable Force Field With Car-Parrinello Molecular Dynamics for QM/MM Simulations. Front Chem 2018; 6:275. [PMID: 30042939 PMCID: PMC6048201 DOI: 10.3389/fchem.2018.00275] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2018] [Accepted: 06/18/2018] [Indexed: 11/13/2022] Open
Abstract
We report a quantum mechanics/polarizable–molecular mechanics (QM/p–MM) potential based molecular dynamics (MD) technique where the core–shell (or the Drude) type polarizable MM force field is interfaced with the plane-wave density functional theory based QM force field which allows Car–Parrinello MD for the QM subsystem. In the QM/p-MM Lagrangian proposed here, the shell (or the Drude) MM variables are treated as extended degrees of freedom along with the Kohn–Sham (KS) orbitals describing the QM wavefunction. The shell and the KS orbital degrees of freedom are then adiabatically decoupled from the nuclear degrees of freedom. In this respect, we also present here the Nosé–Hoover Chain thermostat implementation for the dynamical subsystems. Our approach is then used to investigate the effect of MM polarization on the QM/MM results. Especially, the consequence of MM polarization on reaction free energy barriers, defect formation energy, and structural and dynamical properties are investigated. A low point charge polarizable potential (p–MZHB) for pure siliceous systems is also reported here.
Collapse
Affiliation(s)
- Sudhir K Sahoo
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, India
| | - Nisanth N Nair
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, India
| |
Collapse
|
19
|
Huang J, Simmonett AC, Pickard FC, MacKerell AD, Brooks BR. Mapping the Drude polarizable force field onto a multipole and induced dipole model. J Chem Phys 2018; 147:161702. [PMID: 29096511 DOI: 10.1063/1.4984113] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
The induced dipole and the classical Drude oscillator represent two major approaches for the explicit inclusion of electronic polarizability into force field-based molecular modeling and simulations. In this work, we explore the equivalency of these two models by comparing condensed phase properties computed using the Drude force field and a multipole and induced dipole (MPID) model. Presented is an approach to map the electrostatic model optimized in the context of the Drude force field onto the MPID model. Condensed phase simulations on water and 15 small model compounds show that without any reparametrization, the MPID model yields properties similar to the Drude force field with both models yielding satisfactory reproduction of a range of experimental values and quantum mechanical data. Our results illustrate that the Drude oscillator model and the point induced dipole model are different representations of essentially the same physical model. However, results indicate the presence of small differences between the use of atomic multipoles and off-center charge sites. Additionally, results on the use of dispersion particle mesh Ewald further support its utility for treating long-range Lennard Jones dispersion contributions in the context of polarizable force fields. The main motivation in demonstrating the transferability of parameters between the Drude and MPID models is that the more than 15 years of development of the Drude polarizable force field can now be used with MPID formalism without the need for dual-thermostat integrators nor self-consistent iterations. This opens up a wide range of new methodological opportunities for polarizable models.
Collapse
Affiliation(s)
- Jing Huang
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, 20 Penn St., Baltimore, Maryland 21201, USA
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, 5635 Fishers Lane, Rockville, Maryland 20852, USA
| | - Frank C Pickard
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, 5635 Fishers Lane, Rockville, Maryland 20852, USA
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, 20 Penn St., Baltimore, Maryland 21201, USA
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, 5635 Fishers Lane, Rockville, Maryland 20852, USA
| |
Collapse
|
20
|
Affiliation(s)
- Juan J. Nogueira
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, A-1090 Wien, Austria;,
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, A-1090 Wien, Austria;,
| |
Collapse
|
21
|
Das S, Nam K, Major DT. Rapid Convergence of Energy and Free Energy Profiles with Quantum Mechanical Size in Quantum Mechanical–Molecular Mechanical Simulations of Proton Transfer in DNA. J Chem Theory Comput 2018; 14:1695-1705. [DOI: 10.1021/acs.jctc.7b00964] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Susanta Das
- Department of Chemistry, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - Kwangho Nam
- Department of Chemistry, Umeå University, 901 87 Umeå, Sweden
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, Texas 76019-0065, United States
| | - Dan Thomas Major
- Department of Chemistry, Bar-Ilan University, Ramat-Gan 5290002, Israel
| |
Collapse
|
22
|
Boulanger E, Harvey JN. QM/MM methods for free energies and photochemistry. Curr Opin Struct Biol 2018; 49:72-76. [PMID: 29414514 DOI: 10.1016/j.sbi.2018.01.003] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2017] [Revised: 12/17/2017] [Accepted: 01/02/2018] [Indexed: 11/27/2022]
Abstract
Hybrid computational methods describing a small region of a biomolecular system with quantum mechanics and the bulk with molecular mechanics, referred to as QM/MM methods, are now a central part of computational biochemistry. This review considers developments in the QM/MM approach that make it easier to calculate free energies using accurate QM-based potential energy expressions. We also describe techniques to treat electronic coupling between the core region and the MM environment. Polarizability of the protein matrix is important but so is electronic coupling. Applications of these new methods, especially to photochemistry, are discussed.
Collapse
Affiliation(s)
- Eliot Boulanger
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium.
| | - Jeremy N Harvey
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium.
| |
Collapse
|
23
|
Cerqueira NMFSA, Fernandes PA, Ramos MJ. Protocol for Computational Enzymatic Reactivity Based on Geometry Optimisation. Chemphyschem 2018; 19:669-689. [DOI: 10.1002/cphc.201700339] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2017] [Revised: 09/22/2017] [Indexed: 01/12/2023]
Affiliation(s)
- N. M. F. S. A. Cerqueira
- REQUIMTE-UCIBIO; Departamento de Química e Bioquímica; Faculdade de Ciências; Universidade do Porto; Rua do Campo Alegre s/n 4169-007 Porto Portugal
| | - P. A. Fernandes
- REQUIMTE-UCIBIO; Departamento de Química e Bioquímica; Faculdade de Ciências; Universidade do Porto; Rua do Campo Alegre s/n 4169-007 Porto Portugal
| | - M. J. Ramos
- REQUIMTE-UCIBIO; Departamento de Química e Bioquímica; Faculdade de Ciências; Universidade do Porto; Rua do Campo Alegre s/n 4169-007 Porto Portugal
| |
Collapse
|
24
|
Feldt J, Miranda S, Pratas F, Roma N, Tomás P, Mata RA. Optimization and benchmarking of a perturbative Metropolis Monte Carlo quantum mechanics/molecular mechanics program. J Chem Phys 2017; 147:244105. [PMID: 29289129 DOI: 10.1063/1.5009820] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
In this work, we present an optimized perturbative quantum mechanics/molecular mechanics (QM/MM) method for use in Metropolis Monte Carlo simulations. The model adopted is particularly tailored for the simulation of molecular systems in solution but can be readily extended to other applications, such as catalysis in enzymatic environments. The electrostatic coupling between the QM and MM systems is simplified by applying perturbation theory to estimate the energy changes caused by a movement in the MM system. This approximation, together with the effective use of GPU acceleration, leads to a negligible added computational cost for the sampling of the environment. Benchmark calculations are carried out to evaluate the impact of the approximations applied and the overall computational performance.
Collapse
Affiliation(s)
- Jonas Feldt
- Institut für Physikalische Chemie, Universität Göttingen, Tammannstrasse 6, D-37077 Göttingen, Germany
| | - Sebastião Miranda
- INESC-ID/IST, Universidade de Lisboa, Rua Alves Redol, 1000-029 Lisboa, Portugal
| | - Frederico Pratas
- INESC-ID/IST, Universidade de Lisboa, Rua Alves Redol, 1000-029 Lisboa, Portugal
| | - Nuno Roma
- INESC-ID/IST, Universidade de Lisboa, Rua Alves Redol, 1000-029 Lisboa, Portugal
| | - Pedro Tomás
- INESC-ID/IST, Universidade de Lisboa, Rua Alves Redol, 1000-029 Lisboa, Portugal
| | - Ricardo A Mata
- Institut für Physikalische Chemie, Universität Göttingen, Tammannstrasse 6, D-37077 Göttingen, Germany
| |
Collapse
|
25
|
Gatto P, Lipparini F, Stamm B. Computation of forces arising from the polarizable continuum model within the domain-decomposition paradigm. J Chem Phys 2017; 147:224108. [DOI: 10.1063/1.5008329] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Paolo Gatto
- Mathematics Division, Center for Computational Engineering Science, RWTH Aachen University, Aachen, Germany
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Benjamin Stamm
- Mathematics Division, Center for Computational Engineering Science, RWTH Aachen University, Aachen, Germany
| |
Collapse
|
26
|
Ganguly A, Boulanger E, Thiel W. Importance of MM Polarization in QM/MM Studies of Enzymatic Reactions: Assessment of the QM/MM Drude Oscillator Model. J Chem Theory Comput 2017; 13:2954-2961. [DOI: 10.1021/acs.jctc.7b00016] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Abir Ganguly
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Eliot Boulanger
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
27
|
Mao Y, Shao Y, Dziedzic J, Skylaris CK, Head-Gordon T, Head-Gordon M. Performance of the AMOEBA Water Model in the Vicinity of QM Solutes: A Diagnosis Using Energy Decomposition Analysis. J Chem Theory Comput 2017; 13:1963-1979. [PMID: 28430427 DOI: 10.1021/acs.jctc.7b00089] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The importance of incorporating solvent polarization effects into the modeling of solvation processes has been well-recognized, and therefore a new generation of hybrid quantum mechanics/molecular mechanics (QM/MM) approaches that accounts for this effect is desirable. We present a fully self-consistent, mutually polarizable QM/MM scheme using the AMOEBA force field, in which the total energy of the system is variationally minimized with respect to both the QM electronic density and the MM induced dipoles. This QM/AMOEBA model is implemented through the Q-Chem/LibEFP code interface and then applied to the evaluation of solute-solvent interaction energies for various systems ranging from the water dimer to neutral and ionic solutes (NH3, NH4+, CN-) surrounded by increasing numbers of water molecules (up to 100). In order to analyze the resulting interaction energies, we also utilize an energy decomposition analysis (EDA) scheme which identifies contributions from permanent electrostatics, polarization, and van der Waals (vdW) interaction for the interaction between the QM solute and the solvent molecules described by AMOEBA. This facilitates a component-wise comparison against full QM calculations where the corresponding energy components are obtained via a modified version of the absolutely localized molecular orbitals (ALMO)-EDA. The results show that the present QM/AMOEBA model can yield reasonable solute-solvent interaction energies for neutral and cationic species, while further scrutiny reveals that this accuracy highly relies on the delicate balance between insufficiently favorable permanent electrostatics and softened vdW interaction. For anionic solutes where the charge penetration effect becomes more pronounced, the QM/MM interface turns out to be unbalanced. These results are consistent with and further elucidate our findings in a previous study using a slightly different QM/AMOEBA model ( Dziedzic et al. J. Chem. Phys. 2016 , 145 , 124106 ). The implications of these results for further refinement of this model are also discussed.
Collapse
Affiliation(s)
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma , Norman, Oklahoma 73019, United States
| | - Jacek Dziedzic
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, U.K.,Faculty of Applied Physics and Mathematics, Gdańsk University of Technology , Gdańsk 80-233, Poland
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, U.K
| | | | - Martin Head-Gordon
- Chemical Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
28
|
Cole DJ, Hine NDM. Applications of large-scale density functional theory in biology. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:393001. [PMID: 27494095 DOI: 10.1088/0953-8984/28/39/393001] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Density functional theory (DFT) has become a routine tool for the computation of electronic structure in the physics, materials and chemistry fields. Yet the application of traditional DFT to problems in the biological sciences is hindered, to a large extent, by the unfavourable scaling of the computational effort with system size. Here, we review some of the major software and functionality advances that enable insightful electronic structure calculations to be performed on systems comprising many thousands of atoms. We describe some of the early applications of large-scale DFT to the computation of the electronic properties and structure of biomolecules, as well as to paradigmatic problems in enzymology, metalloproteins, photosynthesis and computer-aided drug design. With this review, we hope to demonstrate that first principles modelling of biological structure-function relationships are approaching a reality.
Collapse
Affiliation(s)
- Daniel J Cole
- Theory of Condensed Matter group, Cavendish Laboratory, 19 JJ Thomson Ave, Cambridge CB3 0HE, UK. School of Chemistry, Newcastle University, Newcastle upon Tyne NE1 7RU, UK
| | | |
Collapse
|
29
|
Dziedzic J, Mao Y, Shao Y, Ponder J, Head-Gordon T, Head-Gordon M, Skylaris CK. TINKTEP: A fully self-consistent, mutually polarizable QM/MM approach based on the AMOEBA force field. J Chem Phys 2016; 145:124106. [DOI: 10.1063/1.4962909] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Jacek Dziedzic
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
- Faculty of Applied Physics and Mathematics, Gdańsk University of Technology, Gdańsk, Poland
| | - Yuezhi Mao
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Yihan Shao
- Q-Chem Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Jay Ponder
- Department of Chemistry, Washington University in St. Louis, St. Louis, Missouri 63130, USA
| | - Teresa Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, USA
- Department of Bioengineering, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
30
|
Bennie SJ, van der Kamp MW, Pennifold RCR, Stella M, Manby FR, Mulholland AJ. A Projector-Embedding Approach for Multiscale Coupled-Cluster Calculations Applied to Citrate Synthase. J Chem Theory Comput 2016; 12:2689-97. [PMID: 27159381 DOI: 10.1021/acs.jctc.6b00285] [Citation(s) in RCA: 50] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Projector-based embedding has recently emerged as a robust multiscale method for the calculation of various electronic molecular properties. We present the coupling of projector embedding with quantum mechanics/molecular mechanics modeling and apply it for the first time to an enzyme-catalyzed reaction. Using projector-based embedding, we combine coupled-cluster theory, density-functional theory (DFT), and molecular mechanics to compute energies for the proton abstraction from acetyl-coenzyme A by citrate synthase. By embedding correlated ab initio methods in DFT we eliminate functional sensitivity and obtain high-accuracy profiles in a procedure that is straightforward to apply.
Collapse
Affiliation(s)
- Simon J Bennie
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| | - Marc W van der Kamp
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
- School of Biochemistry, Biomedical Sciences Building, University Walk, University of Bristol , Bristol BS8 1TD, U.K
| | - Robert C R Pennifold
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| | - Martina Stella
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| | - Frederick R Manby
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| | - Adrian J Mulholland
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| |
Collapse
|
31
|
Kratz EG, Walker AR, Lagardère L, Lipparini F, Piquemal JP, Cisneros GA. LICHEM: A QM/MM program for simulations with multipolar and polarizable force fields. J Comput Chem 2016; 37:1019-29. [PMID: 26781073 PMCID: PMC4808410 DOI: 10.1002/jcc.24295] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2015] [Revised: 09/21/2015] [Accepted: 12/17/2015] [Indexed: 12/23/2022]
Abstract
We introduce an initial implementation of the LICHEM software package. LICHEM can interface with Gaussian, PSI4, NWChem, TINKER, and TINKER-HP to enable QM/MM calculations using multipolar/polarizable force fields. LICHEM extracts forces and energies from unmodified QM and MM software packages to perform geometry optimizations, single-point energy calculations, or Monte Carlo simulations. When the QM and MM regions are connected by covalent bonds, the pseudo-bond approach is employed to smoothly transition between the QM region and the polarizable force field. A series of water clusters and small peptides have been employed to test our initial implementation. The results obtained from these test systems show the capabilities of the new software and highlight the importance of including explicit polarization. © 2016 Wiley Periodicals, Inc.
Collapse
|
32
|
Dixit M, Das S, Mhashal A, Eitan R, Major D. Practical Aspects of Multiscale Classical and Quantum Simulations of Enzyme Reactions. Methods Enzymol 2016; 577:251-86. [DOI: 10.1016/bs.mie.2016.05.046] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
|
33
|
Mei Y, Simmonett AC, Pickard FC, DiStasio RA, Brooks BR, Shao Y. Numerical study on the partitioning of the molecular polarizability into fluctuating charge and induced atomic dipole contributions. J Phys Chem A 2015; 119:5865-82. [PMID: 25945749 DOI: 10.1021/acs.jpca.5b03159] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In order to carry out a detailed analysis of the molecular static polarizability, which is the response of the molecule to a uniform external electric field, the molecular polarizability was computed using the finite-difference method for 21 small molecules, using density functional theory. Within nine charge population schemes (Löwdin, Mulliken, Becke, Hirshfeld, CM5, Hirshfeld-I, NPA, CHELPG, MK-ESP) in common use, the charge fluctuation contribution is found to dominate the molecular polarizability, with its ratio ranging from 59.9% with the Hirshfeld or CM5 scheme to 96.2% with the Mulliken scheme. The Hirshfeld-I scheme is also used to compute the other contribution to the molecular polarizability coming from the induced atomic dipoles, and the atomic polarizabilities in eight small molecules and water pentamer are found to be highly anisotropic for most atoms. Overall, the results suggest that (a) more emphasis probably should be placed on the charge fluctuation terms in future polarizable force field development and (b) an anisotropic polarizability might be more suitable than an isotropic one in polarizable force fields based entirely or partially on the induced atomic dipoles.
Collapse
Affiliation(s)
- Ye Mei
- †State Key Laboratory of Precision Spectroscopy, Department of Physics and Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China.,‡NYU-ECNU Center for Computational Chemistry, NYU Shanghai, Shanghai 200062, China.,⊥Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Andrew C Simmonett
- ⊥Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Frank C Pickard
- ⊥Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Robert A DiStasio
- §Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Bernard R Brooks
- ⊥Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Yihan Shao
- ∥Q-Chem Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, United States
| |
Collapse
|
34
|
Chung LW, Sameera WMC, Ramozzi R, Page AJ, Hatanaka M, Petrova GP, Harris TV, Li X, Ke Z, Liu F, Li HB, Ding L, Morokuma K. The ONIOM Method and Its Applications. Chem Rev 2015; 115:5678-796. [PMID: 25853797 DOI: 10.1021/cr5004419] [Citation(s) in RCA: 743] [Impact Index Per Article: 82.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Lung Wa Chung
- †Department of Chemistry, South University of Science and Technology of China, Shenzhen 518055, China
| | - W M C Sameera
- ‡Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo, Kyoto 606-8103, Japan
| | - Romain Ramozzi
- ‡Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo, Kyoto 606-8103, Japan
| | - Alister J Page
- §Newcastle Institute for Energy and Resources, The University of Newcastle, Callaghan 2308, Australia
| | - Miho Hatanaka
- ‡Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo, Kyoto 606-8103, Japan
| | - Galina P Petrova
- ∥Faculty of Chemistry and Pharmacy, University of Sofia, Bulgaria Boulevard James Bourchier 1, 1164 Sofia, Bulgaria
| | - Travis V Harris
- ‡Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo, Kyoto 606-8103, Japan.,⊥Department of Chemistry, State University of New York at Oswego, Oswego, New York 13126, United States
| | - Xin Li
- #State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Zhuofeng Ke
- ∇School of Chemistry and Chemical Engineering, Sun Yat-sen University, Guangzhou 510275, China
| | - Fengyi Liu
- ○Key Laboratory of Macromolecular Science of Shaanxi Province, School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an, Shaanxi 710119, China
| | - Hai-Bei Li
- ■School of Ocean, Shandong University, Weihai 264209, China
| | - Lina Ding
- ▲School of Pharmaceutical Sciences, Zhengzhou University, 100 Kexue Avenue, Zhengzhou, Henan 450001, China
| | - Keiji Morokuma
- ‡Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo, Kyoto 606-8103, Japan
| |
Collapse
|
35
|
Olsen JMH, Steinmann C, Ruud K, Kongsted J. Polarizable Density Embedding: A New QM/QM/MM-Based Computational Strategy. J Phys Chem A 2015; 119:5344-55. [DOI: 10.1021/jp510138k] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jógvan Magnus Haugaard Olsen
- Laboratory of
Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
- Department
of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense
M, Denmark
| | - Casper Steinmann
- Department
of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense
M, Denmark
| | - Kenneth Ruud
- Centre for Theoretical
and Computational Chemistry, Department of Chemistry, University of Tromsø—The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Jacob Kongsted
- Department
of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense
M, Denmark
| |
Collapse
|
36
|
Caprasecca S, Jurinovich S, Lagardère L, Stamm B, Lipparini F. Achieving Linear Scaling in Computational Cost for a Fully Polarizable MM/Continuum Embedding. J Chem Theory Comput 2015; 11:694-704. [DOI: 10.1021/ct501087m] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Stefano Caprasecca
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe
Moruzzi 3, I-56124 Pisa, Italy
| | - Sandro Jurinovich
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe
Moruzzi 3, I-56124 Pisa, Italy
| | - Louis Lagardère
- Sorbonne Universités, UPMC Univ. Paris 06, Institut du Calcul et de la Simulation, F-75005 Paris, France
| | - Benjamin Stamm
- Sorbonne Universités, UPMC Univ. Paris 06, UMR 7598, Laboratoire Jacques-Louis Lions, F-75005 Paris, France
- CNRS, UMR
7598
and 7616, F-75005 Paris, France
| | - Filippo Lipparini
- Sorbonne Universités, UPMC Univ. Paris 06, Institut du Calcul et de la Simulation, F-75005 Paris, France
- Sorbonne Universités, UPMC Univ. Paris 06, UMR 7598, Laboratoire Jacques-Louis Lions, F-75005 Paris, France
- Sorbonne Universités, UPMC Univ. Paris 06, UMR 7616, Laboratoire de Chimie Théorique, F-75005 Paris, France
| |
Collapse
|
37
|
Riahi S, Rowley CN. The CHARMM-TURBOMOLE interface for efficient and accurate QM/MM molecular dynamics, free energies, and excited state properties. J Comput Chem 2014; 35:2076-86. [PMID: 25178266 DOI: 10.1002/jcc.23716] [Citation(s) in RCA: 45] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2014] [Revised: 08/04/2014] [Accepted: 08/07/2014] [Indexed: 01/05/2023]
Abstract
The quantum mechanical (QM)/molecular mechanical (MM) interface between Chemistry at HARvard Molecular Mechanics (CHARMM) and TURBOMOLE is described. CHARMM provides an extensive set of simulation algorithms, like molecular dynamics (MD) and free energy perturbation, and support for mature nonpolarizable and Drude polarizable force fields. TURBOMOLE provides fast QM calculations using density functional theory or wave function methods and excited state properties. CHARMM-TURBOMOLE is well-suited for extended QM/MM MD simulations using first principles methods with large (triple-ζ) basis sets. We demonstrate these capabilities with a QM/MM simulation of Mg(2+) (aq), where the MM outer sphere water molecules are represented using the SWM4-NDP Drude polarizable force field and the ion and inner coordination sphere are represented using QM PBE, PBE0, and MP2 methods. The relative solvation free energies of Mg(2+) and Zn(2+) were calculated using thermodynamic integration. We also demonstrate the features for excited state properties. We calculate the time-averaged solution absorption spectrum of indole, the emission spectrum of the indole 1La excited state, and the electronic circular dichroism spectrum of an oxacepham.
Collapse
Affiliation(s)
- Saleh Riahi
- Department of Chemistry, Memorial University of Newfoundland, St. John's, Newfoundland, A1B 3X7, Canada
| | | |
Collapse
|
38
|
Blake HL, Robinson D. QM/MM studies of contemporary and novel membrane raft fluorescent probes. Molecules 2014; 19:10230-41. [PMID: 25029071 PMCID: PMC6271554 DOI: 10.3390/molecules190710230] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2014] [Revised: 07/09/2014] [Accepted: 07/10/2014] [Indexed: 11/22/2022] Open
Abstract
We have studied a number of contemporary and novel membrane probes, selected for their structural similarity to membrane raft components, in order to properly anchor themselves within a sphingolipid/cholesterol rich region. A QM/MM approach was adopted in order to understand the structural and electrostatic influences of fluorescence emission shifts of the probes in different lipid and solvation environments. The proposed modifications to the membrane probes have shown encouraging data relating not only to emission shifts within the membrane, but also their ability to anchor within a membrane raft domain and the stability to internalization within a membrane system.
Collapse
Affiliation(s)
- Hannah L Blake
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, UK.
| | - David Robinson
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, UK.
| |
Collapse
|