1
|
Song C, Wang LP. A Polarizable QM/MM Model That Combines the State-Averaged CASSCF and AMOEBA Force Field for Photoreactions in Proteins. J Chem Theory Comput 2024. [PMID: 39088696 DOI: 10.1021/acs.jctc.4c00623] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/03/2024]
Abstract
This study presents the polarizable quantum mechanics/molecular mechanics (QM/MM) embedding of the state-averaged complete active space self-consistent field (SA-CASSCF) in the atomic multipole optimized energetics for biomolecular applications (AMOEBA) force field for the purpose of studying photoreactions in protein environments. We describe two extensions of our previous work that combine SA-CASSCF with AMOEBA water models, allowing it to be generalized to AMOEBA models for proteins and other macromolecules. First, we discuss how our QM/MM model accounts for the discrepancy between the direct and polarization electric fields that arises in the AMOEBA description of intramolecular polarization. A second improvement is the incorporation of link atom schemes to treat instances in which the QM/MM boundary goes through covalent bonds. A single-link atom scheme and double-link atom scheme are considered in this work, and we will discuss how electrostatic interaction, van der Waals interaction, and various kinds of valence terms are treated across the boundary. To test the accuracy of the link atom scheme, we will compare QM/MM with full QM calculations and study how the errors in ground state properties, excited state properties, and excitation energies change when tuning the parameters in the link atom scheme. We will also test the new SA-CASSCF/AMOEBA method on an elementary reaction step in NanoLuc, an artificial bioluminescence luciferase. We will show how the reaction mechanism is different when calculated in the gas phase, in polarizable continuum medium (PCM), versus in protein AMOEBA models.
Collapse
Affiliation(s)
- Chenchen Song
- Department of Chemistry, University of California, Davis, 1 Shields Avenue, Davis, California 95616, United States
| | - Lee-Ping Wang
- Department of Chemistry, University of California, Davis, 1 Shields Avenue, Davis, California 95616, United States
| |
Collapse
|
2
|
Polonius S, Lehrner D, González L, Mai S. Resolving Photoinduced Femtosecond Three-Dimensional Solute-Solvent Dynamics through Surface Hopping Simulations. J Chem Theory Comput 2024; 20:4738-4750. [PMID: 38768386 PMCID: PMC11171268 DOI: 10.1021/acs.jctc.4c00169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Revised: 04/22/2024] [Accepted: 05/07/2024] [Indexed: 05/22/2024]
Abstract
Photoinduced dynamics in solution is governed by mutual solute-solvent interactions, which give rise to phenomena like solvatochromism, the Stokes shift, dual fluorescence, or charge transfer. Understanding these phenomena requires simulating the solute's photoinduced dynamics and simultaneously resolving the three-dimensional solvent distribution dynamics. If using trajectory surface hopping (TSH) to this aim, thousands of trajectories are required to adequately sample the time-dependent three-dimensional solvent distribution functions, and thus resolve the solvent dynamics with sub-Ångstrom and femtosecond accuracy and sufficiently low noise levels. Unfortunately, simulating thousands of trajectories with TSH in the framework of hybrid quantum mechanical/molecular mechanical (QM/MM) can be prohibitively expensive when employing ab initio electronic structure methods. To tackle this challenge, we recently introduced a computationally efficient approach that combines efficient linear vibronic coupling models with molecular mechanics (LVC/MM) via electrostatic embedding [Polonius et al., JCTC 2023, 19, 7171-7186]. This method provides solvent-embedded, nonadiabatically coupled potential energy surfaces while scaling similarly to MM force fields. Here, we employ TSH with LVC/MM to unravel the photoinduced dynamics of two small thiocarbonyl compounds solvated in water. We describe how to estimate the number of trajectories required to produce nearly noise-free three-dimensional solvent distribution functions and present an analysis based on approximately 10,000 trajectories propagated for 3 ps. In the electronic ground state, both molecules exhibit in-plane hydrogen bonds to the sulfur atom. Shortly after excitation, these bonds are broken and reform perpendicular to the molecular plane on timescales that differ by an order of magnitude due to steric effects. We also show that the solvent relaxation dynamics is coupled to the electronic dynamics, including intersystem crossing. These findings are relevant to advance the understanding of the coupled solute-solvent dynamics of solvated photoexcited molecules, e.g., biologically relevant thio-nucleobases.
Collapse
Affiliation(s)
- Severin Polonius
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry (DoSChem), University of Vienna, Währinger Str. 42, 1090 Vienna, Austria
| | - David Lehrner
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| |
Collapse
|
3
|
Bondanza M, Nottoli T, Nottoli M, Cupellini L, Lipparini F, Mennucci B. The OpenMMPol library for polarizable QM/MM calculations of properties and dynamics. J Chem Phys 2024; 160:134106. [PMID: 38557842 DOI: 10.1063/5.0198251] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Accepted: 03/13/2024] [Indexed: 04/04/2024] Open
Abstract
We present a new library designed to provide a simple and straightforward way to implement QM/AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications) and other polarizable QM/MM (Molecular Mechanics) methods based on induced point dipoles. The library, herein referred to as OpenMMPol, is free and open-sourced and is engineered to address the increasing demand for accurate and efficient QM/MM simulations. OpenMMPol is specifically designed to allow polarizable QM/MM calculations of ground state energies and gradients and excitation properties. Key features of OpenMMPol include a modular architecture facilitating extensibility, parallel computing capabilities for enhanced performance on modern cluster architectures, a user-friendly interface for intuitive implementation, and a simple and flexible structure for providing input data. To show the capabilities offered by the library, we present an interface with PySCF to perform QM/AMOEBA molecular dynamics, geometry optimization, and excited-state calculation based on (time-dependent) density functional theory.
Collapse
Affiliation(s)
- Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Tommaso Nottoli
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Michele Nottoli
- Institute of Applied Analysis and Numerical Simulation, Universität Stuttgart, Pfaffenwaldring 57, D-70569 Stuttgart, Germany
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Filippo Lipparini
- 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
|
4
|
Salvadori G, Mazzeo P, Accomasso D, Cupellini L, Mennucci B. Deciphering Photoreceptors Through Atomistic Modeling from Light Absorption to Conformational Response. J Mol Biol 2024; 436:168358. [PMID: 37944793 DOI: 10.1016/j.jmb.2023.168358] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2023] [Revised: 10/28/2023] [Accepted: 11/02/2023] [Indexed: 11/12/2023]
Abstract
In this review, we discuss the successes and challenges of the atomistic modeling of photoreceptors. Throughout our presentation, we integrate explanations of the primary methodological approaches, ranging from quantum mechanical descriptions to classical enhanced sampling methods, all while providing illustrative examples of their practical application to specific systems. To enhance the effectiveness of our analysis, our primary focus has been directed towards the examination of applications across three distinct photoreceptors. These include an example of Blue Light-Using Flavin (BLUF) domains, a bacteriophytochrome, and the orange carotenoid protein (OCP) employed by cyanobacteria for photoprotection. Particular emphasis will be placed on the pivotal role played by the protein matrix in fine-tuning the initial photochemical event within the embedded chromophore. Furthermore, we will investigate how this localized perturbation initiates a cascade of events propagating from the binding pocket throughout the entire protein structure, thanks to the intricate network of interactions between the chromophore and the protein.
Collapse
Affiliation(s)
- Giacomo Salvadori
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Patrizia Mazzeo
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Davide Accomasso
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Lorenzo Cupellini
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| |
Collapse
|
5
|
Pes F, Polack É, Mazzeo P, Dusson G, Stamm B, Lipparini F. A Quasi Time-Reversible Scheme Based on Density Matrix Extrapolation on the Grassmann Manifold for Born-Oppenheimer Molecular Dynamics. J Phys Chem Lett 2023; 14:9720-9726. [PMID: 37879072 PMCID: PMC10626629 DOI: 10.1021/acs.jpclett.3c02098] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Accepted: 10/11/2023] [Indexed: 10/27/2023]
Abstract
This Letter introduces the so-called Quasi Time-Reversible scheme based on Grassmann extrapolation (QTR G-Ext) of density matrices for an accurate calculation of initial guesses in Born-Oppenheimer Molecular Dynamics (BOMD) simulations. The method shows excellent results on four large molecular systems that are representative of real-life production applications, ranging from 21 to 94 atoms simulated with Kohn-Sham (KS) density functional theory surrounded with a classical environment with 6k to 16k atoms. Namely, it clearly reduces the number of self-consistent field iterations while at the same time achieving energy-conserving simulations, resulting in a considerable speed-up of BOMD simulations even when tight convergence of the KS equations is required.
Collapse
Affiliation(s)
- Federica Pes
- Dipartimento
di Chimica e Chimica Industriale, Università
di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Étienne Polack
- CERMICS, École des Ponts and Inria Paris, 6 & 8 avenue Blaise Pascal, 77455 Marne-la-Valée, France
| | - Patrizia Mazzeo
- Dipartimento
di Chimica e Chimica Industriale, Università
di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Geneviève Dusson
- Laboratoire
de Mathématiques de Besançon, UMR CNRS 6623, Université de Franche-Comté, 16 route de Gray, 25030 Besançon, France
| | - Benjamin Stamm
- Institute
of Applied Analysis and Numerical Simulation, University of Stuttgart, 70569 Stuttgart, Germany
| | - Filippo Lipparini
- Dipartimento
di Chimica e Chimica Industriale, Università
di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
6
|
Nakata H, Fedorov DG. Analytic Gradient for Time-Dependent Density Functional Theory Combined with the Fragment Molecular Orbital Method. J Chem Theory Comput 2023; 19:1276-1285. [PMID: 36753486 DOI: 10.1021/acs.jctc.2c01177] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Abstract
The analytic energy gradient of energy with respect to nuclear coordinates is derived for the fragment molecular orbital (FMO) method combined with time-dependent density functional theory (TDDFT). The response terms arising from the use of a polarizable embedding are derived. The obtained analytic FMO-TDDFT gradient is shown to be accurate in comparison to both numerical FMO-TDDFT and unfragmented TDDFT gradients, at the level of two- and three-body expansions. The gradients are used for geometry optimizations, molecular dynamics, vibrational calculations, and simulations of IR and Raman spectra of excited states. The developed method is used to optimize the geometry of the ground and excited electronic states of the photoactive yellow protein (PDB: 2PHY).
Collapse
Affiliation(s)
- Hiroya Nakata
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
7
|
Mazzeo P, Hashem S, Lipparini F, Cupellini L, Mennucci B. Fast Method for Excited-State Dynamics in Complex Systems and Its Application to the Photoactivation of a Blue Light Using Flavin Photoreceptor. J Phys Chem Lett 2023; 14:1222-1229. [PMID: 36716231 PMCID: PMC9923743 DOI: 10.1021/acs.jpclett.2c03797] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2022] [Accepted: 01/20/2023] [Indexed: 06/18/2023]
Abstract
The excited-state dynamics of molecules embedded in complex (bio)matrices is still a challenging goal for quantum chemical models. Hybrid QM/MM models have proven to be an effective strategy, but an optimal combination of accuracy and computational cost still has to be found. Here, we present a method which combines the accuracy of a polarizable embedding QM/MM approach with the computational efficiency of an excited-state self-consistent field method. The newly implemented method is applied to the photoactivation of the blue-light-using flavin (BLUF) domain of the AppA protein. We show that the proton-coupled electron transfer (PCET) process suggested for other BLUF proteins is still valid also for AppA.
Collapse
|
8
|
Song C. State averaged CASSCF in AMOEBA polarizable water model for simulating nonadiabatic molecular dynamics with nonequilibrium solvation effects. J Chem Phys 2023; 158:014101. [PMID: 36610973 DOI: 10.1063/5.0131689] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
This paper presents a state-averaged complete active space self-consistent field (SA-CASSCF) in the atomic multipole optimized energetics for biomolecular application (AMOEBA) polarizable water model, which enables rigorous simulation of non-adiabatic molecular dynamics with nonequilibrium solvation effects. The molecular orbital and configuration interaction coefficients of the solute wavefunction, and the induced dipoles on solvent atoms, are solved by minimizing the state averaged energy variationally. In particular, by formulating AMOEBA water models and the polarizable continuum model (PCM) in a unified way, the algorithms developed for computing SA-CASSCF/PCM energies, analytical gradients, and non-adiabatic couplings in our previous work can be generalized to SA-CASSCF/AMOEBA by properly substituting a specific list of variables. Implementation of this method will be discussed with the emphasis on how the calculations of different terms are partitioned between the quantum chemistry and molecular mechanics codes. We will present and discuss results that demonstrate the accuracy and performance of the implementation. Next, we will discuss results that compare three solvent models that work with SA-CASSCF, i.e., PCM, fixed-charge force fields, and the newly implemented AMOEBA. Finally, the new SA-CASSCF/AMOEBA method has been interfaced with the ab initio multiple spawning method to carry out non-adiabatic molecular dynamics simulations. This method is demonstrated by simulating the photodynamics of the model retinal protonated Schiff base molecule in water.
Collapse
Affiliation(s)
- Chenchen Song
- Department of Chemistry, University of California Davis, Davis, California 95616, USA
| |
Collapse
|
9
|
Bondanza M, Demoulin B, Lipparini F, Barbatti M, Mennucci B. Trajectory Surface Hopping for a Polarizable Embedding QM/MM Formulation. J Phys Chem A 2022; 126:6780-6789. [PMID: 36107729 PMCID: PMC9527758 DOI: 10.1021/acs.jpca.2c04756] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
We present the implementation of trajectory surface-hopping
nonadiabatic
dynamics for a polarizable embedding QM/MM formulation. Time-dependent
density functional theory was used at the quantum mechanical level
of theory, whereas the molecular mechanics description involved the
polarizable AMOEBA force field. This implementation has been obtained
by integrating the surface-hopping program Newton-X NS with an interface
between the Gaussian 16 and the Tinker suites of codes to calculate
QM/AMOEBA energies and forces. The implementation has been tested
on a photoinduced electron-driven proton-transfer reaction involving
pyrimidine and a hydrogen-bonded water surrounded by a small cluster
of water molecules and within a large water droplet.
Collapse
Affiliation(s)
- Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | | | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Mario Barbatti
- Aix Marseille University, CNRS, ICR, 13385 Marseille, France
- Institut Universitaire de France, 75231 Paris, France
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
10
|
Nottoli M, Mazzeo P, Lipparini F, Cupellini L, Mennucci B. A ΔSCF model for excited states within a polarisable embedding. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2089605] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Pisa, Italy
| | - Patrizia Mazzeo
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Pisa, Italy
- Scuola Normale Superiore, Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Pisa, Italy
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Pisa, Italy
| |
Collapse
|
11
|
Song C. State-averaged CASSCF with polarizable continuum model for studying photoreactions in solvents: Energies, analytical nuclear gradients, and non-adiabatic couplings. J Chem Phys 2022; 156:104102. [DOI: 10.1063/5.0085855] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
This paper presents state-averaged complete active space self-consistent field in polarizable continuum model (PCM) for studies of photoreactions in solvents. The wavefunctions of the solute and the PCM surface charges of the solvent are optimized simultaneously such that the state-averaged free energy is variationally minimized. The method supports both fixed weights and dynamic weights where the weights are automatically adjusted based on the energy gaps. The corresponding analytical nuclear gradients and non-adiabatic couplings are also derived. Furthermore, we show how the new method can be entirely formulated in terms of seven basic operations, which allows the implementation to benefit from existing high-performance libraries on graphical processing units. Results demonstrating the accuracy and performance of the implementation are presented and discussed. We also apply the new method to the study of minimal conical intersection search and photoreaction energy pathways in solvents. Effects from the polarity of the solvents and different formulas of dynamic weights are compared and discussed.
Collapse
Affiliation(s)
- Chenchen Song
- Department of Chemistry, University of California Davis, Davis, California 95616, USA
| |
Collapse
|
12
|
Gómez S, Soysal EN, Worth GA. Micro-Solvated DMABN: Excited State Quantum Dynamics and Dual Fluorescence Spectra. Molecules 2021; 26:7247. [PMID: 34885829 PMCID: PMC8658867 DOI: 10.3390/molecules26237247] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2021] [Revised: 11/23/2021] [Accepted: 11/24/2021] [Indexed: 12/04/2022] Open
Abstract
In this work, we report a complete analysis by theoretical and spectroscopic methods of the short-time behaviour of 4-(dimethylamino)benzonitrile (DMABN) in the gas phase as well as in cyclohexane, tetrahydrofuran, acetonitrile, and water solution, after excitation to the La state. The spectroscopic properties of DMABN were investigated experimentally using UV absorption and fluorescence emission spectroscopy. The computational study was developed at different electronic structure levels and using the Polarisable Continuum Model (PCM) and explicit solvent molecules to reproduce the solvent environment. Additionally, excited state quantum dynamics simulations in the diabatic picture using the direct dynamics variational multiconfigurational Gaussian (DD-vMCG) method were performed, the largest quantum dynamics "on-the-fly" simulations performed with this method until now. The comparison with fully converged multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) dynamics on parametrised linear vibronic coupling (LVC) potentials show very similar population decays and evolution of the nuclear wavepacket. The ring C=C stretching and three methyl tilting modes are identified as the responsible motions for the internal conversion from the La to the Lb states. No major differences are observed in the ultrafast initial decay in different solvents, but we show that this effect depends strongly on the level of electronic structure used.
Collapse
Affiliation(s)
- Sandra Gómez
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK;
| | | | - Graham A. Worth
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK;
| |
Collapse
|
13
|
Polack É, Dusson G, Stamm B, Lipparini F. Grassmann Extrapolation of Density Matrices for Born-Oppenheimer Molecular Dynamics. J Chem Theory Comput 2021; 17:6965-6973. [PMID: 34623810 PMCID: PMC8582259 DOI: 10.1021/acs.jctc.1c00751] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2021] [Indexed: 01/16/2023]
Abstract
Born-Oppenheimer molecular dynamics (BOMD) is a powerful but expensive technique. The main bottleneck in a density functional theory BOMD calculation is the solution to the Kohn-Sham (KS) equations that requires an iterative procedure that starts from a guess for the density matrix. Converged densities from previous points in the trajectory can be used to extrapolate a new guess; however, the nonlinear constraint that an idempotent density needs to satisfy makes the direct use of standard linear extrapolation techniques not possible. In this contribution, we introduce a locally bijective map between the manifold where the density is defined and its tangent space so that linear extrapolation can be performed in a vector space while, at the same time, retaining the correct physical properties of the extrapolated density using molecular descriptors. We apply the method to real-life, multiscale, polarizable QM/MM BOMD simulations, showing that sizeable performance gains can be achieved, especially when a tighter convergence to the KS equations is required.
Collapse
Affiliation(s)
- Étienne Polack
- Laboratoire
de Mathématiques de Besançon, UMR CNRS 6623, Université
Bourgogne Franche-Comté, 16 Route de Gray, 25030 Besançon, France
| | - Geneviève Dusson
- Laboratoire
de Mathématiques de Besançon, UMR CNRS 6623, Université
Bourgogne Franche-Comté, 16 Route de Gray, 25030 Besançon, France
| | - Benjamin Stamm
- Department
of Mathematics, RWTH Aachen University, Schinkelstr. 2, 52062 Aachen, Germany
| | - Filippo Lipparini
- Dipartimento
di Chimica e Chimica Industriale, Univeristà
di Pisa, Via G. Moruzzi
13, I-56124 Pisa, Italy
| |
Collapse
|
14
|
Hashem S, Macaluso V, Nottoli M, Lipparini F, Cupellini L, Mennucci B. From crystallographic data to the solution structure of photoreceptors: the case of the AppA BLUF domain. Chem Sci 2021; 12:13331-13342. [PMID: 34777752 PMCID: PMC8528011 DOI: 10.1039/d1sc03000k] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2021] [Accepted: 09/01/2021] [Indexed: 12/28/2022] Open
Abstract
Photoreceptor proteins bind a chromophore, which, upon light absorption, modifies its geometry or its interactions with the protein, finally inducing the structural change needed to switch the protein from an inactive to an active or signaling state. In the Blue Light-Using Flavin (BLUF) family of photoreceptors, the chromophore is a flavin and the changes have been connected with a rearrangement of the hydrogen bond network around it on the basis of spectroscopic changes measured for the dark-to-light conversion. However, the exact conformational change triggered by the photoexcitation is still elusive mainly because a clear consensus on the identity not only of the light activated state but also of the dark one has not been achieved. Here, we present an integrated investigation that combines microsecond MD simulations starting from the two conflicting crystal structures available for the AppA BLUF domain with calculations of NMR, IR and UV-Vis spectra using a polarizable QM/MM approach. Thanks to such a combined analysis of the three different spectroscopic responses, a robust characterization of the structure of the dark state in solution is given together with the uncovering of important flaws of the most popular molecular mechanisms present in the literature for the dark-to-light activation.
Collapse
Affiliation(s)
- Shaima Hashem
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa Via G. Moruzzi 13 56124 Pisa Italy
| | - Veronica Macaluso
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa Via G. Moruzzi 13 56124 Pisa Italy
| | - Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa Via G. Moruzzi 13 56124 Pisa Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa Via G. Moruzzi 13 56124 Pisa Italy
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa Via G. Moruzzi 13 56124 Pisa Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa Via G. Moruzzi 13 56124 Pisa Italy
| |
Collapse
|
15
|
Macaluso V, Hashem S, Nottoli M, Lipparini F, Cupellini L, Mennucci B. Ultrafast Transient Infrared Spectroscopy of Photoreceptors with Polarizable QM/MM Dynamics. J Phys Chem B 2021; 125:10282-10292. [PMID: 34476939 PMCID: PMC8450903 DOI: 10.1021/acs.jpcb.1c05753] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2021] [Revised: 08/04/2021] [Indexed: 01/02/2023]
Abstract
Ultrafast transient infrared (TRIR) spectroscopy is widely used to measure the excitation-induced structural changes of protein-bound chromophores. Here, we design a novel and general strategy to compute TRIR spectra of photoreceptors by combining μs-long MM molecular dynamics with ps-long QM/AMOEBA Born-Oppenheimer molecular dynamics (BOMD) trajectories for both ground and excited electronic states. As a proof of concept, the strategy is here applied to AppA, a blue-light-utilizing flavin (BLUF) protein, found in bacteria. We first analyzed the short-time evolution of the embedded flavin upon excitation revealing that its dynamic Stokes shift is ultrafast and mainly driven by the internal reorganization of the chromophore. A different normal-mode representation was needed to describe ground- and excited-state IR spectra. In this way, we could assign all of the bands observed in the measured transient spectrum. In particular, we could characterize the flavin isoalloxazine-ring region of the spectrum, for which a full and clear description was missing.
Collapse
Affiliation(s)
| | | | - Michele Nottoli
- Dipartimento di Chimica e
Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Filippo Lipparini
- 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
|
16
|
Nottoli M, Nifosì R, Mennucci B, Lipparini F. Energy, Structures, and Response Properties with a Fully Coupled QM/AMOEBA/ddCOSMO Implementation. J Chem Theory Comput 2021; 17:5661-5672. [PMID: 34476941 PMCID: PMC8444335 DOI: 10.1021/acs.jctc.1c00555] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
![]()
We present the implementation
of a fully coupled polarizable QM/MM/continuum
model based on the AMOEBA polarizable force field and the domain decomposition
implementation of the conductor-like screening model. Energies, response
properties, and analytical gradients with respect to both QM and MM
nuclear positions are available, and a generic, atomistic cavity can
be employed. The model is linear scaling in memory requirements and
computational cost with respect to the number of classical atoms and
is therefore suited to model large, complex systems. Using three variants
of the green-fluorescent protein, we investigate the overall computational
cost of such calculations and the effect of the continuum model on
the convergence of the computed properties with respect to the size
of the embedding. We also demonstrate the fundamental role of polarization
effects by comparing polarizable and nonpolarizable embeddings to
fully QM ones.
Collapse
Affiliation(s)
- Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Riccardo Nifosì
- NEST, Istituto Nanoscienze-CNR and Scuola Normale Superiore, Piazza San Silvestro 12, I-56127 Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| |
Collapse
|
17
|
Nottoli M, Bondanza M, Lipparini F, Mennucci B. An enhanced sampling QM/AMOEBA approach: The case of the excited state intramolecular proton transfer in solvated 3-hydroxyflavone. J Chem Phys 2021; 154:184107. [PMID: 34241028 DOI: 10.1063/5.0046844] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
We present an extension of the polarizable quantum mechanical (QM)/AMOEBA approach to enhanced sampling techniques. This is achieved by connecting the enhanced sampling PLUMED library to the machinery based on the interface of Gaussian and Tinker to perform QM/AMOEBA molecular dynamics. As an application, we study the excited state intramolecular proton transfer of 3-hydroxyflavone in two solvents: methanol and methylcyclohexane. By using a combination of molecular dynamics and umbrella sampling, we find an ultrafast component of the transfer, which is common to the two solvents, and a much slower component, which is active in the protic solvent only. The mechanisms of the two components are explained in terms of intramolecular vibrational redistribution and intermolecular hydrogen-bonding, respectively. Ground and excited state free energies along an effective reaction coordinate are finally obtained allowing for a detailed analysis of the solvent mediated mechanism.
Collapse
Affiliation(s)
- Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
18
|
Nottoli M, Cupellini L, Lipparini F, Granucci G, Mennucci B. Multiscale Models for Light-Driven Processes. Annu Rev Phys Chem 2021; 72:489-513. [PMID: 33561359 DOI: 10.1146/annurev-physchem-090419-104031] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Multiscale models combining quantum mechanical and classical descriptions are a very popular strategy to simulate properties and processes of complex systems. Many alternative formulations have been developed, and they are now available in all of the most widely used quantum chemistry packages. Their application to the study of light-driven processes, however, is more recent, and some methodological and numerical problems have yet to be solved. This is especially the case for the polarizable formulation of these models, the recent advances in which we review here. Specifically, we identify and describe the most important specificities that the polarizable formulation introduces into both the simulation of excited-state dynamics and the modeling of excitation energy and electron transfer processes.
Collapse
Affiliation(s)
- Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Giovanni Granucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| |
Collapse
|