1
|
Palmisano VF, Faraji S, Nogueira JJ. Effect of stacking interactions on charge transfer states in photoswitches interacting with ion channels. Phys Chem Chem Phys 2023; 25:8331-8335. [PMID: 36883995 PMCID: PMC10696555 DOI: 10.1039/d2cp05678j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2022] [Accepted: 02/24/2023] [Indexed: 03/06/2023]
Abstract
The activity of ion channels can be reversibly photo-controlled via the binding of molecular photoswitches, often based on an azobenzene scaffold. Those azobenzene derivatives interact with aromatic residues of the protein via stacking interactions. In the present work, the effect of face-to-face and t-shaped stacking interactions on the excited state electronic structure of azobenzene and p-diaminoazobenzene integrated into the NaV1.4 channel is computationally investigated. The formation of a charge transfer state, caused by electron transfer from the protein to the photoswitches, is observed. This state is strongly red shifted when the interaction takes place in a face-to-face orientation and electron donating groups are present on the aromatic ring of the amino acids. The low-energy charge transfer state can interfere with the photoisomerization process after excitation to the bright state by leading to the formation of radical species.
Collapse
Affiliation(s)
- Vito F Palmisano
- Department of Chemistry, Universidad Autónoma de Madrid, 28049, Madrid, Spain.
- Theoretical Chemistry Group, Zernike Institute for Advanced Materials, University of Groningen, Groningen, The Netherlands.
| | - Shirin Faraji
- Theoretical Chemistry Group, Zernike Institute for Advanced Materials, University of Groningen, Groningen, The Netherlands.
| | - Juan J Nogueira
- Department of Chemistry, Universidad Autónoma de Madrid, 28049, Madrid, Spain.
- Institute for Advanced Research in Chemistry (IAdChem), Universidad Autónoma de Madrid, Madrid 28049, Spain
| |
Collapse
|
2
|
Cárdenas G, Lucia‐Tamudo J, Mateo‐delaFuente H, Palmisano VF, Anguita‐Ortiz N, Ruano L, Pérez‐Barcia Á, Díaz‐Tendero S, Mandado M, Nogueira JJ. MoBioTools: A toolkit to setup quantum mechanics/molecular mechanics calculations. J Comput Chem 2023; 44:516-533. [PMID: 36507763 PMCID: PMC10107847 DOI: 10.1002/jcc.27018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2022] [Revised: 09/16/2022] [Accepted: 09/20/2022] [Indexed: 12/15/2022]
Abstract
We present a toolkit that allows for the preparation of QM/MM input files from a conformational ensemble of molecular geometries. The package is currently compatible with trajectory and topology files in Amber, CHARMM, GROMACS and NAMD formats, and has the possibility to generate QM/MM input files for Gaussian (09 and 16), Orca (≥4.0), NWChem and (Open)Molcas. The toolkit can be used in command line, so that no programming experience is required, although it presents some features that can also be employed as a python application programming interface. We apply the toolkit in four situations in which different electronic-structure properties of organic molecules in the presence of a solvent or a complex biological environment are computed: the reduction potential of the nucleobases in acetonitrile, an energy decomposition analysis of tyrosine interacting with water, the absorption spectrum of an azobenzene derivative integrated into a voltage-gated ion channel, and the absorption and emission spectra of the luciferine/luciferase complex. These examples show that the toolkit can be employed in a manifold of situations for both the electronic ground state and electronically excited states. It also allows for the automatic correction of the active space in the case of CASSCF calculations on an ensemble of geometries, as it is shown for the azobenzene derivative photoswitch case.
Collapse
Affiliation(s)
- Gustavo Cárdenas
- Department of ChemistryUniversidad Autónoma de MadridMadridSpain
| | | | | | | | | | - Lorena Ruano
- Department of ChemistryUniversidad Autónoma de MadridMadridSpain
| | | | - Sergio Díaz‐Tendero
- Department of ChemistryUniversidad Autónoma de MadridMadridSpain
- Institute for Advanced Research in Chemistry (IAdChem)Universidad Autónoma de MadridMadridSpain
- Condensed Matter Physics Center (IFIMAC)Universidad Autónoma de MadridMadridSpain
| | - Marcos Mandado
- Department of Physical ChemistryUniversity of VigoVigoSpain
| | - Juan J. Nogueira
- Department of ChemistryUniversidad Autónoma de MadridMadridSpain
- Institute for Advanced Research in Chemistry (IAdChem)Universidad Autónoma de MadridMadridSpain
| |
Collapse
|
3
|
Hégely B, Szirmai ÁB, Mester D, Tajti A, Szalay PG, Kállay M. Performance of Multilevel Methods for Excited States. J Phys Chem A 2022; 126:6548-6557. [PMID: 36095318 PMCID: PMC9511572 DOI: 10.1021/acs.jpca.2c05013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
The performance of multilevel quantum chemical approaches,
which
utilize an atom-based system partitioning scheme to model various
electronic excited states, is studied. The considered techniques include
the mechanical-embedding (ME) of “our own N-layered integrated
molecular orbital and molecular mechanics” (ONIOM) method,
the point charge embedding (PCE), the electronic-embedding (EE) of
ONIOM, the frozen density-embedding (FDE), the projector-based embedding
(PbE), and our local domain-based correlation method. For the investigated
multilevel approaches, the second-order algebraic-diagrammatic construction
[ADC(2)] approach was utilized as the high-level method, which was
embedded in either Hartree–Fock or a density functional environment.
The XH-27 test set of Zech et al. [J. Chem. Theory Comput., 2018, 14, 402829906111] was used for the assessment,
where organic dyes interact with several solvent molecules. With the
selection of the chromophores as active subsystems, we conclude that
the most reliable approach is local domain-based ADC(2) [L-ADC(2)],
and the least robust schemes are ONIOM-ME and ONIOM-EE. The PbE, FDE,
and PCE techniques often approach the accuracy of the L-ADC(2) scheme,
but their precision is far behind. The results suggest that a more
conservative subsystem selection algorithm or the inclusion of subsystem
charge-transfers is required for the atom-based cost-efficient methods
to produce high-accuracy excitation energies.
Collapse
Affiliation(s)
- Bence Hégely
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3, H-1111 Budapest, Hungary.,ELKH-BME Quantum Chemistry Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary
| | - Ádám B Szirmai
- Laboratory of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518 Budapest 112, Hungary
| | - Dávid Mester
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3, H-1111 Budapest, Hungary.,ELKH-BME Quantum Chemistry Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary
| | - Attila Tajti
- Laboratory of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518 Budapest 112, Hungary
| | - Péter G Szalay
- Laboratory of Theoretical Chemistry, Institute of Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518 Budapest 112, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3, H-1111 Budapest, Hungary.,ELKH-BME Quantum Chemistry Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary
| |
Collapse
|
4
|
Lucia-Tamudo J, Cárdenas G, Anguita-Ortiz N, Díaz-Tendero S, Nogueira JJ. Computation of Oxidation Potentials of Solvated Nucleobases by Static and Dynamic Multilayer Approaches. J Chem Inf Model 2022; 62:3365-3380. [PMID: 35771991 PMCID: PMC9326891 DOI: 10.1021/acs.jcim.2c00234] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
The determination
of the redox properties of nucleobases is of
paramount importance to get insight into the charge-transfer processes
in which they are involved, such as those occurring in DNA-inspired
biosensors. Although many theoretical and experimental studies have
been conducted, the value of the one-electron oxidation potentials
of nucleobases is not well-defined. Moreover, the most appropriate
theoretical protocol to model the redox properties has not been established
yet. In this work, we have implemented and evaluated different static
and dynamic approaches to compute the one-electron oxidation potentials
of solvated nucleobases. In the static framework, two thermodynamic
cycles have been tested to assess their accuracy against the direct
determination of oxidation potentials from the adiabatic ionization
energies. Then, the introduction of vibrational sampling, the effect
of implicit and explicit solvation models, and the application of
the Marcus theory have been analyzed through dynamic methods. The
results revealed that the static direct determination provides more
accurate results than thermodynamic cycles. Moreover, the effect of
sampling has not shown to be relevant, and the results are improved
within the dynamic framework when the Marcus theory is applied, especially
in explicit solvent, with respect to the direct approach. Finally,
the presence of different tautomers in water does not affect significantly
the one-electron oxidation potentials.
Collapse
Affiliation(s)
- Jesús Lucia-Tamudo
- Department of Chemistry, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Gustavo Cárdenas
- Department of Chemistry, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Nuria Anguita-Ortiz
- Department of Chemistry, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Sergio Díaz-Tendero
- Department of Chemistry, Universidad Autónoma de Madrid, 28049 Madrid, Spain.,Institute for Advanced Research in Chemistry (IAdChem), Universidad Autónoma de Madrid, 28049 Madrid, Spain.,Condensed Matter Physics Center (IFIMAC), Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Juan J Nogueira
- Department of Chemistry, Universidad Autónoma de Madrid, 28049 Madrid, Spain.,Institute for Advanced Research in Chemistry (IAdChem), Universidad Autónoma de Madrid, 28049 Madrid, Spain
| |
Collapse
|
5
|
Avagliano D, Bonfanti M, Garavelli M, González L. QM/MM Nonadiabatic Dynamics: the SHARC/COBRAMM Approach. J Chem Theory Comput 2021; 17:4639-4647. [PMID: 34114454 DOI: 10.1021/acs.jctc.1c00318] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present the SHARC/COBRAMM approach to enable easy and efficient excited-state dynamics simulations at different levels of electronic structure theory in the presence of complex environments using a quantum mechanics/molecular mechanics (QM/MM) setup. SHARC is a trajectory surface-hoping method that can incorporate the simultaneous effects of nonadiabatic and spin-orbit couplings in the excited-state dynamics of molecular systems. COBRAMM allows ground- and excited-state QM/MM calculations using a subtractive scheme, with electrostatic embedding and a hydrogen link-atom approach. The combination of both free and open-source program packages provides a modular and extensive framework to model nonadiabatic processes after light irradiation from the atomistic scale to the nano-scale. As an example, the relaxation of acrolein from S1 to T1 in solution is provided.
Collapse
Affiliation(s)
- Davide Avagliano
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1180 Vienna, Austria
| | - Matteo Bonfanti
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Viale Del Risorgimento, 4, I-40136 Bologna, Italy
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Viale Del Risorgimento, 4, I-40136 Bologna, Italy
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1180 Vienna, Austria.,Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Straße 17, A-1180 Vienna, Austria
| |
Collapse
|
6
|
Kostjukova LO, Leontieva SV, Kostjukov VV. Vibronic absorption spectrum and electronic properties of methylene blue in aqueous solution: TD-DFT study. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.116369] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
|
7
|
Yan Z, Li X, Chung LW. Multiscale Quantum Refinement Approaches for Metalloproteins. J Chem Theory Comput 2021; 17:3783-3796. [PMID: 34032440 DOI: 10.1021/acs.jctc.1c00148] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Biomolecules with metal ion(s) (e.g., metalloproteins) play many important biological roles. However, accurate structural determination of metalloproteins, particularly those containing transition metal ion(s), is challenging due to their complicated electronic structure, complex bonding of metal ions, and high number of conformations in biomolecules. Quantum refinement, which was proposed to combine crystallographic data with computational chemistry methods by several groups, can improve the local structures of some proteins. In this study, a quantum refinement method combining several multiscale computational schemes with experimental (X-ray diffraction) information was developed for metalloproteins. Various quantum refinement approaches using different ONIOM (our own N-layered integrated molecular orbital and molecular mechanics) combinations of quantum mechanics (QM), semiempirical (SE), and molecular mechanics (MM) methods were conducted to assess the performance and reliability on the refined local structure in two metalloproteins. The structures for two (Cu- or Zn-containing) metalloproteins were refined by combining two-layer ONIOM2(QM1/QM2) and ONIOM2(QM/MM) and three-layer ONIOM3(QM1/QM2/MM) schemes with experimental data. The accuracy of the quantum-refined metal binding sites was also examined and compared in these multiscale quantum refinement calculations. ONIOM3(QM/SE/MM) schemes were found to give good results with lower computational costs and were proposed to be a good choice for the multiscale computational scheme for quantum refinement calculations of metal binding site(s) in metalloproteins with high efficiency. Additionally, a two-center ONIOM approach was employed to speed up the quantum refinement calculations for the Zn metalloprotein with two remote active sites/ligands. Moreover, a recent quantum-embedding wavefunction-in-density functional theory (WF-in-DFT) method was also adopted as the high-level method in unprecedented ONIOM2(CCSD-in-B3LYP/MM) and ONIOM3(CCSD-in-B3LYP/SE/MM) calculations, which can be regarded as novel pseudo-three- and pseudo-four-layer ONIOM methods, respectively, to refine the key Zn binding site at the coupled-cluster singles and doubles (CCSD) level. These refined results indicate that multiscale quantum refinement schemes can be used to improve the structural accuracy obtained for local metal binding site(s) in metalloproteins with high efficiency.
Collapse
Affiliation(s)
- Zeyin Yan
- Shenzhen Grubbs Institute, Department of Chemistry and Guangdong Provincial Key Laboratory of Catalysis, Southern University of Science and Technology, Shenzhen 518055, China
| | - Xin Li
- Shenzhen Grubbs Institute, Department of Chemistry and Guangdong Provincial Key Laboratory of Catalysis, Southern University of Science and Technology, Shenzhen 518055, China
| | - Lung Wa Chung
- Shenzhen Grubbs Institute, Department of Chemistry and Guangdong Provincial Key Laboratory of Catalysis, Southern University of Science and Technology, Shenzhen 518055, China
| |
Collapse
|
8
|
Borrego-Sánchez A, Zemmouche M, Carmona-García J, Francés-Monerris A, Mulet P, Navizet I, Roca-Sanjuán D. Multiconfigurational Quantum Chemistry Determinations of Absorption Cross Sections (σ) in the Gas Phase and Molar Extinction Coefficients (ε) in Aqueous Solution and Air-Water Interface. J Chem Theory Comput 2021; 17:3571-3582. [PMID: 33974417 PMCID: PMC8444339 DOI: 10.1021/acs.jctc.0c01083] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2020] [Indexed: 11/29/2022]
Abstract
Theoretical determinations of absorption cross sections (σ) in the gas phase and molar extinction coefficients (ε) in condensed phases (water solution, interfaces or surfaces, protein or nucleic acids embeddings, etc.) are of interest when rates of photochemical processes, J = ∫ ϕ(λ) σ(λ) I(λ) dλ, are needed, where ϕ(λ) and I(λ) are the quantum yield of the process and the irradiance of the light source, respectively, as functions of the wavelength λ. Efficient computational strategies based on single-reference quantum-chemistry methods have been developed enabling determinations of line shapes or, in some cases, achieving rovibrational resolution. Developments are however lacking for strongly correlated problems, with many excited states, high-order excitations, and/or near degeneracies between states of the same and different spin multiplicities. In this work, we define and compare the performance of distinct computational strategies using multiconfigurational quantum chemistry, nuclear sampling of the chromophore (by means of molecular dynamics, ab initio molecular dynamics, or Wigner sampling), and conformational and statistical sampling of the environment (by means of molecular dynamics). A new mathematical approach revisiting previous absolute orientation algorithms is also developed to improve alignments of geometries. These approaches are benchmarked through the nπ* band of acrolein not only in the gas phase and water solution but also in a gas-phase/water interface, a common situation for instance in atmospheric chemistry. Subsequently, the best strategy is used to compute the absorption band for the adduct formed upon addition of an OH radical to the C6 position of uracil and compared with the available experimental data. Overall, quantum Wigner sampling of the chromophore with molecular dynamics sampling of the environment with CASPT2 electronic-structure determinations arise as a powerful methodology to predict meaningful σ(λ) and ε(λ) band line shapes with accurate absolute intensities.
Collapse
Affiliation(s)
- Ana Borrego-Sánchez
- Instituto
Andaluz de Ciencias de la Tierra, CSIC-University
of Granada, Av. de las
Palmeras 4, 18100 Armilla, Granada, Spain
| | - Madjid Zemmouche
- MSME,
Univ Gustave Eiffel, CNRS UMR 8208, Univ Paris-Est Créteil 8208, F-77454 Marne-la-Vallée, France
| | - Javier Carmona-García
- Instituto
de Ciencia Molecular, Universitat de València, P.O. Box 22085, València, Spain
| | - Antonio Francés-Monerris
- Université
de Lorraine and CNRS, LPCT UMR 7019, F-54000 Nancy, France
- Departamento
de Química Física, Universitat
de València, C/Dr.
Moliner 50, 46100 Burjassot, Spain
| | - Pep Mulet
- Departamento
de Matemáticas Área de Matemática Aplicada Facultad
de Matemáticas C/Dr. Moliner, 50 46100 Burjassot, Spain
| | - Isabelle Navizet
- MSME,
Univ Gustave Eiffel, CNRS UMR 8208, Univ Paris-Est Créteil 8208, F-77454 Marne-la-Vallée, France
| | - Daniel Roca-Sanjuán
- Instituto
de Ciencia Molecular, Universitat de València, P.O. Box 22085, València, Spain
| |
Collapse
|
9
|
Cárdenas G, Nogueira JJ. Stacking Effects on Anthraquinone/DNA Charge-Transfer Electronically Excited States. Molecules 2020; 25:E5927. [PMID: 33333751 PMCID: PMC7765225 DOI: 10.3390/molecules25245927] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2020] [Revised: 12/03/2020] [Accepted: 12/13/2020] [Indexed: 12/14/2022] Open
Abstract
The design of more efficient photosensitizers is a matter of great importance in the field of cancer treatment by means of photodynamic therapy. One of the main processes involved in the activation of apoptosis in cancer cells is the oxidative stress on DNA once a photosensitizer is excited by light. As a consequence, it is very relevant to investigate in detail the binding modes of the chromophore with DNA, and the nature of the electronically excited states that participate in the induction of DNA damage, for example, charge-transfer states. In this work, we investigate the electronic structure of the anthraquinone photosensitizer intercalated into a double-stranded poly(dG-dC) decamer model of DNA. First, the different geometric configurations are analyzed by means of classical molecular dynamics simulations. Then, the excited states for the most relevant poses of anthraquinone inside the binding pocket are computed by an electrostatic-embedding quantum mechanics/molecular mechanics approach, where anthraquinone and one of the nearby guanine residues are described quantum mechanically to take into account intermolecular charge-transfer states. The excited states are characterized as monomer, exciton, excimer, and charge-transfer states based on the analysis of the transition density matrix, and each of these contributions to the total density of states and absorption spectrum is discussed in terms of the stacking interactions. These results are relevant as they represent the footing for future studies on the reactivity of anthraquinone derivatives with DNA and give insights on possible geometrical configurations that potentially favor the oxidative stress of DNA.
Collapse
Affiliation(s)
- Gustavo Cárdenas
- Chemistry Department, Universidad Autónoma de Madrid, Calle Francisco Tomás y Valiente, 7, 28049 Madrid, Spain;
| | - Juan J. Nogueira
- Chemistry Department, Universidad Autónoma de Madrid, Calle Francisco Tomás y Valiente, 7, 28049 Madrid, Spain;
- IADCHEM, Institute for Advanced Research in Chemistry, Universidad Autónoma de Madrid, Calle Francisco Tomás y Valiente, 7, 28049 Madrid, Spain
| |
Collapse
|
10
|
Mai S, González L. Molecular Photochemistry: Recent Developments in Theory. Angew Chem Int Ed Engl 2020; 59:16832-16846. [PMID: 32052547 PMCID: PMC7540682 DOI: 10.1002/anie.201916381] [Citation(s) in RCA: 77] [Impact Index Per Article: 19.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2019] [Revised: 02/12/2020] [Indexed: 12/16/2022]
Abstract
Photochemistry is a fascinating branch of chemistry that is concerned with molecules and light. However, the importance of simulating light-induced processes is reflected also in fields as diverse as biology, material science, and medicine. This Minireview highlights recent progress achieved in theoretical chemistry to calculate electronically excited states of molecules and simulate their photoinduced dynamics, with the aim of reaching experimental accuracy. We focus on emergent methods and give selected examples that illustrate the progress in recent years towards predicting complex electronic structures with strong correlation, calculations on large molecules, describing multichromophoric systems, and simulating non-adiabatic molecular dynamics over long time scales, for molecules in the gas phase or in complex biological environments.
Collapse
Affiliation(s)
- Sebastian Mai
- Photonics InstituteVienna University of TechnologyGusshausstrasse 27–291040ViennaAustria
| | - Leticia González
- Institute of Theoretical ChemistryFaculty of ChemistryUniversity of ViennaWähringer Strasse 171090ViennaAustria
| |
Collapse
|
11
|
Mališ M, Luber S. Trajectory Surface Hopping Nonadiabatic Molecular Dynamics with Kohn–Sham ΔSCF for Condensed-Phase Systems. J Chem Theory Comput 2020; 16:4071-4086. [DOI: 10.1021/acs.jctc.0c00372] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Momir Mališ
- University of Zurich, Department of Chemistry, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Sandra Luber
- University of Zurich, Department of Chemistry, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| |
Collapse
|
12
|
Mai S, González L. Molekulare Photochemie: Moderne Entwicklungen in der theoretischen Chemie. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.201916381] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Sebastian Mai
- Institut für Photonik Technische Universität Wien Gußhausstraße 27–29 1040 Wien Österreich
| | - Leticia González
- Institut für theoretische Chemie Fakultät für Chemie Universität Wien Währinger Straße 17 1090 Wien Österreich
| |
Collapse
|
13
|
Assessing Configurational Sampling in the Quantum Mechanics/Molecular Mechanics Calculation of Temoporfin Absorption Spectrum and Triplet Density of States. Molecules 2018; 23:molecules23112932. [PMID: 30424014 PMCID: PMC6278509 DOI: 10.3390/molecules23112932] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2018] [Revised: 11/05/2018] [Accepted: 11/06/2018] [Indexed: 12/20/2022] Open
Abstract
The absorption properties of Temoporfin, a second-generation photosensitizer employed in photodynamic therapy, are calculated with an electrostatic-embedding quantum mechanics/molecular mechanics (QM/MM) scheme in methanol. The suitability of several ensembles of geometries generated by different sampling techniques, namely classical-molecular-dynamics (MD) and QM/MM-MD thermal sampling, Wigner quantum sampling and a hybrid protocol, which combines the thermal and quantum approaches, is assessed. It is found that a QM description of the chromophore during the sampling is needed in order to achieve a good agreement with respect to the experimental spectrum. Such a good agreement is obtained with both QM/MM-MD and Wigner samplings, demonstrating that a proper description of the anharmonic motions of the chromophore is not relevant in the computation of the absorption properties. In addition, it is also found that solvent organization is a rather fast process and a long sampling is not required. Finally, it is also demonstrated that the same exchange-correlation functional should be employed in the sampling and in the computation of the excited states properties to avoid unphysical triplet states with relative energies close or below 0 eV.
Collapse
|