1
|
Li Manni G, Fdez. Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry. J Chem Theory Comput 2023; 19:6933-6991. [PMID: 37216210 PMCID: PMC10601490 DOI: 10.1021/acs.jctc.3c00182] [Citation(s) in RCA: 42] [Impact Index Per Article: 42.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Indexed: 05/24/2023]
Abstract
The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.
Collapse
Affiliation(s)
- Giovanni Li Manni
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Ignacio Fdez. Galván
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Ali Alavi
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
- Yusuf Hamied
Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Flavia Aleotti
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Francesco Aquilante
- Theory and
Simulation of Materials (THEOS) and National Centre for Computational
Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jochen Autschbach
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
| | - Davide Avagliano
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Alberto Baiardi
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jie J. Bao
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Stefano Battaglia
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Letitia Birnoschi
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | - Alejandro Blanco-González
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - Sergey I. Bokarev
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
- Chemistry
Department, School of Natural Sciences, Technical University of Munich, Lichtenbergstr. 4, 85748 Garching, Germany
| | - Ria Broer
- Theoretical
Chemistry, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747AG Groningen, The Netherlands
| | - Roberto Cacciari
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Paul B. Calio
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Rebecca K. Carlson
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Rafael Carvalho Couto
- Division
of Theoretical Chemistry and Biology, School of Engineering Sciences
in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Luis Cerdán
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
- Instituto
de Óptica (IO−CSIC), Consejo
Superior de Investigaciones Científicas, 28006, Madrid, Spain
| | - Liviu F. Chibotaru
- Department
of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Nicholas F. Chilton
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | | | - Irene Conti
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Sonia Coriani
- Department
of Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, 2800 Kongens Lyngby, Denmark
| | - Juliana Cuéllar-Zuquin
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Razan E. Daoud
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Nike Dattani
- HPQC Labs, Waterloo, N2T 2K9 Ontario Canada
- HPQC College, Waterloo, N2T 2K9 Ontario Canada
| | - Piero Decleva
- Istituto
Officina dei Materiali IOM-CNR and Dipartimento di Scienze Chimiche
e Farmaceutiche, Università degli
Studi di Trieste, I-34121 Trieste, Italy
| | - Coen de Graaf
- Department
of Physical and Inorganic Chemistry, Universitat
Rovira i Virgili, Tarragona 43007, Spain
- ICREA, Pg. Lluís
Companys 23, 08010 Barcelona, Spain
| | - Mickaël
G. Delcey
- Division
of Theoretical Chemistry and Biology, School of Engineering Sciences
in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Luca De Vico
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Werner Dobrautz
- Chalmers
University of Technology, Department of Chemistry
and Chemical Engineering, 41296 Gothenburg, Sweden
| | - Sijia S. Dong
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry and Chemical Biology, Department of Physics, and Department
of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Rulin Feng
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
- Department
of Chemistry, Fudan University, Shanghai 200433, China
| | - Nicolas Ferré
- Institut
de Chimie Radicalaire (UMR-7273), Aix-Marseille
Univ, CNRS, ICR 13013 Marseille, France
| | | | - Laura Gagliardi
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Marco Garavelli
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Yafu Guan
- State Key
Laboratory of Molecular Reaction Dynamics and Center for Theoretical
Computational Chemistry, Dalian Institute
of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| | - Meiyuan Guo
- SSRL, SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Matthew R. Hennefarth
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Matthew R. Hermes
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Chad E. Hoyer
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Miquel Huix-Rotllant
- Institut
de Chimie Radicalaire (UMR-7273), Aix-Marseille
Univ, CNRS, ICR 13013 Marseille, France
| | - Vishal Kumar Jaiswal
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Andy Kaiser
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Danil S. Kaliakin
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - Marjan Khamesian
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Daniel S. King
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Vladislav Kochetov
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Marek Krośnicki
- Institute
of Theoretical Physics and Astrophysics, Faculty of Mathematics, Physics
and Informatics, University of Gdańsk, ul Wita Stwosza 57, 80-952, Gdańsk, Poland
| | | | - Ernst D. Larsson
- Division
of Theoretical Chemistry, Chemical Centre, Lund University, P.O. Box 124, SE-22100, Lund, Sweden
| | - Susi Lehtola
- Molecular
Sciences Software Institute, Blacksburg, Virginia 24061, United States
- Department
of Chemistry, University of Helsinki, P.O. Box 55, FI-00014 University of Helsinki, Finland
| | - Marie-Bernadette Lepetit
- Condensed
Matter Theory Group, Institut Néel, CNRS UPR 2940, 38042 Grenoble, France
- Theory
Group, Institut Laue Langevin, 38042 Grenoble, France
| | - Hans Lischka
- Department
of Chemistry and Biochemistry, Texas Tech
University, Lubbock, Texas 79409-1061, United States
| | - Pablo López Ríos
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Marcus Lundberg
- Department
of Chemistry − Ångström Laboratory, Uppsala University, SE-75120 Uppsala, Sweden
| | - Dongxia Ma
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | | | - Francesco Montorsi
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Maximilian Mörchen
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Artur Nenov
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Vu Ha Anh Nguyen
- Department
of Chemistry, National University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Yoshio Nishimoto
- Graduate
School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - Meagan S. Oakley
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Massimo Olivucci
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Markus Oppel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Daniele Padula
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Riddhish Pandharkar
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Quan Manh Phung
- Department
of Chemistry, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
- Institute
of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8601, Japan
| | - Felix Plasser
- Department
of Chemistry, Loughborough University, Loughborough, LE11 3TU, U.K.
| | - Gerardo Raggi
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
- Quantum
Materials and Software LTD, 128 City Road, London, EC1V 2NX, United Kingdom
| | - Elisa Rebolini
- Scientific
Computing Group, Institut Laue Langevin, 38042 Grenoble, France
| | - Markus Reiher
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Ivan Rivalta
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Daniel Roca-Sanjuán
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Thies Romig
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Arta Anushirwan Safari
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Aitor Sánchez-Mansilla
- Department
of Physical and Inorganic Chemistry, Universitat
Rovira i Virgili, Tarragona 43007, Spain
| | - Andrew M. Sand
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry and Biochemistry, Butler University, Indianapolis, Indiana 46208, United States
| | - Igor Schapiro
- Institute
of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Thais R. Scott
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| | - Javier Segarra-Martí
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Francesco Segatta
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Dumitru-Claudiu Sergentu
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
- Laboratory
RA-03, RECENT AIR, A. I. Cuza University of Iaşi, RA-03 Laboratory (RECENT AIR), Iaşi 700506, Romania
| | - Prachi Sharma
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Ron Shepard
- Chemical
Sciences and Engineering Division, Argonne
National Laboratory, Lemont, Illinois 60439, USA
| | - Yinan Shu
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Jakob K. Staab
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | - Tjerk P. Straatsma
- National
Center for Computational Sciences, Oak Ridge
National Laboratory, Oak Ridge, Tennessee 37831-6373, United States
- Department
of Chemistry and Biochemistry, University
of Alabama, Tuscaloosa, Alabama 35487-0336, United States
| | | | - Bruno Nunes Cabral Tenorio
- Department
of Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, 2800 Kongens Lyngby, Denmark
| | - Donald G. Truhlar
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Liviu Ungur
- Department
of Chemistry, National University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Morgane Vacher
- Nantes
Université, CNRS, CEISAM, UMR 6230, F-44000 Nantes, France
| | - Valera Veryazov
- Division
of Theoretical Chemistry, Chemical Centre, Lund University, P.O. Box 124, SE-22100, Lund, Sweden
| | - Torben Arne Voß
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Oskar Weser
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Dihua Wu
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Xuchun Yang
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - David Yarkony
- Department
of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Chen Zhou
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - J. Patrick Zobel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Roland Lindh
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
- Uppsala
Center for Computational Chemistry (UC3), Uppsala University, PO Box 576, SE-751 23 Uppsala. Sweden
| |
Collapse
|
2
|
Kaiser A, Daoud RE, Aquilante F, Kühn O, De Vico L, Bokarev SI. A Multiconfigurational Wave Function Implementation of the Frenkel Exciton Model for Molecular Aggregates. J Chem Theory Comput 2023; 19:2918-2928. [PMID: 37115036 DOI: 10.1021/acs.jctc.3c00185] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
We present an implementation of the Frenkel exciton model into the OpenMolcas program package enabling calculations of collective electronic excited states of molecular aggregates based on a multiconfigurational wave function description of the individual monomers. The computational protocol avoids using diabatization schemes and, thus, supermolecule calculations. Additionally, the use of the Cholesky decomposition of the two-electron integrals entering pair interactions enhances the efficiency of the computational scheme. The application of the method is exemplified for two test systems, that is, a formaldehyde oxime and a bacteriochlorophyll-like dimer. For the sake of comparison with the dipole approximation, we restrict our considerations to situations where intermonomer exchange can be neglected. The protocol is expected to be beneficial for aggregates composed of molecules with extended π systems, unpaired electrons such as radicals or transition metal centers, where it should outperform widely used methods based on time-dependent density functional theory.
Collapse
Affiliation(s)
- Andy Kaiser
- Institut für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock. Germany
| | - Razan E Daoud
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Francesco Aquilante
- Theory and Simulation of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Oliver Kühn
- Institut für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock. Germany
| | - Luca De Vico
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Sergey I Bokarev
- Institut für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock. Germany
- Chemistry Department, Technical University of Munich, Lichtenbergstr. 4, 85748 Garching, Germany
| |
Collapse
|
3
|
Pedraza-González L, Barneschi L, Marszałek M, Padula D, De Vico L, Olivucci M. Automated QM/MM Screening of Rhodopsin Variants with Enhanced Fluorescence. J Chem Theory Comput 2023; 19:293-310. [PMID: 36516450 DOI: 10.1021/acs.jctc.2c00928] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
We present a computational protocol for the fast and automated screening of excited-state hybrid quantum mechanics/molecular mechanics (QM/MM) models of rhodopsins to be used as fluorescent probes based on the automatic rhodopsin modeling protocol (a-ARM). Such "a-ARM fluorescence screening protocol" is implemented through a general Python-based driver, PyARM, that is also proposed here. The implementation and performance of the protocol are benchmarked using different sets of rhodopsin variants whose absorption and, more relevantly, emission spectra have been experimentally measured. We show that, despite important limitations that make unsafe to use it as a black-box tool, the protocol reproduces the observed trends in fluorescence and it is capable of selecting novel potentially fluorescent rhodopsins. We also show that the protocol can be used in mechanistic investigations to discern fluorescence enhancement effects associated with a near degeneracy of the S1/S2 states or, alternatively, with a barrier generated via coupling of the S0/S1 wave functions.
Collapse
Affiliation(s)
- Laura Pedraza-González
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Via A. Moro 2, I-53100 Siena, Italy
| | - Leonardo Barneschi
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Via A. Moro 2, I-53100 Siena, Italy
| | - Michał Marszałek
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Via A. Moro 2, I-53100 Siena, Italy.,Faculty of Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiaǹskiego 27, 50-370 Wrocław, Poland
| | - Daniele Padula
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Via A. Moro 2, I-53100 Siena, Italy
| | - Luca De Vico
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Via A. Moro 2, I-53100 Siena, Italy
| | - Massimo Olivucci
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Via A. Moro 2, I-53100 Siena, Italy.,Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
| |
Collapse
|
4
|
Barneschi L, Marsili E, Pedraza-González L, Padula D, De Vico L, Kaliakin D, Blanco-González A, Ferré N, Huix-Rotllant M, Filatov M, Olivucci M. On the fluorescence enhancement of arch neuronal optogenetic reporters. Nat Commun 2022; 13:6432. [PMID: 36307417 PMCID: PMC9616920 DOI: 10.1038/s41467-022-33993-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2022] [Accepted: 10/07/2022] [Indexed: 12/25/2022] Open
Abstract
The lack of a theory capable of connecting the amino acid sequence of a light-absorbing protein with its fluorescence brightness is hampering the development of tools for understanding neuronal communications. Here we demonstrate that a theory can be established by constructing quantum chemical models of a set of Archaerhodopsin reporters in their electronically excited state. We found that the experimentally observed increase in fluorescence quantum yield is proportional to the computed decrease in energy difference between the fluorescent state and a nearby photoisomerization channel leading to an exotic diradical of the protein chromophore. This finding will ultimately support the development of technologies for searching novel fluorescent rhodopsin variants and unveil electrostatic changes that make light emission brighter and brighter.
Collapse
Affiliation(s)
- Leonardo Barneschi
- grid.9024.f0000 0004 1757 4641Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro 2, I-53100 Siena, Italy
| | - Emanuele Marsili
- grid.9024.f0000 0004 1757 4641Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro 2, I-53100 Siena, Italy ,grid.8250.f0000 0000 8700 0572University of Durham, Department of Chemistry, South Road, Durham, DH1 3LE United Kingdom ,grid.5337.20000 0004 1936 7603Present Address: Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, United Kingdom
| | - Laura Pedraza-González
- grid.9024.f0000 0004 1757 4641Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro 2, I-53100 Siena, Italy ,grid.5395.a0000 0004 1757 3729Present Address: Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe Moruzzi, 13, I-56124 Pisa, Italy
| | - Daniele Padula
- grid.9024.f0000 0004 1757 4641Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro 2, I-53100 Siena, Italy
| | - Luca De Vico
- grid.9024.f0000 0004 1757 4641Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro 2, I-53100 Siena, Italy
| | - Danil Kaliakin
- grid.253248.a0000 0001 0661 0035Department of Chemistry, Bowling Green State University, Bowling Green, OH 43403 USA
| | - Alejandro Blanco-González
- grid.253248.a0000 0001 0661 0035Department of Chemistry, Bowling Green State University, Bowling Green, OH 43403 USA
| | - Nicolas Ferré
- grid.462456.70000 0004 4902 8637Institut de Chimie Radicalaire (UMR-7273), Aix-Marseille Université, CNRS, 13397 Marseille, Cedex 20 France
| | - Miquel Huix-Rotllant
- grid.462456.70000 0004 4902 8637Institut de Chimie Radicalaire (UMR-7273), Aix-Marseille Université, CNRS, 13397 Marseille, Cedex 20 France
| | - Michael Filatov
- grid.258803.40000 0001 0661 1556Department of Chemistry, Kyungpook National University, Daegu, 702-701 South Korea
| | - Massimo Olivucci
- grid.9024.f0000 0004 1757 4641Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro 2, I-53100 Siena, Italy ,grid.253248.a0000 0001 0661 0035Department of Chemistry, Bowling Green State University, Bowling Green, OH 43403 USA ,grid.11843.3f0000 0001 2157 9291University of Strasbourg Institute for Advanced Studies, 5, alleé duGeń eŕ al Rouvillois, F-67083 Strasbourg, France
| |
Collapse
|
5
|
Pedraza-González L, Barneschi L, Padula D, De Vico L, Olivucci M. Evolution of the Automatic Rhodopsin Modeling (ARM) Protocol. Top Curr Chem (Cham) 2022; 380:21. [PMID: 35291019 PMCID: PMC8924150 DOI: 10.1007/s41061-022-00374-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2021] [Accepted: 01/29/2022] [Indexed: 10/27/2022]
Abstract
In recent years, photoactive proteins such as rhodopsins have become a common target for cutting-edge research in the field of optogenetics. Alongside wet-lab research, computational methods are also developing rapidly to provide the necessary tools to analyze and rationalize experimental results and, most of all, drive the design of novel systems. The Automatic Rhodopsin Modeling (ARM) protocol is focused on providing exactly the necessary computational tools to study rhodopsins, those being either natural or resulting from mutations. The code has evolved along the years to finally provide results that are reproducible by any user, accurate and reliable so as to replicate experimental trends. Furthermore, the code is efficient in terms of necessary computing resources and time, and scalable in terms of both number of concurrent calculations as well as features. In this review, we will show how the code underlying ARM achieved each of these properties.
Collapse
Affiliation(s)
- Laura Pedraza-González
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via Aldo Moro 2, 53100, Siena, Italy. .,Department of Chemistry and Industrial Chemistry, University of Pisa, Via Moruzzi 13, 56124, Pisa, Italy.
| | - Leonardo Barneschi
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via Aldo Moro 2, 53100, Siena, Italy
| | - Daniele Padula
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via Aldo Moro 2, 53100, Siena, Italy
| | - Luca De Vico
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via Aldo Moro 2, 53100, Siena, Italy.
| | - Massimo Olivucci
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via Aldo Moro 2, 53100, Siena, Italy. .,Department of Chemistry, Bowling Green State University, Bowling Green, OH, 43403, USA.
| |
Collapse
|
6
|
Mroginski MA, Adam S, Amoyal GS, Barnoy A, Bondar AN, Borin VA, Church JR, Domratcheva T, Ensing B, Fanelli F, Ferré N, Filiba O, Pedraza-González L, González R, González-Espinoza CE, Kar RK, Kemmler L, Kim SS, Kongsted J, Krylov AI, Lahav Y, Lazaratos M, NasserEddin Q, Navizet I, Nemukhin A, Olivucci M, Olsen JMH, Pérez de Alba Ortíz A, Pieri E, Rao AG, Rhee YM, Ricardi N, Sen S, Solov'yov IA, De Vico L, Wesolowski TA, Wiebeler C, Yang X, Schapiro I. Frontiers in Multiscale Modeling of Photoreceptor Proteins. Photochem Photobiol 2021; 97:243-269. [PMID: 33369749 DOI: 10.1111/php.13372] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Accepted: 11/01/2020] [Indexed: 02/06/2023]
Abstract
This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.
Collapse
Affiliation(s)
| | - Suliman Adam
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Gil S Amoyal
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Avishai Barnoy
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Ana-Nicoleta Bondar
- Freie Universität Berlin, Department of Physics, Theoretical Molecular Biophysics Group, Berlin, Germany
| | - Veniamin A Borin
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Jonathan R Church
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Tatiana Domratcheva
- Department of Chemistry, Lomonosov Moscow State University, Moscow, Russia.,Department Biomolecular Mechanisms, Max Planck Institute for Medical Research, Heidelberg, Germany
| | - Bernd Ensing
- Van 't Hoff Institute for Molecular Science and Amsterdam Center for Multiscale Modeling, University of Amsterdam, Amsterdam, The Netherlands
| | - Francesca Fanelli
- Department of Life Sciences, Center for Neuroscience and Neurotechnology, Università degli Studi di Modena e Reggio Emilia, Modena, Italy
| | | | - Ofer Filiba
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Laura Pedraza-González
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Siena, Italy
| | - Ronald González
- Institut für Chemie, Technische Universität Berlin, Berlin, Germany
| | | | - Rajiv K Kar
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Lukas Kemmler
- Freie Universität Berlin, Department of Physics, Theoretical Molecular Biophysics Group, Berlin, Germany
| | - Seung Soo Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon, Korea
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Odense, Denmark
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles, CA, USA
| | - Yigal Lahav
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel.,MIGAL - Galilee Research Institute, S. Industrial Zone, Kiryat Shmona, Israel
| | - Michalis Lazaratos
- Freie Universität Berlin, Department of Physics, Theoretical Molecular Biophysics Group, Berlin, Germany
| | - Qays NasserEddin
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Isabelle Navizet
- MSME, Univ Gustave Eiffel, CNRS UMR 8208, Univ Paris Est Creteil, Marne-la-Vallée, France
| | - Alexander Nemukhin
- Department of Chemistry, Lomonosov Moscow State University, Moscow, Russia.,Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Moscow, Russia
| | - Massimo Olivucci
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Siena, Italy.,Chemistry Department, Bowling Green State University, Bowling Green, OH, USA
| | - Jógvan Magnus Haugaard Olsen
- Department of Chemistry, Aarhus University, Aarhus, Denmark.,Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, Tromsø, Norway
| | - Alberto Pérez de Alba Ortíz
- Van 't Hoff Institute for Molecular Science and Amsterdam Center for Multiscale Modeling, University of Amsterdam, Amsterdam, The Netherlands
| | - Elisa Pieri
- Aix-Marseille Univ, CNRS, ICR, Marseille, France
| | - Aditya G Rao
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon, Korea
| | - Niccolò Ricardi
- Département de Chimie Physique, Université de Genève, Genève, Switzerland
| | - Saumik Sen
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Ilia A Solov'yov
- Department of Physics, Carl von Ossietzky University Oldenburg, Oldenburg, Germany
| | - Luca De Vico
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Siena, Italy
| | | | - Christian Wiebeler
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Xuchun Yang
- Chemistry Department, Bowling Green State University, Bowling Green, OH, USA
| | - Igor Schapiro
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| |
Collapse
|
7
|
Aquilante F, Autschbach J, Baiardi A, Battaglia S, Borin VA, Chibotaru LF, Conti I, De Vico L, Delcey M, Fdez Galván I, Ferré N, Freitag L, Garavelli M, Gong X, Knecht S, Larsson ED, Lindh R, Lundberg M, Malmqvist PÅ, Nenov A, Norell J, Odelius M, Olivucci M, Pedersen TB, Pedraza-González L, Phung QM, Pierloot K, Reiher M, Schapiro I, Segarra-Martí J, Segatta F, Seijo L, Sen S, Sergentu DC, Stein CJ, Ungur L, Vacher M, Valentini A, Veryazov V. Modern quantum chemistry with [Open]Molcas. J Chem Phys 2020; 152:214117. [PMID: 32505150 DOI: 10.1063/5.0004835] [Citation(s) in RCA: 212] [Impact Index Per Article: 53.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.
Collapse
Affiliation(s)
- Francesco Aquilante
- Theory and Simulation of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, Buffalo, New York 14260-3000, USA
| | - Alberto Baiardi
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Stefano Battaglia
- Department of Chemistry - BMC, Uppsala University, P.O. Box 576, SE-751 23 Uppsala, Sweden
| | - Veniamin A Borin
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | - Liviu F Chibotaru
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Irene Conti
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, Bologna I-40136, Italy
| | - Luca De Vico
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, via Aldo Moro 2, 53100 Siena, Italy
| | - Mickaël Delcey
- Department of Chemistry - Ångström Laboratory, Uppsala University, SE-751 21 Uppsala, Sweden
| | - Ignacio Fdez Galván
- Department of Chemistry - BMC, Uppsala University, P.O. Box 576, SE-751 23 Uppsala, Sweden
| | - Nicolas Ferré
- Aix-Marseille University, CNRS, Institut Chimie Radicalaire, Marseille, France
| | - Leon Freitag
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, Bologna I-40136, Italy
| | - Xuejun Gong
- Department of Chemistry, University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Stefan Knecht
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Ernst D Larsson
- Division of Theoretical Chemistry, Lund University, P.O. Box 124, Lund 22100, Sweden
| | - Roland Lindh
- Department of Chemistry - BMC, Uppsala University, P.O. Box 576, SE-751 23 Uppsala, Sweden
| | - Marcus Lundberg
- Department of Chemistry - Ångström Laboratory, Uppsala University, SE-751 21 Uppsala, Sweden
| | - Per Åke Malmqvist
- Division of Theoretical Chemistry, Lund University, P.O. Box 124, Lund 22100, Sweden
| | - Artur Nenov
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, Bologna I-40136, Italy
| | - Jesper Norell
- Department of Physics, AlbaNova University Center, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Michael Odelius
- Department of Physics, AlbaNova University Center, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Massimo Olivucci
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, via Aldo Moro 2, 53100 Siena, Italy
| | - Thomas B Pedersen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Laura Pedraza-González
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, via Aldo Moro 2, 53100 Siena, Italy
| | - Quan M Phung
- Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Chikusa, Nagoya 464-8602, Japan
| | - Kristine Pierloot
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Igor Schapiro
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | - Javier Segarra-Martí
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, 80 Wood Lane, London W12 0BZ, United Kingdom
| | - Francesco Segatta
- Dipartimento di Chimica Industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, Bologna I-40136, Italy
| | - Luis Seijo
- Departamento de Química, Instituto Universitario de Ciencia de Materiales Nicolás Cabrera, and Condensed Matter Physics Center (IFIMAC), Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Saumik Sen
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | | | - Christopher J Stein
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Liviu Ungur
- Department of Chemistry, University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Morgane Vacher
- Laboratoire CEISAM - UMR CNRS 6230, Université de Nantes, 44300 Nantes, France
| | - Alessio Valentini
- Theoretical Physical Chemistry, Research Unit MolSys, Université de Liège, Allée du 6 Août, 11, 4000 Liège, Belgium
| | - Valera Veryazov
- Division of Theoretical Chemistry, Lund University, P.O. Box 124, Lund 22100, Sweden
| |
Collapse
|
8
|
Pedraza-González L, Marín MDC, Jorge AN, Ruck TD, Yang X, Valentini A, Olivucci M, De Vico L. Web-ARM: A Web-Based Interface for the Automatic Construction of QM/MM Models of Rhodopsins. J Chem Inf Model 2020; 60:1481-1493. [PMID: 31909998 PMCID: PMC7101466 DOI: 10.1021/acs.jcim.9b00615] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
This article introduces Web-ARM, a specialized tool, online available, designed to build quantum mechanical/molecular mechanical models of rhodopsins, a widely spread family of light-responsive proteins. Web-ARM allows the rapidly building of models of rhodopsins with a documented quality and the prediction of trends in UV-vis absorption maximum wavelengths, based on their excitation energies computed at the CASPT2//CASSCF/Amber level of theory. Web-ARM builds upon the recently reported, python-based a-ARM protocol [J. Chem. Theory Comput., 2019, 15, 3134-3152] and, as such, necessitates only a crystallographic structure or a comparative model in PDB format and a very basic knowledge of the studied rhodopsin system. The user-friendly web interface uses such input to generate congruous, gas-phase models of rhodopsins and, if requested, their mutants. We present two possible applications of Web-ARM, which showcase how the interface can be employed to assist both research and educational activities in fields at the interface between chemistry and biology. The first application shows how, through Web-ARM, research projects (e.g., rhodopsin and rhodopsin mutant screening) can be carried out in significantly less time with respect to using the required computational photochemistry tools via a command line. The second application documents the use of Web-ARM in a real-life educational/training activity, through a hands-on experience illustrating the concepts of rhodopsin color tuning.
Collapse
Affiliation(s)
- Laura Pedraza-González
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Via A. Moro 2, I-53100 Siena, Italy
| | - María Del Carmen Marín
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Via A. Moro 2, I-53100 Siena, Italy
| | - Alejandro N Jorge
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States of America
| | - Tyler D Ruck
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States of America
| | - Xuchun Yang
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States of America
| | - Alessio Valentini
- Theoretical Physical Chemistry, Research Unit MolSys, Université de Liège, Allée du 6 Août, 11, 4000 Liège, Belgium
| | - Massimo Olivucci
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Via A. Moro 2, I-53100 Siena, Italy
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States of America
| | - Luca De Vico
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Via A. Moro 2, I-53100 Siena, Italy
| |
Collapse
|
9
|
Affiliation(s)
- André Anda
- Chemical and Quantum Physics, School of Science, RMIT University, Melbourne, VIC 3001, Australia
- ARC Centre of Excellence in Exciton Science, School of Science, RMIT University, Melbourne, VIC 3001, Australia
| | - Thorsten Hansen
- Department of Chemistry, Copenhagen University, Universitetsparken 5, DK-2100, Copenhagen Ø, Denmark
| | - Luca De Vico
- Department of Biotechnologies, Chemistry and Pharmacy, University of Siena, Via Aldo Moro 2, I-53100, Siena, Italy
| |
Collapse
|
10
|
Pedraza-González L, De Vico L, del Carmen Marín M, Fanelli F, Olivucci M. a-ARM: Automatic Rhodopsin Modeling with Chromophore Cavity Generation, Ionization State Selection, and External Counterion Placement. J Chem Theory Comput 2019; 15:3134-3152. [PMID: 30916955 PMCID: PMC7141608 DOI: 10.1021/acs.jctc.9b00061] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The Automatic Rhodopsin Modeling (ARM) protocol has recently been proposed as a tool for the fast and parallel generation of basic hybrid quantum mechanics/molecular mechanics (QM/MM) models of wild type and mutant rhodopsins. However, in its present version, input preparation requires a few hours long user's manipulation of the template protein structure, which also impairs the reproducibility of the generated models. This limitation, which makes model building semiautomatic rather than fully automatic, comprises four tasks: definition of the retinal chromophore cavity, assignment of protonation states of the ionizable residues, neutralization of the protein with external counterions, and finally congruous generation of single or multiple mutations. In this work, we show that the automation of the original ARM protocol can be extended to a level suitable for performing the above tasks without user's manipulation and with an input preparation time of minutes. The new protocol, called a-ARM, delivers fully reproducible (i.e., user independent) rhodopsin QM/MM models as well as an improved model quality. More specifically, we show that the trend in vertical excitation energies observed for a set of 25 wild type and 14 mutant rhodopsins is predicted by the new protocol better than when using the original. Such an agreement is reflected by an estimated (relative to the probed set) trend deviation of 0.7 ± 0.5 kcal mol-1 (0.03 ± 0.02 eV) and mean absolute error of 1.0 kcal mol-1 (0.04 eV).
Collapse
Affiliation(s)
- Laura Pedraza-González
- Department of Biotechnologies, Chemistry and Pharmacy, Università degli Studi di Siena, via A. Moro 2, I-53100 Siena, Italy
| | - Luca De Vico
- Department of Biotechnologies, Chemistry and Pharmacy, Università degli Studi di Siena, via A. Moro 2, I-53100 Siena, Italy
| | - María del Carmen Marín
- Department of Biotechnologies, Chemistry and Pharmacy, Università degli Studi di Siena, via A. Moro 2, I-53100 Siena, Italy
| | - Francesca Fanelli
- Department of Life Sciences, Center for Neuroscience and Neurotechnology, Università degli Studi di Modena e Reggio Emilia, I-41125 Modena, Italy
| | - Massimo Olivucci
- Department of Biotechnologies, Chemistry and Pharmacy, Università degli Studi di Siena, via A. Moro 2, I-53100 Siena, Italy
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
| |
Collapse
|
11
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
12
|
Marsili E, Farag MH, Yang X, De Vico L, Olivucci M. Two-State, Three-Mode Parametrization of the Force Field of a Retinal Chromophore Model. J Phys Chem A 2019; 123:1710-1719. [PMID: 30753077 DOI: 10.1021/acs.jpca.8b10010] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In recent years, the potential energy surfaces of the penta-2,4-dieniminium cation have been investigated using several electronic structure methods. The resulting pool of geometrical, electronic, and energy data provides a suitable basis for the construction of a topographically correct analytical model of the molecule force field and, therefore, for a better understanding of this class of molecules, which includes the chromophore of visual pigments. In the present contribution, we report the construction of such a model for regions of the force field that drive the photochemical and thermal isomerization of the central double bound of the cation. While previous models included only two modes, it is here shown that the proposed three-mode model and corresponding set of parameters are able to reproduce the complex topographical and electronic structure features seen in electronically correlated data obtained at the XMCQDPT2//CASSCF/6-31G* level of theory.
Collapse
Affiliation(s)
- Emanuele Marsili
- Dipartimento di Biotecnologie, Chimica e Farmacia , Università di Siena , via A. Moro 2 , I-53100 Siena , Italy
| | - Marwa H Farag
- Department of Chemistry , University of Southern California , Los Angeles , California 90089-0482 , United States
| | - Xuchun Yang
- Department of Chemistry , Bowling Green State University , Bowling Green , Ohio 43403 , United States and
| | - Luca De Vico
- Dipartimento di Biotecnologie, Chimica e Farmacia , Università di Siena , via A. Moro 2 , I-53100 Siena , Italy
| | - Massimo Olivucci
- Dipartimento di Biotecnologie, Chimica e Farmacia , Università di Siena , via A. Moro 2 , I-53100 Siena , Italy.,Department of Chemistry , Bowling Green State University , Bowling Green , Ohio 43403 , United States and.,Institut de Physique et Chimie des Matériaux de Strasbourg, UMR 7504 , Université de Strasbourg-CNRS , F-67034 Strasbourg , France
| |
Collapse
|
13
|
Anda A, De Vico L, Hansen T. Intermolecular Modes between LH2 Bacteriochlorophylls and Protein Residues: The Effect on the Excitation Energies. J Phys Chem B 2017; 121:5499-5508. [PMID: 28485594 DOI: 10.1021/acs.jpcb.7b02071] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Light-harvesting system 2 (LH2) executes the primary processes of photosynthesis in purple bacteria; photon absorption, and energy transportation to the reaction center. A detailed mechanistic insight into these operations is obscured by the complexity of the light-harvesting systems, particularly by the chromophore-environment interaction. In this work, we focus on the effects of the protein residues that are ligated to the bacteriochlorophylls (BChls) and construct potential energy surfaces of the ground and first optically excited state for the various BChl-residue systems where we in each case consider two degrees of freedom in the intermolecular region. We find that the excitation energies are only slightly affected by the considered modes. In addition, we see that axial ligands and hydrogen-bonded residues have opposite effects on both excitation energies and oscillator strengths by comparing to the isolated BChls. Our results indicate that only a small part of the chromophore-environment interaction can be associated with the intermolecular region between a BChl and an adjacent residue, but that it may be possible to selectively raise or lower the excitation energy at the axial and planar residue positions, respectively.
Collapse
Affiliation(s)
- André Anda
- Department of Chemistry, H. C. Ørsted Institute, University of Copenhagen , Universitetsparken 5, DK-2100 Copenhagen, Denmark
| | - Luca De Vico
- Department of Chemistry, H. C. Ørsted Institute, University of Copenhagen , Universitetsparken 5, DK-2100 Copenhagen, Denmark.,Department of Biotechnology, Chemistry and Pharmacy, University of Siena , via Aldo Moro 2, 53100 Siena, Italy
| | - Thorsten Hansen
- Department of Chemistry, H. C. Ørsted Institute, University of Copenhagen , Universitetsparken 5, DK-2100 Copenhagen, Denmark
| |
Collapse
|
14
|
Affiliation(s)
- Alessandro Pirrotta
- Nano-Science Center and Department of Chemistry, University of Copenhagen, 2100 Copenhagen Ø, Denmark
| | - Luca De Vico
- Department of Chemistry, University of Copenhagen, 2100 Copenhagen Ø, Denmark
| | - Gemma C. Solomon
- Nano-Science Center and Department of Chemistry, University of Copenhagen, 2100 Copenhagen Ø, Denmark
| | - Ignacio Franco
- Departments of Chemistry and Physics, University of Rochester, Rochester, New York 14627-0216, USA
| |
Collapse
|
15
|
Bogh S, Simmermacher M, Westberg M, Bregnhøj M, Rosenberg M, De Vico L, Veiga M, Laursen BW, Ogilby PR, Sauer SPA, Sørensen TJ. Azadioxatriangulenium and Diazaoxatriangulenium: Quantum Yields and Fundamental Photophysical Properties. ACS Omega 2017; 2:193-203. [PMID: 31457221 PMCID: PMC6641101 DOI: 10.1021/acsomega.6b00211] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/26/2016] [Accepted: 01/06/2017] [Indexed: 05/24/2023]
Abstract
Over the last decade, we have investigated and exploited the photophysical properties of triangulenium dyes. Azadioxatriangulenium (ADOTA) and diazaoxatriangulenium (DAOTA), in particular, have features that make them useful in various fluorescence-based technologies (e.g., bioimaging). Through our work with ADOTA and DAOTA, we became aware that the reported fluorescence quantum yields (ϕfl) for these dyes are lower than their actual values. We thus set out to further investigate the fundamental structure-property relationships in these unique conjugated cationic systems. The nonradiative processes in the systems were explored using transient absorption spectroscopy and time-resolved emission spectroscopy in combination with computational chemistry. The influence of molecular oxygen on the fluorescence properties was explored, and the singlet oxygen sensitization efficiencies of ADOTA and DAOTA were determined. We conclude that, for these dyes, the amount of nonradiative deactivation of the first excited singlet state (S1) of the azaoxa-triangulenium fluorophores is low, that the rate of such deactivation is slower than what is observed in common cationic dyes, that there are no observable radiative transitions occurring from the first excited triplet state (T1) of these dyes, and that the efficiency of sensitized singlet oxygen production is low (ϕΔ ≤ 10%). These photophysical results provide a solid base upon which technological applications of these fluorescent dyes can be built.
Collapse
Affiliation(s)
- Sidsel
A. Bogh
- Nano-Science
Center & Department of Chemistry, University
of Copenhagen, Universitetsparken
5, 2100 København
Ø, Denmark
| | - Mats Simmermacher
- Nano-Science
Center & Department of Chemistry, University
of Copenhagen, Universitetsparken
5, 2100 København
Ø, Denmark
| | - Michael Westberg
- Department
of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Mikkel Bregnhøj
- Department
of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Martin Rosenberg
- Nano-Science
Center & Department of Chemistry, University
of Copenhagen, Universitetsparken
5, 2100 København
Ø, Denmark
| | - Luca De Vico
- Nano-Science
Center & Department of Chemistry, University
of Copenhagen, Universitetsparken
5, 2100 København
Ø, Denmark
| | - Manoel Veiga
- PicoQuant
GmbH, Rudower Chaussee
29, 12489 Berlin, Germany
| | - Bo W. Laursen
- Nano-Science
Center & Department of Chemistry, University
of Copenhagen, Universitetsparken
5, 2100 København
Ø, Denmark
| | - Peter R. Ogilby
- Department
of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Stephan P. A. Sauer
- Nano-Science
Center & Department of Chemistry, University
of Copenhagen, Universitetsparken
5, 2100 København
Ø, Denmark
| | - Thomas Just Sørensen
- Nano-Science
Center & Department of Chemistry, University
of Copenhagen, Universitetsparken
5, 2100 København
Ø, Denmark
| |
Collapse
|
16
|
Governa P, Giachetti D, Biagi M, Manetti F, De Vico L. Hypothesis on Serenoa repens (Bartram) small extract inhibition of prostatic 5 α-reductase through an in silico approach on 5 β-reductase x-ray structure. PeerJ 2016; 4:e2698. [PMID: 27904805 PMCID: PMC5126621 DOI: 10.7717/peerj.2698] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2016] [Accepted: 10/18/2016] [Indexed: 11/20/2022] Open
Abstract
Benign prostatic hyperplasia is a common disease in men aged over 50 years old, with an incidence increasing to more than 80% over the age of 70, that is increasingly going to attract pharmaceutical interest. Within conventional therapies, such as α-adrenoreceptor antagonists and 5α-reductase inhibitor, there is a large requirement for treatments with less adverse events on, e.g., blood pressure and sexual function: phytotherapy may be the right way to fill this need. Serenoa repens standardized extract has been widely studied and its ability to reduce lower urinary tract symptoms related to benign prostatic hyperplasia is comprehensively described in literature. An innovative investigation on the mechanism of inhibition of 5α-reductase by Serenoa repens extract active principles is proposed in this work through computational methods, performing molecular docking simulations on the crystal structure of human liver 5β-reductase. The results confirm that both sterols and fatty acids can play a role in the inhibition of the enzyme, thus, suggesting a competitive mechanism of inhibition. This work proposes a further confirmation for the rational use of herbal products in the management of benign prostatic hyperplasia, and suggests computational methods as an innovative, low cost, and non-invasive process for the study of phytocomplex activity toward proteic targets.
Collapse
Affiliation(s)
- Paolo Governa
- Department of Physical Sciences, Earth and Environment, University of Siena, Siena, Italy; Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Daniela Giachetti
- Department of Physical Sciences, Earth and Environment, University of Siena , Siena , Italy
| | - Marco Biagi
- Department of Physical Sciences, Earth and Environment, University of Siena , Siena , Italy
| | - Fabrizio Manetti
- Department of Biotechnology, Chemistry and Pharmacy, University of Siena , Siena , Italy
| | - Luca De Vico
- Department of Chemistry, University of Copenhagen , Copenhagen , Denmark
| |
Collapse
|
17
|
Affiliation(s)
- André Anda
- Department
of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100, København Ø, Denmark
| | - Luca De Vico
- Department
of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100, København Ø, Denmark
| | - Thorsten Hansen
- Department
of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100, København Ø, Denmark
| | - Darius Abramavičius
- Department
of Theoretical Physics, Vilnius University, Sauletekio ave. 9, Build. 3, LT-10222 Vilnius, Lithuania
| |
Collapse
|
18
|
Storm FE, Olsen ST, Hansen T, De Vico L, Jackson NE, Ratner MA, Mikkelsen KV. Boron Subphthalocyanine Based Molecular Triad Systems for the Capture of Solar Energy. J Phys Chem A 2016; 120:7694-7703. [DOI: 10.1021/acs.jpca.6b05518] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Freja E. Storm
- Department
of Chemistry, H. C. Ørsted Institute, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen, Denmark
| | - Stine T. Olsen
- Department
of Chemistry, H. C. Ørsted Institute, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen, Denmark
| | - Thorsten Hansen
- Department
of Chemistry, H. C. Ørsted Institute, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen, Denmark
| | - Luca De Vico
- Department
of Chemistry, H. C. Ørsted Institute, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen, Denmark
| | - Nicholas E. Jackson
- Department
of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Mark A. Ratner
- Department
of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Kurt V. Mikkelsen
- Department
of Chemistry, H. C. Ørsted Institute, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen, Denmark
| |
Collapse
|
19
|
Anda A, Hansen T, De Vico L. Multireference Excitation Energies for Bacteriochlorophylls A within Light Harvesting System 2. J Chem Theory Comput 2016; 12:1305-13. [DOI: 10.1021/acs.jctc.5b01104] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- André Anda
- Department of Chemistry, Copenhagen University, Universitetsparken 5, DK-2100, Copenhagen Ø, Denmark
| | - Thorsten Hansen
- Department of Chemistry, Copenhagen University, Universitetsparken 5, DK-2100, Copenhagen Ø, Denmark
| | - Luca De Vico
- Department of Chemistry, Copenhagen University, Universitetsparken 5, DK-2100, Copenhagen Ø, Denmark
| |
Collapse
|
20
|
Aquilante F, Autschbach J, Carlson RK, Chibotaru LF, Delcey MG, De Vico L, Fdez Galván I, Ferré N, Frutos LM, Gagliardi L, Garavelli M, Giussani A, Hoyer CE, Li Manni G, Lischka H, Ma D, Malmqvist PÅ, Müller T, Nenov A, Olivucci M, Pedersen TB, Peng D, Plasser F, Pritchard B, Reiher M, Rivalta I, Schapiro I, Segarra-Martí J, Stenrup M, Truhlar DG, Ungur L, Valentini A, Vancoillie S, Veryazov V, Vysotskiy VP, Weingart O, Zapata F, Lindh R. Molcas 8: New capabilities for multiconfigurational quantum chemical calculations across the periodic table. J Comput Chem 2015; 37:506-41. [PMID: 26561362 DOI: 10.1002/jcc.24221] [Citation(s) in RCA: 1083] [Impact Index Per Article: 120.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2015] [Revised: 09/07/2015] [Accepted: 09/09/2015] [Indexed: 12/17/2022]
Abstract
In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas-Kroll-Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC-PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large-scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.
Collapse
Affiliation(s)
- Francesco Aquilante
- Department of Chemistry - Ångström, The Theoretical Chemistry Programme, Uppsala University, Box 518, Uppsala, 751 20, Sweden.,Dipartimento di Chimica "G. Ciamician", Università di Bologna, via Selmi 2, IT-40126, Bologna, Italy
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York, 14260-3000, USA
| | - Rebecca K Carlson
- Department of Chemistry, Supercomputing Institute, and Chemical Theory Center, University of Minnesota, Minneapolis, Minnesota, 55455-0431, USA
| | - Liviu F Chibotaru
- Division of Quantum and Physical Chemistry, and INPAC, Institute for Nanoscale Physics and Chemistry, Katholieke Universiteit Leuven Celestijnenlaan, 200F, 3001, Belgium
| | - Mickaël G Delcey
- Department of Chemistry - Ångström, The Theoretical Chemistry Programme, Uppsala University, Box 518, Uppsala, 751 20, Sweden
| | - Luca De Vico
- Department of Chemistry, Copenhagen University, Universitetsparken 5, Copenhagen Ø, 2100, Denmark
| | - Ignacio Fdez Galván
- Department of Chemistry - Ångström, The Theoretical Chemistry Programme, Uppsala University, Box 518, Uppsala, 751 20, Sweden.,Uppsala Center for Computational Chemistry - UC3, Uppsala University, Box 518, Uppsala, 751 20, Sweden
| | - Nicolas Ferré
- Université d'Aix-Marseille, CNRS, Institut de Chimie Radicalaire, Campus Étoile/Saint-Jérôme Case 521, Avenue Esc. Normandie Niemen, Marseille Cedex 20, 13397, France
| | - Luis Manuel Frutos
- Unidad Docente de Química Física, Universidad de Alcalá, E-28871 Alcalá de Henares, Madrid, Spain
| | - Laura Gagliardi
- Department of Chemistry, Supercomputing Institute, and Chemical Theory Center, University of Minnesota, Minneapolis, Minnesota, 55455-0431, USA
| | - Marco Garavelli
- Dipartimento di Chimica "G. Ciamician", Università di Bologna, via Selmi 2, IT-40126, Bologna, Italy.,Université de Lyon, CNRS, École Normale Supérieure de Lyon, 46 Allée d'Italie, Lyon Cedex 07, F-69364, France
| | - Angelo Giussani
- Dipartimento di Chimica "G. Ciamician", Università di Bologna, via Selmi 2, IT-40126, Bologna, Italy
| | - Chad E Hoyer
- Department of Chemistry, Supercomputing Institute, and Chemical Theory Center, University of Minnesota, Minneapolis, Minnesota, 55455-0431, USA
| | - Giovanni Li Manni
- Department of Chemistry, Supercomputing Institute, and Chemical Theory Center, University of Minnesota, Minneapolis, Minnesota, 55455-0431, USA.,Max Planck Institut für Festkörperforschung, Heisenbergstraße 1, Stuttgart, 70569, Germany
| | - Hans Lischka
- Department of Chemistry and Biochemistry, Texas Tech University, Memorial Circle and Boston, Lubbock, Texas, 79409-1061, USA.,Institute for Theoretical Chemistry, University of Vienna, Währingerstraße 17, Vienna, A-1090, Austria
| | - Dongxia Ma
- Department of Chemistry, Supercomputing Institute, and Chemical Theory Center, University of Minnesota, Minneapolis, Minnesota, 55455-0431, USA.,Max Planck Institut für Festkörperforschung, Heisenbergstraße 1, Stuttgart, 70569, Germany
| | - Per Åke Malmqvist
- Department of Theoretical Chemistry, Lund University, Chemical Center, P.O.B 124 S-221 00, Lund, Sweden
| | - Thomas Müller
- Jülich Supercomputing Centre (JSC), Forschungszentrum Jülich GmbH, Institute for Advanced Simulation (IAS), Wilhelm-Johnen-Straße, Jülich, 52425, Germany
| | - Artur Nenov
- Dipartimento di Chimica "G. Ciamician", Università di Bologna, via Selmi 2, IT-40126, Bologna, Italy
| | - Massimo Olivucci
- Department of Biotechnology, Chemistry and Pharmacy, University of Siena, via Aldo Moro 2, Siena, 53100, Italy.,Chemistry Department, Bowling Green State University, 141 Overman Hall, Bowling Green, Ohio, 43403, USA.,Institut de Physique et Chimie des Matériaux de Strasbourg & Labex NIE, Université de Strasbourg, CNRS UMR 7504, 23 Rue du Loess, Strasbourg, 67034, France
| | - Thomas Bondo Pedersen
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, Oslo, 0315, Norway
| | - Daoling Peng
- College of Chemistry and Environment, South China Normal University, Guangzhou, 510006, China
| | - Felix Plasser
- Institute for Theoretical Chemistry, University of Vienna, Währingerstraße 17, Vienna, A-1090, Austria
| | - Ben Pritchard
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York, 14260-3000, USA
| | - Markus Reiher
- ETH Zurich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, Zurich, CH-8093, Switzerland
| | - Ivan Rivalta
- Université de Lyon, CNRS, École Normale Supérieure de Lyon, 46 Allée d'Italie, Lyon Cedex 07, F-69364, France
| | - Igor Schapiro
- Institut de Physique et Chimie des Matériaux de Strasbourg & Labex NIE, Université de Strasbourg, CNRS UMR 7504, 23 Rue du Loess, Strasbourg, 67034, France.,Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, 91904, Israel
| | - Javier Segarra-Martí
- Dipartimento di Chimica "G. Ciamician", Università di Bologna, via Selmi 2, IT-40126, Bologna, Italy
| | - Michael Stenrup
- Department of Chemistry - Ångström, The Theoretical Chemistry Programme, Uppsala University, Box 518, Uppsala, 751 20, Sweden.,Uppsala Center for Computational Chemistry - UC3, Uppsala University, Box 518, Uppsala, 751 20, Sweden
| | - Donald G Truhlar
- Department of Chemistry, Supercomputing Institute, and Chemical Theory Center, University of Minnesota, Minneapolis, Minnesota, 55455-0431, USA
| | - Liviu Ungur
- Division of Quantum and Physical Chemistry, and INPAC, Institute for Nanoscale Physics and Chemistry, Katholieke Universiteit Leuven Celestijnenlaan, 200F, 3001, Belgium
| | - Alessio Valentini
- Unidad Docente de Química Física, Universidad de Alcalá, E-28871 Alcalá de Henares, Madrid, Spain.,Department of Biotechnology, Chemistry and Pharmacy, University of Siena, via Aldo Moro 2, Siena, 53100, Italy
| | - Steven Vancoillie
- Department of Theoretical Chemistry, Lund University, Chemical Center, P.O.B 124 S-221 00, Lund, Sweden
| | - Valera Veryazov
- Department of Theoretical Chemistry, Lund University, Chemical Center, P.O.B 124 S-221 00, Lund, Sweden
| | - Victor P Vysotskiy
- Department of Theoretical Chemistry, Lund University, Chemical Center, P.O.B 124 S-221 00, Lund, Sweden
| | - Oliver Weingart
- Institut für Theoretische Chemie und Computerchemie, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, Düsseldorf, 40225, Germany
| | - Felipe Zapata
- Unidad Docente de Química Física, Universidad de Alcalá, E-28871 Alcalá de Henares, Madrid, Spain
| | - Roland Lindh
- Department of Chemistry - Ångström, The Theoretical Chemistry Programme, Uppsala University, Box 518, Uppsala, 751 20, Sweden.,Uppsala Center for Computational Chemistry - UC3, Uppsala University, Box 518, Uppsala, 751 20, Sweden
| |
Collapse
|
21
|
Jensen JH, Willemoës M, Winther JR, De Vico L. In silico prediction of mutant HIV-1 proteases cleaving a target sequence. PLoS One 2014; 9:e95833. [PMID: 24796579 PMCID: PMC4010418 DOI: 10.1371/journal.pone.0095833] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2013] [Accepted: 03/31/2014] [Indexed: 11/17/2022] Open
Abstract
HIV-1 protease represents an appealing system for directed enzyme re-design, since it has various different endogenous targets, a relatively simple structure and it is well studied. Recently Chaudhury and Gray (Structure (2009) 17: 1636–1648) published a computational algorithm to discern the specificity determining residues of HIV-1 protease. In this paper we present two computational tools aimed at re-designing HIV-1 protease, derived from the algorithm of Chaudhuri and Gray. First, we present an energy-only based methodology to discriminate cleavable and non cleavable peptides for HIV-1 proteases, both wild type and mutant. Secondly, we show an algorithm we developed to predict mutant HIV-1 proteases capable of cleaving a new target substrate peptide, different from the natural targets of HIV-1 protease. The obtained in silico mutant enzymes were analyzed in terms of cleavability and specificity towards the target peptide using the energy-only methodology. We found two mutant proteases as best candidates for specificity and cleavability towards the target sequence.
Collapse
Affiliation(s)
- Jan H Jensen
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Martin Willemoës
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Jakob R Winther
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Luca De Vico
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
22
|
Hediger MR, De Vico L, Rannes JB, Jäckel C, Besenmatter W, Svendsen A, Jensen JH. In silico screening of 393 mutants facilitates enzyme engineering of amidase activity in CalB. PeerJ 2013; 1:e145. [PMID: 24010022 PMCID: PMC3757469 DOI: 10.7717/peerj.145] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2013] [Accepted: 08/06/2013] [Indexed: 11/25/2022] Open
Abstract
Our previously presented method for high throughput computational screening of mutant activity (Hediger et al., 2012) is benchmarked against experimentally measured amidase activity for 22 mutants of Candida antarctica lipase B (CalB). Using an appropriate cutoff criterion for the computed barriers, the qualitative activity of 15 out of 22 mutants is correctly predicted. The method identifies four of the six most active mutants with ≥3-fold wild type activity and seven out of the eight least active mutants with ≤0.5-fold wild type activity. The method is further used to screen all sterically possible (386) double-, triple- and quadruple-mutants constructed from the most active single mutants. Based on the benchmark test at least 20 new promising mutants are identified.
Collapse
Affiliation(s)
- Martin R Hediger
- Department of Chemistry, University of Copenhagen , Copenhagen , Denmark
| | | | | | | | | | | | | |
Collapse
|
23
|
Hediger MR, Steinmann C, De Vico L, Jensen JH. A computational method for the systematic screening of reaction barriers in enzymes: searching for Bacillus circulans xylanase mutants with greater activity towards a synthetic substrate. PeerJ 2013; 1:e111. [PMID: 23904990 PMCID: PMC3728886 DOI: 10.7717/peerj.111] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2013] [Accepted: 07/03/2013] [Indexed: 12/04/2022] Open
Abstract
We present a semi-empirical (PM6-based) computational method for systematically estimating the effect of all possible single mutants, within a certain radius of the active site, on the barrier height of an enzymatic reaction. The intent of this method is not a quantitative prediction of the barrier heights, but rather to identify promising mutants for further computational or experimental study. The method is applied to identify promising single and double mutants of Bacillus circulans xylanase (BCX) with increased hydrolytic activity for the artificial substrate ortho-nitrophenyl β-xylobioside (ONPX2). The estimated reaction barrier for wild-type (WT) BCX is 18.5 kcal/mol, which is in good agreement with the experimental activation free energy value of 17.0 kcal/mol extracted from the observed k cat using transition state theory (Joshi et al., 2001). The PM6 reaction profiles for eight single point mutations are recomputed using FMO-MP2/PCM/6-31G(d) single points. PM6 predicts an increase in barrier height for all eight mutants while FMO predicts an increase for six of the eight mutants. Both methods predict that the largest change in barrier occurs for N35F, where PM6 and FMO predict a 9.0 and 15.8 kcal/mol increase, respectively. We thus conclude that PM6 is sufficiently accurate to identify promising mutants for further study. We prepared a set of all theoretically possible (342) single mutants in which every amino acid of the active site (except for the catalytically active residues E78 and E172) was mutated to every other amino acid. Based on results from the single mutants we construct a set of 111 double mutants consisting of all possible pairs of single mutants with the lowest barrier for a particular position and compute their reaction profile. None of the mutants have, to our knowledge, been prepared experimentally and therefore present experimentally testable predictions.
Collapse
Affiliation(s)
- Martin R. Hediger
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Casper Steinmann
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Odense, Denmark
| | - Luca De Vico
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Jan H. Jensen
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
24
|
Lloret N, Frederiksen RS, Møller TC, Rieben NI, Upadhyay S, De Vico L, Jensen JH, Nygård J, Martinez KL. Effects of buffer composition and dilution on nanowire field-effect biosensors. Nanotechnology 2013; 24:035501. [PMID: 23263553 DOI: 10.1088/0957-4484/24/3/035501] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
Nanowire-based field-effect transistors (FETs) can be used as ultra-sensitive and label-free biosensors for detecting protein-protein interactions. A way to increase the performance of such sensors is to dilute the sensing buffer drastically. However, we show here that this can have an important effect on the function of the proteins. Moreover, it is demonstrated that this dilution significantly affects the pH stability of the sensing buffer, which consequently impacts the charge of the protein and thus the response and signal-to-noise ratio in the sensing experiments. Three model systems are investigated experimentally to illustrate the impact on ligand-protein and protein-protein interactions. Simulations are performed to illustrate the effect on the performance of the sensors. Combining various parameters, the current study provides a means for evaluating and selecting the most appropriate buffer composition for bioFET measurements.
Collapse
Affiliation(s)
- Noémie Lloret
- Bio-Nanotechnology and Nanomedicine Laboratory, Department of Chemistry & Nano-Science Center, University of Copenhagen, Copenhagen, Denmark
| | | | | | | | | | | | | | | | | |
Collapse
|
25
|
Hediger MR, De Vico L, Svendsen A, Besenmatter W, Jensen JH. A computational methodology to screen activities of enzyme variants. PLoS One 2012; 7:e49849. [PMID: 23284627 PMCID: PMC3524253 DOI: 10.1371/journal.pone.0049849] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2012] [Accepted: 10/14/2012] [Indexed: 11/22/2022] Open
Abstract
We present a fast computational method to efficiently screen enzyme activity. In the presented method, the effect of mutations on the barrier height of an enzyme-catalysed reaction can be computed within 24 hours on roughly 10 processors. The methodology is based on the PM6 and MOZYME methods as implemented in MOPAC2009, and is tested on the first step of the amide hydrolysis reaction catalyzed by the Candida Antarctica lipase B (CalB) enzyme. The barrier heights are estimated using adiabatic mapping and shown to give barrier heights to within 3 kcal/mol of B3LYP/6-31G(d)//RHF/3-21G results for a small model system. Relatively strict convergence criteria (0.5 kcal/(molÅ)), long NDDO cutoff distances within the MOZYME method (15 Å) and single point evaluations using conventional PM6 are needed for reliable results. The generation of mutant structures and subsequent setup of the semiempirical calculations are automated so that the effect on barrier heights can be estimated for hundreds of mutants in a matter of weeks using high performance computing.
Collapse
Affiliation(s)
- Martin R. Hediger
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Luca De Vico
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | | | | | - Jan H. Jensen
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
- * E-mail:
| |
Collapse
|
26
|
De Vico L, Iversen L, Sørensen MH, Brandbyge M, Nygård J, Martinez KL, Jensen JH. Predicting and rationalizing the effect of surface charge distribution and orientation on nano-wire based FET bio-sensors. Nanoscale 2011; 3:3635-3640. [PMID: 21811738 DOI: 10.1039/c1nr10316d] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
A single charge screening model of surface charge sensors in liquids (De Vico et al., Nanoscale, 2011, 3, 706-717) is extended to multiple charges to model the effect of the charge distributions of analyte proteins on FET sensor response. With this model we show that counter-intuitive signal changes (e.g. a positive signal change due to a net positive protein binding to a p-type conductor) can occur for certain combinations of charge distributions and Debye lengths. The new method is applied to interpret published experimental data on Streptavidin (Ishikawa et al., ACS Nano, 2009, 3, 3969-3976) and Nucleocapsid protein (Ishikawa et al., ACS Nano, 2009, 3, 1219-1224).
Collapse
Affiliation(s)
- Luca De Vico
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100, Copenhagen, Denmark.
| | | | | | | | | | | | | |
Collapse
|
27
|
De Vico L, Sørensen MH, Iversen L, Rogers DM, Sørensen BS, Brandbyge M, Nygård J, Martinez KL, Jensen JH. Quantifying signal changes in nano-wire based biosensors. Nanoscale 2011; 3:706-717. [PMID: 21173975 DOI: 10.1039/c0nr00442a] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
In this work, we present a computational methodology for predicting the change in signal (conductance sensitivity) of a nano-BioFET sensor (a sensor based on a biomolecule binding another biomolecule attached to a nano-wire field effect transistor) upon binding its target molecule. The methodology is a combination of the screening model of surface charge sensors in liquids developed by Brandbyge and co-workers [Sørensen et al., Appl. Phys. Lett., 2007, 91, 102105], with the PROPKA method for predicting the pH-dependent charge of proteins and protein-ligand complexes, developed by Jensen and co-workers [Li et al., Proteins: Struct., Funct., Bioinf., 2005, 61, 704-721, Bas et al., Proteins: Struct., Funct., Bioinf., 2008, 73, 765-783]. The predicted change in conductance sensitivity based on this methodology is compared to previously published data on nano-BioFET sensors obtained by other groups. In addition, the conductance sensitivity dependence from various parameters is explored for a standard wire, representative of a typical experimental setup. In general, the experimental data can be reproduced with sufficient accuracy to help interpret them. The method has the potential for even more quantitative predictions when key experimental parameters (such as the charge carrier density of the nano-wire or receptor density on the device surface) can be determined (and reported) more accurately.
Collapse
Affiliation(s)
- Luca De Vico
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100, Denmark
| | | | | | | | | | | | | | | | | |
Collapse
|
28
|
Aquilante F, De Vico L, Ferré N, Ghigo G, Malmqvist PA, Neogrády P, Pedersen TB, Pitonák M, Reiher M, Roos BO, Serrano-Andrés L, Urban M, Veryazov V, Lindh R. MOLCAS 7: the next generation. J Comput Chem 2010; 31:224-47. [PMID: 19499541 DOI: 10.1002/jcc.21318] [Citation(s) in RCA: 1301] [Impact Index Per Article: 92.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Some of the new unique features of the MOLCAS quantum chemistry package version 7 are presented in this report. In particular, the Cholesky decomposition method applied to some quantum chemical methods is described. This approach is used both in the context of a straight forward approximation of the two-electron integrals and in the generation of so-called auxiliary basis sets. The article describes how the method is implemented for most known wave functions models: self-consistent field, density functional theory, 2nd order perturbation theory, complete-active space self-consistent field multiconfigurational reference 2nd order perturbation theory, and coupled-cluster methods. The report further elaborates on the implementation of a restricted-active space self-consistent field reference function in conjunction with 2nd order perturbation theory. The average atomic natural orbital basis for relativistic calculations, covering the whole periodic table, are described and associated unique properties are demonstrated. Furthermore, the use of the arbitrary order Douglas-Kroll-Hess transformation for one-component relativistic calculations and its implementation are discussed. This section especially focuses on the implementation of the so-called picture-change-free atomic orbital property integrals. Moreover, the ElectroStatic Potential Fitted scheme, a version of a quantum mechanics/molecular mechanics hybrid method implemented in MOLCAS, is described and discussed. Finally, the report discusses the use of the MOLCAS package for advanced studies of photo chemical phenomena and the usefulness of the algorithms for constrained geometry optimization in MOLCAS in association with such studies.
Collapse
Affiliation(s)
- Francesco Aquilante
- Department of Physical Chemistry, University of Geneva, 30 Quai Ernest Ansermet, CH-1211 Geneva, Switzerland
| | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
29
|
Liu F, Liu Y, Vico LD, Lindh R. A CASSCF/CASPT2 approach to the decomposition of thiazole-substituted dioxetanone: Substitution effects and charge-transfer induced electron excitation. Chem Phys Lett 2009. [DOI: 10.1016/j.cplett.2009.11.009] [Citation(s) in RCA: 42] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
30
|
Affiliation(s)
- Fengyi Liu
- Department of Theoretical Chemistry, Lund University, Chemical Center, P.O. Box 124, SE-221 00 Lund, Sweden, and College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Yajun Liu
- Department of Theoretical Chemistry, Lund University, Chemical Center, P.O. Box 124, SE-221 00 Lund, Sweden, and College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Luca De Vico
- Department of Theoretical Chemistry, Lund University, Chemical Center, P.O. Box 124, SE-221 00 Lund, Sweden, and College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Roland Lindh
- Department of Theoretical Chemistry, Lund University, Chemical Center, P.O. Box 124, SE-221 00 Lund, Sweden, and College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| |
Collapse
|
31
|
Affiliation(s)
- Luca De Vico
- Department of Theoretical Chemistry, Lund University, P.O. Box 124, SE-22100 Lund, Sweden
| | - Roland Lindh
- Department of Theoretical Chemistry, Lund University, P.O. Box 124, SE-22100 Lund, Sweden
| |
Collapse
|
32
|
De Vico L, Pegado L, Heimdal J, Söderhjelm P, Roos BO. The ozone ring closure as a test for multi-state multi-configurational second order perturbation theory (MS-CASPT2). Chem Phys Lett 2008. [DOI: 10.1016/j.cplett.2008.06.065] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
33
|
Abstract
The thermal decomposition of 1,2-dioxetane and the associated production of chemiluminescent products, model for a wide range of chemiluminescent reactions, has been studied at the multistate multiconfigurational second-order perturbation level of theory. This study is in qualitative and quantitative agreement with experimental observations with respect to the activation energy and the observed increase of triplet and singlet excited products as substituents are added to the parent molecule. The, previously incomplete, reaction mechanism of the chemiluminescence of 1,2-dioxetane is now rationalized and described as mainly due to a particular form of entropic trapping.
Collapse
Affiliation(s)
- Luca De Vico
- Department of Theoretical Chemistry and Chemical Physics, Lund University, Chemical Centre, P. O. Box 188, S-22100 Lund, Sweden.
| | | | | | | |
Collapse
|
34
|
Abstract
The UV photodissociation (<5 eV) of diiodomethane (CH(2)I(2)) is investigated by spin-orbit ab initio calculations. The experimentally observed photodissociation channels in the gas and condensed phases are clearly assigned by multi-state second-order multiconfigurational perturbation theory in conjunction with spin-orbit interaction through complete active space-state interaction potential energy curves. The calculated results indicate that the fast dissociations of the first two singlet states of CH(2)I(2) and CH(2)I--I lead to geminate-radical products, CH(2)I (.)+I((2)P(3/2)) or CH(2)I (.)+ I*((2)P(1/2)). The recombination process from CH(2)I--I to CH(2)I(2) is explained by an isomerization process and a secondary photodissociation reaction of CH(2)I--I. Finally, the study reveals that spin-orbits effects are significant in the quantitative analysis of the electronic spectrum of the CH(2)I--I species.
Collapse
Affiliation(s)
- Ya-Jun Liu
- College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China.
| | | | | | | |
Collapse
|
35
|
|
36
|
Affiliation(s)
- Luca De Vico
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, S-221 00 Lund, Sweden
| | - Massimo Olivucci
- Dipartimento di Chimica, Università di Siena, via A. De Gasperi 2, I-53100 Siena, Italy, and Centro per lo Studio dei Sistemi Complessi, Via Tommaso Pendola 37, I-53100 Siena, Italy
| | - Roland Lindh
- Department of Chemical Physics, Lund University, Chemical Centre, P.O. Box 124, S-221 00 Lund, Sweden
| |
Collapse
|
37
|
De Vico L, Garavelli M, Bernardi F, Olivucci M. Photoisomerization Mechanism of 11-cis-Locked Artificial Retinal Chromophores: Acceleration and Primary Photoproduct Assignment. J Am Chem Soc 2005; 127:2433-42. [PMID: 15724998 DOI: 10.1021/ja045747u] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
CASPT2//CASSCF/6-31G photochemical reaction path computations for two 4-cis-nona-2,4,6,8-tetraeniminium cation derivatives, with the 4-cis double bond embedded in a seven- and eight-member ring, are carried out to model the reactivity of the corresponding ring-locked retinal chromophores. The comparison of the excited state branches of the two reaction paths with that of the native chromophore, is used to unveil the factors responsible for the remarkably short (60 fs) excited state (S(1)) lifetime observed when an artificial rhodopsin containing an eight member ring-locked retinal is photoexcited. Indeed, it is shown that the strain imposed by the eight-member ring on the chromophore backbone leads to a dramatic change in the shape of the S(1) energy surface. Our models are also used to investigate the nature of the primary photoproducts observed in different artificial rhodopsins. It is seen that only the eight member ring-locked retinal model can access a shallow energy minimum on the ground state. This result implies that the primary, photorhodopsin-like, transient observed in artificial rhodopsins could correspond to a shallow excited state minimum. Similarly, the second, bathorhodopsin-like, transient species could be assigned to a ground state structure displaying a nearly all-trans conformation.
Collapse
Affiliation(s)
- Luca De Vico
- Dipartimento di Chimica, Università di Siena, via De Gasperi 2 Siena, I-53100 Italy
| | | | | | | |
Collapse
|
38
|
Helbing J, Bregy H, Bredenbeck J, Pfister R, Hamm P, Huber R, Wachtveitl J, De Vico L, Olivucci M. A Fast Photoswitch for Minimally Perturbed Peptides: Investigation of the trans → cis Photoisomerization of N-Methylthioacetamide. J Am Chem Soc 2004; 126:8823-34. [PMID: 15250736 DOI: 10.1021/ja049227a] [Citation(s) in RCA: 74] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Thio amino acids can be integrated into the backbone of peptides without significantly perturbing their structure. In this contribution we use ultrafast infrared and visible spectroscopy as well as state-of-the-art ab initio computations to investigate the photoisomerization of the trans form of N-methylthioacetamide (NMTAA) as a model conformational photoswitch. Following the S2 excitation of trans-NMTAA in water, the return of the molecule into the trans ground state and the formation of the cis isomer is observed on a dual time scale, with a fast component of 8-9 ps and a slow time constant of approximately 250 ps. On both time scales the probability of isomerization to the cis form is found to be 30-40%, independently of excitation wavelength. Ab initio CASPT2//CASSCF photochemical reaction path calculations indicate that, in vacuo, the trans-->cis isomerization event takes place on the S1 and/or T1 triplet potential energy surfaces and is controlled by very small energy barriers, in agreement with the experimentally observed picosecond time scale. Furthermore, the calculations identify one S2/S1 and four nearly isoenergetic S1/S0 conical intersection decay channels. In line with the observed isomerization probability, only one of the S1/S0 conical intersections yields the cis conformation upon S1-->S0 decay. A substantially equivalent excited-state relaxation results from four T1/S0 intersystem crossing points.
Collapse
Affiliation(s)
- Jan Helbing
- Physikalisch-Chemisches Institut, Winterthurerstrasse 190, Universität Zürich, CH-8057 Zürich, Switzerland
| | | | | | | | | | | | | | | | | |
Collapse
|
39
|
Abstract
The results of new and recently reported CASSCF/6-31G* photoisomerization path computations of a series of models of the 11-cis retinal chromophore of the visual pigment rhodopsin are discussed. The results indicate that, with respect to the chromophore in vacuo, certain structural, intramolecular and environmental factors are capable of speeding up the excited-state decay associated with the cis --> trans isomerization motion. Using suitable protonated Schiff-base models, it is shown that three structural factors can potentially speed up the isomerization: (i) reducing the length of the conjugated chain, (ii) twisting of the hydrocarbon end of the conjugated chain with respect to the protonated Schiff-base end and (iii) ring locking of the conjugated chain with an eight-membered ring. All these factors operate through increasing the slope of the excited-state energy surface and enhancing the coupling between stretching and torsional modes. We argue that the protein catalysis seen in rhodopsin may, at least partly, exploit the same principles.
Collapse
Affiliation(s)
- Adalgisa Sinicropi
- Dipartimento di Chimica, Università di Siena, via Aldo Moro, I-53100 Siena, Italy
| | | | | | | |
Collapse
|
40
|
De Vico L, Page CS, Garavelli M, Bernardi F, Basosi R, Olivucci M. Reaction path analysis of the "tunable" photoisomerization selectivity of free and locked retinal chromophores. J Am Chem Soc 2002; 124:4124-34. [PMID: 11942852 DOI: 10.1021/ja017502c] [Citation(s) in RCA: 63] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Multiconfigurational second-order perturbation theory computations and reaction path mapping for the retinal protonated Schiff base models all-trans-nona-2,4,6,8-tetraeniminium and 2-cis-nona-2,4,6,8-tetraeniminium cation demonstrate that, in isolated conditions, retinal chromophores exhibit at least three competing excited-state double bond isomerization paths. These paths are associated with the photoisomerization of the double bonds in positions 9, 11, and 13, respectively, and are controlled by barriers that favor the position 11. The computations provide a basis for the understanding of the observed excited-state lifetime in both naturally occurring and synthetic chromophores in solution and, tentatively, in the protein environment. In particular, we provide a rationalization of the excited-state lifetimes observed for a group of locked retinal chromophores which suggests that photoisomerization in bacteriorhodopsin is the result of simultaneous specific "catalysis" (all-trans --> 13-cis path) accompanied by specific "inhibition" (all-trans --> 11-cis path). The nature of the S(1) --> S(0) decay channel associated with the three paths has also been investigated at the CASSCF level of theory. It is shown that the energy surfaces in the vicinity of the conical intersection for the photoisomerization about the central double bond of retinal (position 11) and the two corresponding lateral double bonds (positions 9 and 13) are structurally different.
Collapse
Affiliation(s)
- Luca De Vico
- Dipartimento di Chimica, Università di Siena, via Aldo Moro, Siena, I-53100 Italy
| | | | | | | | | | | |
Collapse
|