1
|
Sarkar R, Loos PF, Boggio-Pasqua M, Jacquemin D. Assessing the Performances of CASPT2 and NEVPT2 for Vertical Excitation Energies. J Chem Theory Comput 2022; 18:2418-2436. [PMID: 35333060 DOI: 10.1021/acs.jctc.1c01197] [Citation(s) in RCA: 25] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Methods able to simultaneously account for both static and dynamic electron correlations have often been employed, not only to model photochemical events but also to provide reference values for vertical transition energies, hence allowing benchmarking of lower-order models. In this category, both the complete-active-space second-order perturbation theory (CASPT2) and the N-electron valence state second-order perturbation theory (NEVPT2) are certainly popular, the latter presenting the advantage of not requiring the application of the empirical ionization-potential-electron-affinity (IPEA) and level shifts. However, the actual accuracy of these multiconfigurational approaches is not settled yet. In this context, to assess the performances of these approaches, the present work relies on highly accurate (±0.03 eV) aug-cc-pVTZ vertical transition energies for 284 excited states of diverse character (174 singlet, 110 triplet, 206 valence, 78 Rydberg, 78 n → π*, 119 π → π*, and 9 double excitations) determined in 35 small- to medium-sized organic molecules containing from three to six non-hydrogen atoms. The CASPT2 calculations are performed with and without IPEA shift and compared to the partially contracted (PC) and strongly contracted (SC) variants of NEVPT2. We find that both CASPT2 with IPEA shift and PC-NEVPT2 provide fairly reliable vertical transition energy estimates, with slight overestimations and mean absolute errors of 0.11 and 0.13 eV, respectively. These values are found to be rather uniform for the various subgroups of transitions. The present work completes our previous benchmarks focused on single-reference wave function methods ( J. Chem. Theory Comput. 2018, 14, 4360; J. Chem. Theory Comput. 2020, 16, 1711), hence allowing for a fair comparison between various families of electronic structure methods. In particular, we show that ADC(2), CCSD, and CASPT2 deliver similar accuracies for excited states with a dominant single-excitation character.
Collapse
Affiliation(s)
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, CNRS, UPS, Université de Toulouse, Toulouse 31062, France
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques, CNRS, UPS, Université de Toulouse, Toulouse 31062, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
2
|
Giuliani G, Melaccio F, Gozem S, Cappelli A, Olivucci M. QM/MM Investigation of the Spectroscopic Properties of the Fluorophore of Bacterial Luciferase. J Chem Theory Comput 2021; 17:605-613. [PMID: 33449693 DOI: 10.1021/acs.jctc.0c01078] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We employ replica-exchange molecular dynamics (REMD) and a hybrid ab initio multiconfigurational quantum mechanics/molecular mechanics (QM/MM) approach to model the absorption and fluorescence properties of bacterial luciferin-luciferase. Specifically, we employ complete active space perturbation theory (CASPT2) and study the effect of active space, basis set, and IPEA shift on the computed energies. We discuss the effect of the protein environment on the fluorophore's excited-state potential energy surface and the role that the protein plays in enhancing the fluorescence quantum yield in bacterial bioluminescence.
Collapse
Affiliation(s)
- Germano Giuliani
- Department of Biotechnology, Chemistry and Pharmacy, University of Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Federico Melaccio
- Department of Biotechnology, Chemistry and Pharmacy, University of Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Samer Gozem
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States
| | - Andrea Cappelli
- Department of Biotechnology, Chemistry and Pharmacy, University of Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Massimo Olivucci
- Department of Biotechnology, Chemistry and Pharmacy, University of Siena, Via A. Moro 2, 53100 Siena, Italy.,Department of Chemistry, Bowling Green State University, Bowing Green, Ohio 43403, United States
| |
Collapse
|
3
|
Nikolaev DM, Shtyrov AA, Mereshchenko AS, Panov MS, Tveryanovich YS, Ryazantsev MN. An assessment of water placement algorithms in quantum mechanics/molecular mechanics modeling: the case of rhodopsins' first spectral absorption band maxima. Phys Chem Chem Phys 2020; 22:18114-18123. [PMID: 32761024 DOI: 10.1039/d0cp02638g] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
Quantum mechanics/molecular mechanics (QM/MM) models are a widely used tool to obtain detailed insight into the properties and functioning of proteins. The outcome of QM/MM studies heavily depends on the quality of the applied QM/MM model. Prediction and right placement of internal water molecules in protein cavities is one of the critical parts of any QM/MM model construction. Herein, we performed a systematic study of four protein hydration algorithms. We tested these algorithms for their ability to predict X-ray-resolved water molecules for a set of membrane photosensitive rhodopsin proteins, as well as the influence of the applied water placement algorithms on the QM/MM calculated absorption maxima (λmax) of these proteins. We used 49 rhodopsins and their intermediates with available X-ray structures as the test set. We found that a proper choice of hydration algorithms and setups is needed to predict functionally important water molecules in the chromophore-binding cavity of rhodopsins, such as the water cluster in the N-H region of bacteriorhodopsin or two water molecules in the binding pocket of bovine visual rhodopsin. The QM/MM calculated λmax of rhodopsins is also quite sensitive to the applied protein hydration protocols. The best methodology allows obtaining an 18.0 nm average value for the absolute deviation of the calculated λmax from the experimental λmax. Although the major effect of water molecules on λmax originates from the water molecules located in the binding pocket, the water molecules outside the binding pocket also affect the calculated λmax mainly by causing a reorganization of the protein structure. The results reported in this study can be used for the evaluation and further development of hydration methodologies, in general, and rhodopsin QM/MM models, in particular.
Collapse
Affiliation(s)
- Dmitrii M Nikolaev
- Nanotechnology Research and Education Centre RAS, Saint Petersburg Academic University, 8/3 Khlopina Street, St. Petersburg 194021, Russia.
| | | | | | | | | | | |
Collapse
|
4
|
Marín MDC, De Vico L, Dong SS, Gagliardi L, Truhlar DG, Olivucci M. Assessment of MC-PDFT Excitation Energies for a Set of QM/MM Models of Rhodopsins. J Chem Theory Comput 2019; 15:1915-1923. [PMID: 30721054 PMCID: PMC7096677 DOI: 10.1021/acs.jctc.8b01069] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
A methodology for the automatic production of quantum mechanical/molecular mechanical (QM/MM) models of retinal-binding rhodopsin proteins and subsequent prediction of their spectroscopic properties has been proposed recently by some of the authors. The technology employed for the evaluation of the excitation energies is called Automatic Rhodopsin Modeling (ARM), and it involves the use of the complete active space self-consistent field (CASSCF) method followed by a multiconfiguration second-order perturbation theory (in particular, CASPT2) calculation of external correlation energies. Although it was shown that ARM is capable of successfully reproducing and predicting spectroscopic property trends in chromophore-embedding protein sets, practical applications of such technology are limited by the high computational costs of the multiconfiguration perturbation theory calculations. In the present work we benchmark the more affordable multiconfiguration pair-density functional theory (MC-PDFT) method whose accuracy has been recently validated for retinal chromophores in the gas phase, indicating that MC-PDFT could potentially be used to analyze large (e.g., few hundreds) sets of rhodopsin proteins. Here, we test this theory for a set of rhodopsin QM/MM models whose experimental absorption maxima (λ a max) have been measured. The results indicate that MC-PDFT may be employed to calculate λ a max values for this important class of photoresponsive proteins.
Collapse
Affiliation(s)
- María Del Carmen Marín
- Department of Biotechnologies, Chemistry and Pharmacy , University of Siena , 53100 Siena , Italy
- Chemistry Department , Bowling Green State University , Bowling Green , 43403 Ohio , United States
| | - Luca De Vico
- Department of Biotechnologies, Chemistry and Pharmacy , University of Siena , 53100 Siena , Italy
| | - Sijia S Dong
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , Minneapolis , Minnesota 55455-0431 , United States
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , Minneapolis , Minnesota 55455-0431 , United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , Minneapolis , Minnesota 55455-0431 , United States
| | - Massimo Olivucci
- Department of Biotechnologies, Chemistry and Pharmacy , University of Siena , 53100 Siena , Italy
- Chemistry Department , Bowling Green State University , Bowling Green , 43403 Ohio , United States
- USIAS Institut d'Études Avanceés , Université de Strasbourg , 67083 Strasbourg , France
| |
Collapse
|
5
|
Gil ES, de Araújo BB, Gonçalves PFB. CASPT2, CASSCF and non-adiabatic molecular dynamics (NAMD) studies on the low-lying electronic states of 1 H-1,2,3-triazole photolysis. Phys Chem Chem Phys 2019; 21:25809-25819. [DOI: 10.1039/c9cp04313f] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The photolysis mechanisms of 1H-1,2,3-triazole and 1H-1,2,3-benzotriazole were elucidated by employing multiconfigurational methods (CASSCF and CASPT2) and non-adiabatic molecular dynamics.
Collapse
Affiliation(s)
- Eduarda Sangiogo Gil
- Grupo de Química Teórica
- Universidade Federal do Rio Grande do Sul – Instituto de Química
- Porto Alegre
- Brazil
| | - Bruno Bercini de Araújo
- Grupo de Química Teórica
- Universidade Federal do Rio Grande do Sul – Instituto de Química
- Porto Alegre
- Brazil
| | - Paulo F. B. Gonçalves
- Grupo de Química Teórica
- Universidade Federal do Rio Grande do Sul – Instituto de Química
- Porto Alegre
- Brazil
| |
Collapse
|
6
|
Ai Y, Zhao C, Xing J, Liu Y, Wang Z, Jin J, Xia S, Cui G, Wang X. Excited-State Decay Pathways of Flavin Molecules in Five Redox Forms: The Role of Conical Intersections. J Phys Chem A 2018; 122:7954-7961. [DOI: 10.1021/acs.jpca.8b07582] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Yuejie Ai
- College of Environmental Science and Engineering, North China Electric Power University, Beijing 102206, P.R. China
| | - Chaofeng Zhao
- College of Environmental Science and Engineering, North China Electric Power University, Beijing 102206, P.R. China
| | - Jinlu Xing
- College of Environmental Science and Engineering, North China Electric Power University, Beijing 102206, P.R. China
| | - Yang Liu
- College of Environmental Science and Engineering, North China Electric Power University, Beijing 102206, P.R. China
| | - Zhangxia Wang
- College of Environmental Science and Engineering, North China Electric Power University, Beijing 102206, P.R. China
| | - Jiaren Jin
- College of Environmental Science and Engineering, North China Electric Power University, Beijing 102206, P.R. China
| | - Shuhua Xia
- College of Life and Environmental Science, Minzu University of China, Beijing 100081, P.R. China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P.R. China
| | - Xiangke Wang
- College of Environmental Science and Engineering, North China Electric Power University, Beijing 102206, P.R. China
| |
Collapse
|
7
|
Nikolaev D, Shtyrov AA, Panov MS, Jamal A, Chakchir OB, Kochemirovsky VA, Olivucci M, Ryazantsev MN. A Comparative Study of Modern Homology Modeling Algorithms for Rhodopsin Structure Prediction. ACS OMEGA 2018; 3:7555-7566. [PMID: 30087916 PMCID: PMC6068592 DOI: 10.1021/acsomega.8b00721] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/13/2018] [Accepted: 06/21/2018] [Indexed: 06/08/2023]
Abstract
Rhodopsins are seven α-helical membrane proteins that are of great importance in chemistry, biology, and modern biotechnology. Any in silico study on rhodopsin properties and functioning requires a high-quality three-dimensional structure. Due to particular difficulties with obtaining membrane protein structures from the experiment, in silico prediction of the three-dimensional rhodopsin structure based only on its primary sequence is an especially important task. For the last few years, significant progress was made in the field of protein structure prediction, especially for methods based on comparative modeling. However, the majority of this progress was made for soluble proteins and further investigations are needed to achieve similar progress for membrane proteins. In this paper, we evaluate the performance of modern protein structure prediction methodologies (implemented in the Medeller, I-TASSER, and Rosetta packages) for their ability to predict rhodopsin structures. Three widely used methodologies were considered: two general methodologies that are commonly applied to soluble proteins and a methodology that uses constraints that are specific for membrane proteins. The test pool consisted of 36 target-template pairs with different sequence similarities that was constructed on the basis of 24 experimental rhodopsin structures taken from the RCSB database. As a result, we showed that all three considered methodologies allow obtaining rhodopsin structures with the quality that is close to the crystallographic one (root mean square deviation (RMSD) of the predicted structure from the corresponding X-ray structure up to 1.5 Å) if the target-template sequence identity is higher than 40%. Moreover, all considered methodologies provided structures of average quality (RMSD < 4.0 Å) if the target-template sequence identity is higher than 20%. Such structures can be subsequently used for further investigation of molecular mechanisms of protein functioning and for the development of modern protein-based biotechnologies.
Collapse
Affiliation(s)
- Dmitrii
M. Nikolaev
- Nanotechnology
Research and Education Centre RAS, Saint-Petersburg
Academic University, 8/3 Khlopina Street, St. Petersburg 194021, Russia
| | - Andrey A. Shtyrov
- Nanotechnology
Research and Education Centre RAS, Saint-Petersburg
Academic University, 8/3 Khlopina Street, St. Petersburg 194021, Russia
| | - Maxim S. Panov
- Institute
of Chemistry, Saint Petersburg State University, 7/9 Universitetskaya emb., St. Petersburg 199034, Russia
| | - Adeel Jamal
- Department
of Chemical Engineering, Massachusetts Institute
of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Oleg B. Chakchir
- Nanotechnology
Research and Education Centre RAS, Saint-Petersburg
Academic University, 8/3 Khlopina Street, St. Petersburg 194021, Russia
| | - Vladimir A. Kochemirovsky
- Institute
of Chemistry, Saint Petersburg State University, 7/9 Universitetskaya emb., St. Petersburg 199034, Russia
| | - Massimo Olivucci
- Department
of Biotechnology, Chemistry and Pharmacy, Università di Siena, via A. Moro 2, Siena I-53100, Italy
| | - Mikhail N. Ryazantsev
- Institute
of Chemistry, Saint Petersburg State University, 7/9 Universitetskaya emb., St. Petersburg 199034, Russia
- Institute
of Macromolecular Compounds of the Russian Academy of Sciences, 31 Bolshoy pr., St. Petersburg 199004, Russia
| |
Collapse
|
8
|
Wolański Ł, Grabarek D, Andruniów T. Is the choice of a standard zeroth-order hamiltonian in CASPT2 ansatz optimal in calculations of excitation energies in protonated and unprotonated schiff bases of retinal? J Comput Chem 2018; 39:1470-1480. [PMID: 29635695 DOI: 10.1002/jcc.25217] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2017] [Revised: 03/09/2018] [Accepted: 03/10/2018] [Indexed: 11/07/2022]
Abstract
To account for systematic error of CASPT2 method empirical modification of the zeroth-order Hamiltonian with Ionization Potential-Electron Affinity (IPEA) shift was introduced. The optimized IPEA value (0.25 a.u.), called standard IPEA (S-IPEA), was recommended but due to its unsatisfactory performance in multiple metallic and organic compounds it has been questioned lately as a general parameter working properly for all molecules under CASPT2 study. As we are interested in Schiff bases of retinal, an important question emerging from this conflict of choice, to use or not to use S-IPEA, is whether the introduction of the modified zeroth-order Hamiltonian into CASPT2 ansatz does really improve their energetics. To achieve this goal, we assessed an impact of the IPEA shift value, in a range of 0-0.35 a.u., on vertical excitation energies to low-lying singlet states of two protonated (RPSBs) and two unprotonated (RSBs) Schiff bases of retinal for which experimental data in gas phase are available. In addition, an effect of geometry, basis set, and active space on computed VEEs is also reported. We find, that for these systems, the choice of S-IPEA significantly overestimates both S0 →S1 and S0 →S2 energies and the best theoretical estimate, in reference to the experimental data, is provided with either unmodified zeroth-order Hamiltonian or small value of the IPEA shift in a range of 0.05-0.15 a.u., depending on active space and basis set size, equilibrium geometry, and character of the excited state. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Łukasz Wolański
- Advanced Materials Engineering and Modelling Group, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| | - Dawid Grabarek
- Advanced Materials Engineering and Modelling Group, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| | - Tadeusz Andruniów
- Advanced Materials Engineering and Modelling Group, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, Wroclaw, 50-370, Poland
| |
Collapse
|
9
|
Dokukina I, Marian CM, Weingart O. New Perspectives on an Old Issue: A Comparative MS-CASPT2 and OM2-MRCI Study of Polyenes and Protonated Schiff Bases. Photochem Photobiol 2017; 93:1345-1355. [PMID: 28833170 DOI: 10.1111/php.12833] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2017] [Accepted: 08/13/2017] [Indexed: 11/29/2022]
Abstract
Polyenic systems are involved in light perception of numerous living organisms. Although a π-conjugated backbone is a common feature of all polyenes, their photophysics may vary. We provide a comparative quantum mechanical study of low-lying S1 and S2 excited states in short (3-5 double bonds) symmetric all-trans linear polyenes and corresponding protonated Schiff bases. In our investigation, we use the well-established ab initio multireference CASPT2 approach and benchmark the efficient semiempirical OM2-MRCI approach against it. For all protonated Schiff bases, MS-CASPT2 results in two distinct S1 minima with inverted and noninverted bond length pattern, respectively. We find that OM2-MRCI is a computationally affordable and reliable alternative to MS-CASPT2 for investigations of polyenic systems, particularly when highly demanding calculations (e.g. excited-state dynamics) need to be performed.
Collapse
Affiliation(s)
- Irina Dokukina
- Institut für Theoretische Chemie und Computerchemie, Heinrich-Heine Universität Düsseldorf, Düsseldorf, Germany
| | - Christel M Marian
- Institut für Theoretische Chemie und Computerchemie, Heinrich-Heine Universität Düsseldorf, Düsseldorf, Germany
| | - Oliver Weingart
- Institut für Theoretische Chemie und Computerchemie, Heinrich-Heine Universität Düsseldorf, Düsseldorf, Germany
| |
Collapse
|
10
|
Nikolaev DM, Emelyanov A, Boitsov VM, Panov MS, Ryazantsev MN. A voltage-dependent fluorescent indicator for optogenetic applications, archaerhodopsin-3: Structure and optical properties from in silico modeling. F1000Res 2017; 6:33. [PMID: 28435665 PMCID: PMC5381632 DOI: 10.12688/f1000research.10541.3] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 11/23/2017] [Indexed: 12/25/2022] Open
Abstract
It was demonstrated in recent studies that some rhodopsins can be used in optogenetics as fluorescent indicators of membrane voltage. One of the promising candidates for these applications is archaerhodopsin-3. While it has already shown encouraging results, there is still a large room for improvement. One of possible directions is increasing the intensity of the protein's fluorescent signal. Rational design of mutants with an improved signal is an important task, which requires both experimental and theoretical studies. Herein, we used a homology-based computational approach to predict the three-dimensional structure of archaerhodopsin-3, and a Quantum Mechanics/Molecular Mechanics (QM/MM) hybrid approach with high-level multireference ab initio methodology (SORCI+Q/AMBER) to model optical properties of this protein. We demonstrated that this methodology allows for reliable prediction of structure and spectral properties of archaerhodopsin-3. The results of this study can be utilized for computational molecular design of efficient fluorescent indicators of membrane voltage for modern optogenetics on the basis of archaerhodopsin-3.
Collapse
Affiliation(s)
- Dmitrii M Nikolaev
- Saint-Petersburg National Research Academic University of the Russian Academy of Science, St. Petersburg, Russian Federation
| | - Anton Emelyanov
- Saint-Petersburg National Research Academic University of the Russian Academy of Science, St. Petersburg, Russian Federation
| | - Vitaly M Boitsov
- Saint-Petersburg National Research Academic University of the Russian Academy of Science, St. Petersburg, Russian Federation
| | - Maxim S Panov
- Saint-Petersburg State University, St. Petersburg, Russian Federation
| | - Mikhail N Ryazantsev
- Saint-Petersburg State University, St. Petersburg, Russian Federation.,Saint-Petersburg Scientific Center of the Russian Academy of Sciences, St. Petersburg, Russian Federation
| |
Collapse
|
11
|
Melaccio F, del Carmen Marín M, Valentini A, Montisci F, Rinaldi S, Cherubini M, Yang X, Kato Y, Stenrup M, Orozco-Gonzalez Y, Ferré N, Luk HL, Kandori H, Olivucci M. Toward Automatic Rhodopsin Modeling as a Tool for High-Throughput Computational Photobiology. J Chem Theory Comput 2016; 12:6020-6034. [DOI: 10.1021/acs.jctc.6b00367] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Affiliation(s)
- Federico Melaccio
- Department
of Biotechnology, Chemistry e Pharmacy, Università di Siena, via A. Moro 2, I-53100 Siena, Italy
| | - María del Carmen Marín
- Department
of Biotechnology, Chemistry e Pharmacy, Università di Siena, via A. Moro 2, I-53100 Siena, Italy
- Department
of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
| | - Alessio Valentini
- Department
of Biotechnology, Chemistry e Pharmacy, Università di Siena, via A. Moro 2, I-53100 Siena, Italy
| | - Fabio Montisci
- Department
of Biotechnology, Chemistry e Pharmacy, Università di Siena, via A. Moro 2, I-53100 Siena, Italy
| | - Silvia Rinaldi
- Department
of Biotechnology, Chemistry e Pharmacy, Università di Siena, via A. Moro 2, I-53100 Siena, Italy
| | - Marco Cherubini
- Department
of Biotechnology, Chemistry e Pharmacy, Università di Siena, via A. Moro 2, I-53100 Siena, Italy
| | - Xuchun Yang
- Department
of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
| | - Yoshitaka Kato
- Department
of Frontier Materials, Graduate School of Engineering, Nagoya Institute of Technology, Gokiso-cho, Showa-ku, Nagoya, Japan
| | - Michael Stenrup
- Aix-Marseille Université, CNRS, ICR, 13284 Marseille, France
| | - Yoelvis Orozco-Gonzalez
- Department
of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
- Institut
de Physique et Chimie des Matériaux de Strasbourg, UMR 7504 Université de Strasbourg-CNRS, F-67034 Strasbourg, France
- USIAS
Institut d’Études Avancées, Université de Strasbourg, 5 allée du Général Rouvillois, F-67083 Strasbourg, France
| | - Nicolas Ferré
- Aix-Marseille Université, CNRS, ICR, 13284 Marseille, France
| | - Hoi Ling Luk
- Department
of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
| | - Hideki Kandori
- Department
of Frontier Materials, Graduate School of Engineering, Nagoya Institute of Technology, Gokiso-cho, Showa-ku, Nagoya, Japan
| | - Massimo Olivucci
- Department
of Biotechnology, Chemistry e Pharmacy, Università di Siena, via A. Moro 2, I-53100 Siena, Italy
- Department
of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
- Institut
de Physique et Chimie des Matériaux de Strasbourg, UMR 7504 Université de Strasbourg-CNRS, F-67034 Strasbourg, France
- USIAS
Institut d’Études Avancées, Université de Strasbourg, 5 allée du Général Rouvillois, F-67083 Strasbourg, France
| |
Collapse
|
12
|
Grabarek D, Walczak E, Andruniów T. Assessing the Accuracy of Various Ab Initio Methods for Geometries and Excitation Energies of Retinal Chromophore Minimal Model by Comparison with CASPT3 Results. J Chem Theory Comput 2016; 12:2346-56. [DOI: 10.1021/acs.jctc.6b00108] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Dawid Grabarek
- Advanced Materials Engineering
and Modelling Group, Wroclaw University of Technology, Wyb. Wyspianskiego
27, 50-370 Wroclaw, Poland
| | - Elżbieta Walczak
- Advanced Materials Engineering
and Modelling Group, Wroclaw University of Technology, Wyb. Wyspianskiego
27, 50-370 Wroclaw, Poland
| | - Tadeusz Andruniów
- Advanced Materials Engineering
and Modelling Group, Wroclaw University of Technology, Wyb. Wyspianskiego
27, 50-370 Wroclaw, Poland
| |
Collapse
|
13
|
McConnell S, McKenzie RH, Olsen S. Valence-bond non-equilibrium solvation model for a twisting monomethine cyanine. J Chem Phys 2015; 142:084502. [PMID: 25725740 DOI: 10.1063/1.4907758] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
We propose and analyze a two-state valence-bond model of non-equilibrium solvation effects on the excited-state twisting reaction of monomethine cyanines. Suppression of this reaction is thought responsible for environment-dependent fluorescence yield enhancement in these dyes. Fluorescence is quenched because twisting is accompanied via the formation of dark twisted intramolecular charge-transfer (TICT) states. For monomethine cyanines, where the ground state is a superposition of structures with different bond and charge localizations, there are two possible twisting pathways with different charge localizations in the excited state. For parameters corresponding to symmetric monomethines, the model predicts two low-energy twisting channels on the excited-state surface, which leads to a manifold of TICT states. For typical monomethines, twisting on the excited state surface will occur with a small barrier or no barrier. Changes in the solvation configuration can differentially stabilize TICT states in channels corresponding to different bonds, and that the position of a conical intersection between adiabatic states moves in response to solvation to stabilize either one channel or the other. There is a conical intersection seam that grows along the bottom of the excited-state potential with increasing solvent polarity. For monomethine cyanines with modest-sized terminal groups in moderately polar solution, the bottom of the excited-state potential surface is completely spanned by a conical intersection seam.
Collapse
Affiliation(s)
- Sean McConnell
- School of Mathematics and Physics, The University of Queensland, Brisbane 4072, Australia
| | - Ross H McKenzie
- School of Mathematics and Physics, The University of Queensland, Brisbane 4072, Australia
| | - Seth Olsen
- School of Mathematics and Physics, The University of Queensland, Brisbane 4072, Australia
| |
Collapse
|
14
|
Brunk E, Rothlisberger U. Mixed Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations of Biological Systems in Ground and Electronically Excited States. Chem Rev 2015; 115:6217-63. [PMID: 25880693 DOI: 10.1021/cr500628b] [Citation(s) in RCA: 299] [Impact Index Per Article: 33.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Elizabeth Brunk
- †Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland.,‡Joint BioEnergy Institute, Lawrence Berkeley National Laboratory, Emeryville, California 94618, United States
| | - Ursula Rothlisberger
- †Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland.,§National Competence Center of Research (NCCR) MARVEL-Materials' Revolution: Computational Design and Discovery of Novel Materials, 1015 Lausanne, Switzerland
| |
Collapse
|
15
|
Olsen S. Canonical-ensemble state-averaged complete active space self-consistent field (SA-CASSCF) strategy for problems with more diabatic than adiabatic states: Charge-bond resonance in monomethine cyanines. J Chem Phys 2015; 142:044116. [DOI: 10.1063/1.4904298] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
|
16
|
Schwabe T, Beerepoot MTP, Olsen JMH, Kongsted J. Analysis of computational models for an accurate study of electronic excitations in GFP. Phys Chem Chem Phys 2015; 17:2582-8. [DOI: 10.1039/c4cp04524f] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
The accuracy of PERI-CC2 is successfully assessed against RVS-CC2 for GFP model clusters and subsequently applied to a whole protein model.
Collapse
Affiliation(s)
- Tobias Schwabe
- Center for Bioinformatics and Physical Chemistry Institute
- D-22148 Hamburg
- Germany
| | - Maarten T. P. Beerepoot
- Centre for Theoretical and Computational Chemistry
- Department of Chemistry
- University of Tromsø – The Arctic University of Norway
- N-9037 Tromsø
- Norway
| | - 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
| | - Jacob Kongsted
- Department of Physics
- Chemistry and Pharmacy
- University of Southern Denmark
- DK-5230 Odense
- Denmark
| |
Collapse
|
17
|
Walczak E, Andruniów T. Impacts of retinal polyene (de)methylation on the photoisomerization mechanism and photon energy storage of rhodopsin. Phys Chem Chem Phys 2015; 17:17169-81. [DOI: 10.1039/c5cp01939g] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Similar to native rhodopsin, a two-mode space-saving isomerization mechanism drives the photoreaction in (de)methylated rhodopsin analogues.
Collapse
Affiliation(s)
- Elżbieta Walczak
- Department of Chemistry
- Wroclaw University of Technology
- 50-370 Wroclaw
- Poland
| | - Tadeusz Andruniów
- Department of Chemistry
- Wroclaw University of Technology
- 50-370 Wroclaw
- Poland
| |
Collapse
|
18
|
Olsen S. Locally-Excited (LE) versus Charge-Transfer (CT) Excited State Competition in a Series of Para-Substituted Neutral Green Fluorescent Protein (GFP) Chromophore Models. J Phys Chem B 2014; 119:2566-75. [DOI: 10.1021/jp508723d] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Seth Olsen
- School of Mathematics and
Physics, The University of Queensland, Brisbane, QLD 4072, Australia
| |
Collapse
|
19
|
Gozem S, Mirzakulova E, Schapiro I, Melaccio F, Glusac KD, Olivucci M. A Conical Intersection Controls the Deactivation of the Bacterial Luciferase Fluorophore. Angew Chem Int Ed Engl 2014; 53:9870-5. [DOI: 10.1002/anie.201404011] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2014] [Indexed: 01/24/2023]
|
20
|
Gozem S, Mirzakulova E, Schapiro I, Melaccio F, Glusac KD, Olivucci M. A Conical Intersection Controls the Deactivation of the Bacterial Luciferase Fluorophore. Angew Chem Int Ed Engl 2014. [DOI: 10.1002/ange.201404011] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
|
21
|
Gozem S, Melaccio F, Valentini A, Filatov M, Huix-Rotllant M, Ferré N, Frutos LM, Angeli C, Krylov AI, Granovsky AA, Lindh R, Olivucci M. Shape of Multireference, Equation-of-Motion Coupled-Cluster, and Density Functional Theory Potential Energy Surfaces at a Conical Intersection. J Chem Theory Comput 2014; 10:3074-84. [PMID: 26588278 DOI: 10.1021/ct500154k] [Citation(s) in RCA: 143] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
We report and characterize ground-state and excited-state potential energy profiles using a variety of electronic structure methods along a loop lying on the branching plane associated with a conical intersection (CI) of a reduced retinal model, the penta-2,4-dieniminium cation (PSB3). Whereas the performance of the equation-of-motion coupled-cluster, density functional theory, and multireference methods had been tested along the excited- and ground-state paths of PSB3 in our earlier work, the ability of these methods to correctly describe the potential energy surface shape along a CI branching plane has not yet been investigated. This is the focus of the present contribution. We find, in agreement with earlier studies by others, that standard time-dependent DFT (TDDFT) does not yield the correct two-dimensional (i.e., conical) crossing along the branching plane but rather a one-dimensional (i.e., linear) crossing along the same plane. The same type of behavior is found for SS-CASPT2(IPEA=0), SS-CASPT2(IPEA=0.25), spin-projected SF-TDDFT, EOM-SF-CCSD, and, finally, for the reference MRCISD+Q method. In contrast, we found that MRCISD, CASSCF, MS-CASPT2(IPEA=0), MS-CASPT2(IPEA=0.25), XMCQDPT2, QD-NEVPT2, non-spin-projected SF-TDDFT, and SI-SA-REKS yield the expected conical crossing. To assess the effect of the different crossing topologies (i.e., linear or conical) on the PSB3 photoisomerization efficiency, we discuss the results of 100 semiclassical trajectories computed by CASSCF and SS-CASPT2(IPEA=0.25) for a PSB3 derivative. We show that for the same initial conditions, the two methods yield similar dynamics leading to isomerization quantum yields that differ by only a few percent.
Collapse
Affiliation(s)
- Samer Gozem
- Department of Chemistry, Bowling Green State University , Bowling Green, Ohio 43403, United States
| | - Federico Melaccio
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena , via A. Moro 2, I-53100 Siena, Italy
| | - Alessio Valentini
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena , via A. Moro 2, I-53100 Siena, Italy.,Departamento de Química Física, Universidad de Alcalá , E-28871 Alcalá de Henares, Madrid, Spain
| | - Michael Filatov
- Institut für Physikalische und Theoretische Chemie, Universität Bonn , Beringstrasse 4, 53115 Bonn, Germany
| | - Miquel Huix-Rotllant
- Aix-Marseille Université, CNRS, Institut de Chimie Radicalaire , Marseille, France
| | - Nicolas Ferré
- Aix-Marseille Université, CNRS, Institut de Chimie Radicalaire , Marseille, France
| | - Luis Manuel Frutos
- Departamento de Química Física, Universidad de Alcalá , E-28871 Alcalá de Henares, Madrid, Spain
| | - Celestino Angeli
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università di Ferrara , via Fossato di Mortara 17, I-44121 Ferrara, Italy
| | - Anna I Krylov
- Department of Chemistry, University of Southern California , Los Angeles, California 90089-0482, United States
| | | | - Roland Lindh
- Department of Chemistry - Ångström, the Theoretical Chemistry Programme, POB 518, SE-751 20 Uppsala, Sweden
| | - Massimo Olivucci
- Department of Chemistry, Bowling Green State University , Bowling Green, Ohio 43403, United States.,Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena , via A. Moro 2, I-53100 Siena, Italy
| |
Collapse
|
22
|
Robb MA. In This Molecule There Must Be a Conical Intersection. ADVANCES IN PHYSICAL ORGANIC CHEMISTRY 2014. [DOI: 10.1016/b978-0-12-800256-8.00003-5] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/04/2022]
|
23
|
Abstract
Molecular dynamics simulations of biomolecules have matured into powerful tools of structural biology. In addition to the commonly used empirical force field potentials, quantum mechanical descriptions are gaining popularity for structure optimization and dynamic simulations of peptides and proteins. In this chapter, we introduce methodological developments such as the QM/MM framework and linear-scaling QM that make efficient calculations on large biomolecules possible. We identify the most common scenarios in which quantum descriptions of peptides and proteins are employed, such as structural refinement, force field development, treatment of unusual residues, and predicting spectroscopic and exited state properties. The benefits and shortcomings of QM potentials, in comparison to classical force fields, are discussed, with special emphasis on the sampling problems of protein conformational space. Finally, recent examples of QM/MM calculations in light-sensitive membrane proteins illustrate typical applications of the reviewed methods.
Collapse
Affiliation(s)
- Thomas Steinbrecher
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | | |
Collapse
|
24
|
Ab initio, density functional theory, and semi-empirical calculations. METHODS IN MOLECULAR BIOLOGY (CLIFTON, N.J.) 2012; 924:3-27. [PMID: 23034743 DOI: 10.1007/978-1-62703-017-5_1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Abstract
This chapter introduces the theory and applications of commonly used methods of electronic structure calculation, with particular emphasis on methods applicable for modelling biomolecular systems. This chapter is sectioned as follows. We start by presenting ab initio methods, followed by a treatment of density functional theory (DFT) and some recent advances in semi-empirical methods. Treatment of excited states as well as basis sets are also presented.
Collapse
|
25
|
Gozem S, Huntress M, Schapiro I, Lindh R, Granovsky AA, Angeli C, Olivucci M. Dynamic Electron Correlation Effects on the Ground State Potential Energy Surface of a Retinal Chromophore Model. J Chem Theory Comput 2012; 8:4069-80. [DOI: 10.1021/ct3003139] [Citation(s) in RCA: 129] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Samer Gozem
- Department of Chemistry, Bowling
Green State University, Bowling Green, Ohio 43403, United States
| | - Mark Huntress
- Department of Chemistry, Bowling
Green State University, Bowling Green, Ohio 43403, United States
| | - Igor Schapiro
- Department of Chemistry, Bowling
Green State University, Bowling Green, Ohio 43403, United States
| | - Roland Lindh
- Department of Chemistry −
Ångström, the Theoretical Chemistry Programme, POB 518,
SE-751 20 Uppsala, Sweden
| | | | - Celestino Angeli
- Dipartimento di
Chimica, Università
di Ferrara, via Borsari 46, I-44121 Ferrara, Italy
| | - Massimo Olivucci
- Department of Chemistry, Bowling
Green State University, Bowling Green, Ohio 43403, United States
- Dipartimento di Chimica, Università
di Siena, via De Gasperi 2, I-53100 Siena, Italy
| |
Collapse
|
26
|
Ryazantsev MN, Altun A, Morokuma K. Color Tuning in rhodopsins: the origin of the spectral shift between the chloride-bound and anion-free forms of halorhodopsin. J Am Chem Soc 2012; 134:5520-3. [PMID: 22397521 DOI: 10.1021/ja3009117] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Detailed knowledge of the molecular mechanisms that control the spectral properties in the rhodopsin protein family is important for understanding the functions of these photoreceptors and for the rational design of artificial photosensitive proteins. Here we used a high-level ab initio QM/MM method to investigate the mechanism of spectral tuning in the chloride-bound and anion-free forms of halorhodopsin from Natronobacterium pharaonis (phR) and the interprotein spectral shift between them. We demonstrate that the chloride ion tunes the spectral properties of phR via two distinct mechanisms: (i) electrostatic interaction with the chromophore, which results in a 95 nm difference between the absorption maxima of the two forms, and (ii) induction of a structural reorganization in the protein, which changes the positions of charged and polar residues and reduces this difference to 29 nm. The present study expands our knowledge concerning the role of the reorganization of the internal H-bond network for color tuning in general and provides a detailed investigation of the tuning mechanism in phR in particular.
Collapse
Affiliation(s)
- Mikhail N Ryazantsev
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| | | | | |
Collapse
|
27
|
Roca-Sanjuán D, Aquilante F, Lindh R. Multiconfiguration second-order perturbation theory approach to strong electron correlation in chemistry and photochemistry. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2011. [DOI: 10.1002/wcms.97] [Citation(s) in RCA: 124] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
28
|
González L, Escudero D, Serrano-Andrés L. Progress and Challenges in the Calculation of Electronic Excited States. Chemphyschem 2011; 13:28-51. [DOI: 10.1002/cphc.201100200] [Citation(s) in RCA: 310] [Impact Index Per Article: 23.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2011] [Revised: 08/05/2011] [Indexed: 11/09/2022]
|
29
|
Jonasson G, Teuler JM, Vallverdu G, Mérola F, Ridard J, Lévy B, Demachy I. Excited State Dynamics of the Green Fluorescent Protein on the Nanosecond Time Scale. J Chem Theory Comput 2011; 7:1990-7. [DOI: 10.1021/ct200150r] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Gabriella Jonasson
- Laboratoire de Chimie-Physique, UMR8000, CNRS, Université Paris-Sud 11, 91405 Orsay Cedex, France
| | - Jean-Marie Teuler
- Laboratoire de Chimie-Physique, UMR8000, CNRS, Université Paris-Sud 11, 91405 Orsay Cedex, France
| | - Germain Vallverdu
- Equipe de Chimie-Physique, IPREM, UMR5254, Université de Pau et des Pays de l′Adour, 64000 Pau, France
| | - Fabienne Mérola
- Laboratoire de Chimie-Physique, UMR8000, CNRS, Université Paris-Sud 11, 91405 Orsay Cedex, France
| | - Jacqueline Ridard
- Laboratoire de Chimie-Physique, UMR8000, CNRS, Université Paris-Sud 11, 91405 Orsay Cedex, France
| | - Bernard Lévy
- Laboratoire de Chimie-Physique, UMR8000, CNRS, Université Paris-Sud 11, 91405 Orsay Cedex, France
| | - Isabelle Demachy
- Laboratoire de Chimie-Physique, UMR8000, CNRS, Université Paris-Sud 11, 91405 Orsay Cedex, France
| |
Collapse
|
30
|
Schapiro I, Melaccio F, Laricheva EN, Olivucci M. Using the computer to understand the chemistry of conical intersections. Photochem Photobiol Sci 2011; 10:867-86. [DOI: 10.1039/c0pp00290a] [Citation(s) in RCA: 51] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|