1
|
Iuzzolino G, Perrella F, Valadan M, Petrone A, Altucci C, Rega N. Photophysics of a nucleic acid-protein crosslinking model strongly depends on solvation dynamics: an experimental and theoretical study. Phys Chem Chem Phys 2024; 26:11755-11769. [PMID: 38563904 DOI: 10.1039/d3cp06254f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
We present a combined experimental and theoretical study of the photophysics of 5-benzyluracil (5BU) in methanol, which is a model system for interactions between nucleic acids and proteins. A molecular dynamics study of 5BU in solution through efficient DFT-based hybrid ab initio potentials revealed a remarkable conformational flexibility - allowing the population of two main conformers - as well as specific solute-solvent interactions, which both appear as relevant factors for the observed 5BU optical absorption properties. The simulated absorption spectrum, calculated on such an ensemble, enabled a molecular interpretation of the experimental UV-Vis lowest energy band, which is also involved in the induced photo-reactivity upon irradiation. In particular, the first two excited states (mainly involving the uracil moiety) both contribute to the 5BU lowest energy absorption. Moreover, as a key finding, the nature and brightness of such electronic transitions are strongly influenced by 5BU conformation and the microsolvation of its heteroatoms.
Collapse
Affiliation(s)
- Gabriele Iuzzolino
- Scuola Superiore Meridionale, Largo San Marcellino 10, Napoli I-80138, Italy
- Dipartimento di Scienze Chimiche, Università degli Studi di Napoli Federico II, via Cintia 21, Napoli I-80126, Italy.
| | - Fulvio Perrella
- Scuola Superiore Meridionale, Largo San Marcellino 10, Napoli I-80138, Italy
| | - Mohammadhassan Valadan
- Dipartimento di Scienze Biomediche Avanzate, Università degli Studi di Napoli Federico II, via Pansini 5, Napoli I-80131, Italy
- Istituto Nazionale di Fisica Nucleare, Unità di Napoli, via Cintia 21, Napoli I-80126, Italy
| | - Alessio Petrone
- Scuola Superiore Meridionale, Largo San Marcellino 10, Napoli I-80138, Italy
- Dipartimento di Scienze Chimiche, Università degli Studi di Napoli Federico II, via Cintia 21, Napoli I-80126, Italy.
- Istituto Nazionale di Fisica Nucleare, Unità di Napoli, via Cintia 21, Napoli I-80126, Italy
| | - Carlo Altucci
- Dipartimento di Scienze Biomediche Avanzate, Università degli Studi di Napoli Federico II, via Pansini 5, Napoli I-80131, Italy
- Istituto Nazionale di Fisica Nucleare, Unità di Napoli, via Cintia 21, Napoli I-80126, Italy
- Istituto di Scienze Applicate e Sistemi Intelligenti "Eduardo Caianiello", URT UNINA, via Cintia 21, Napoli I-80126, Italy
| | - Nadia Rega
- Scuola Superiore Meridionale, Largo San Marcellino 10, Napoli I-80138, Italy
- Dipartimento di Scienze Chimiche, Università degli Studi di Napoli Federico II, via Cintia 21, Napoli I-80126, Italy.
- Istituto Nazionale di Fisica Nucleare, Unità di Napoli, via Cintia 21, Napoli I-80126, Italy
| |
Collapse
|
2
|
Cerezo J, Santoro F. FCclasses3: Vibrationally-resolved spectra simulated at the edge of the harmonic approximation. J Comput Chem 2023; 44:626-643. [PMID: 36380723 PMCID: PMC10100349 DOI: 10.1002/jcc.27027] [Citation(s) in RCA: 38] [Impact Index Per Article: 38.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Revised: 10/03/2022] [Accepted: 10/16/2022] [Indexed: 11/17/2022]
Abstract
We introduce FCclasses3, a code to carry out vibronic simulations of electronic spectra and nonradiative rates, based on the harmonic approximation. Key new features are: implementation of the full family of vertical and adiabatic harmonic models, vibrational analysis in curvilinear coordinates, extension to several electronic spectroscopies and implementation of time-dependent approaches. The use of curvilinear valence internal coordinates allows the adoption of quadratic model potential energy surfaces (PES) of the initial and final states expanded at arbitrary configurations. Moreover, the implementation of suitable projectors provides a robust framework for defining reduced-dimensionality models by sorting flexible coordinates out of the harmonic subset, so that they can then be treated at anharmonic level, or with mixed quantum classical approaches. A set of tools to facilitate input preparation and output analysis is also provided. We show the program at work in the simulation of different spectra (one and two-photon absorption, emission and resonance Raman) and internal conversion rate of a typical rigid molecule, anthracene. Then, we focus on absorption and emission spectra of a series of flexible polyphenyl molecules, highlighting the relevance of some of the newly implemented features. The code is freely available at http://www.iccom.cnr.it/en/fcclasses/.
Collapse
Affiliation(s)
- Javier Cerezo
- Departamento de Química and Institute for Advanced Research in Chemical Sciences (IAdChem)Universidad Autónoma de MadridMadridSpain
- Consiglio Nazionale delle RicercheIstituto di Chimica dei Composti Organo Metallici (ICCOM‐CNR)PisaItaly
| | - Fabrizio Santoro
- Consiglio Nazionale delle RicercheIstituto di Chimica dei Composti Organo Metallici (ICCOM‐CNR)PisaItaly
| |
Collapse
|
3
|
Aranda D, García-González F, Avila Ferrer FJ, López-Tocón I, Soto J, Otero JC. Computational Model for Electrochemical Surface-Enhanced Raman Scattering: Key Role of the Surface Charges and Synergy between Electromagnetic and Charge-Transfer Enhancement Mechanisms. J Chem Theory Comput 2022; 18:6802-6815. [PMID: 36222738 DOI: 10.1021/acs.jctc.2c00633] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a computational model for electrochemical surface-enhanced Raman scattering (EC-SERS). The surface excess of charge induced by the electrode potential (Vel) was introduced by applying an external electric field to a set of clusters [Agn]q with (n, q) of (19, ±1) or (20, 0) on which a molecule adsorbs. Using DFT/TD-DFT calculations, these metal-molecule complexes were classified by the adsorbate partial charge, and the main Vel-dependent properties were simultaneously studied with the aid of vibronic resonance Raman computations, namely, changes on the vibrational wavenumbers, relative intensities, and enhancement factors (EFs) for all SERS mechanisms: chemical or nonresonant, and resonance Raman with bright states of the adsorbate, charge-transfer (CT) states, and plasmon-like excitations on the metal cluster. We selected two molecules to test our model, pyridine, for which Vel has a remarkable effect, and 9,10-bis((E)-2-(pyridin-4-yl)vinyl)anthracene, which is almost insensitive to the applied bias. The results nicely reproduced most of the experimental observations, while the limitations of our approach were critically evaluated. We detected that accounting explicitly for the surface charges is key for EC-SERS models and that the highest calculated EFs, up to 107 to 108, are obtained by interstate coupling of bright local excitations of the metal cluster and CT states. These results highlight the importance of nonadiabatic effects in SERS and the capabilities of EC-SERS as a technique with potential to study excited-state coupling by tuning the CT and plasmon-like states by manipulating Vel.
Collapse
Affiliation(s)
- Daniel Aranda
- Instituto de Ciencia Molecular (ICMol), Universidad de Valencia, Catedrático J. Beltrán 2, 46980Paterna, Valencia, Spain.,Andalucía Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071Málaga, Spain
| | - Francisco García-González
- Andalucía Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071Málaga, Spain
| | - Francisco José Avila Ferrer
- Andalucía Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071Málaga, Spain
| | - Isabel López-Tocón
- Andalucía Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071Málaga, Spain
| | - Juan Soto
- Andalucía Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071Málaga, Spain
| | - Juan Carlos Otero
- Andalucía Tech, Facultad de Ciencias, Departamento de Química Física, Universidad de Málaga, 29071Málaga, Spain
| |
Collapse
|
4
|
Sun Y, Yuan CM, Xu SY, Li Y, Yang XB, Wen TC, Zhou K. An Unprecedented Tetracyclic Diterpenoid with Anti-inflammatory Activity from Papiliomyces sp. Tetrahedron Lett 2022. [DOI: 10.1016/j.tetlet.2022.153680] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
|
5
|
Segalina A, Cerezo J, Prampolini G, Santoro F, Pastore M. Accounting for Vibronic Features through a Mixed Quantum-Classical Scheme: Structure, Dynamics, and Absorption Spectra of a Perylene Diimide Dye in Solution. J Chem Theory Comput 2020; 16:7061-7077. [DOI: 10.1021/acs.jctc.0c00919] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Affiliation(s)
- Alekos Segalina
- Université de Lorraine & CNRS, LPCT, UMR 7019, F-54000 Nancy, France
| | - Javier Cerezo
- Departamento de Química and Institute for Advanced Research in Chemical Sciences (IAdChem), Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti Organo Metallici, Consiglio Nazionale delle Ricerche, (ICCOM-CNR) SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Fabrizio Santoro
- Istituto di Chimica dei Composti Organo Metallici, Consiglio Nazionale delle Ricerche, (ICCOM-CNR) SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | | |
Collapse
|
6
|
Giovannini T, Egidi F, Cappelli C. Theory and algorithms for chiroptical properties and spectroscopies of aqueous systems. Phys Chem Chem Phys 2020; 22:22864-22879. [PMID: 33043930 DOI: 10.1039/d0cp04027d] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Chiroptical properties and spectroscopies are valuable tools to study chiral molecules and assign absolute configurations. The spectra that result from chiroptical measurements may be very rich and complex, and hide much of their information content. For this reason, the interplay between experiments and calculations is especially useful, provided that all relevant physico-chemical interactions that are present in the experimental sample are accurately modelled. The inherent difficulty associated to the calculation of chiral signals of systems in aqueous solutions requires the development of specific tools, able to account for the peculiarities of water-solute interactions, and especially its ability to form hydrogen bonds. In this perspective we discuss a multiscale approach, which we have developed and challenged to model the most used chiroptical techniques.
Collapse
Affiliation(s)
- Tommaso Giovannini
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), 7491 Trondheim, Norway
| | | | | |
Collapse
|
7
|
Raucci U, Perrella F, Donati G, Zoppi M, Petrone A, Rega N. Ab-initio molecular dynamics and hybrid explicit-implicit solvation model for aqueous and nonaqueous solvents: GFP chromophore in water and methanol solution as case study. J Comput Chem 2020; 41:2228-2239. [PMID: 32770577 DOI: 10.1002/jcc.26384] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2020] [Revised: 06/21/2020] [Accepted: 06/27/2020] [Indexed: 12/15/2022]
Abstract
Solute-solvent interactions are proxies for understanding how the electronic density of a chromophore interacts with the environment in a more exhaustive way. The subtle balance between polarization, electrostatic, and non-bonded interactions need to be accurately described to obtain good agreement between simulations and experiments. First principles approaches providing accurate configurational sampling through molecular dynamics may be a suitable choice to describe solvent effects on solute chemical-physical properties and spectroscopic features, such as optical absorption of dyes. In this context, accurate energy potentials, obtained by hybrid implicit/explicit solvation methods along with employing nonperiodic boundary conditions, are required to represent bulk solvent around a large solute-solvent cluster. In this work, a novel strategy to simulate methanol solutions is proposed combining ab initio molecular dynamics, a hybrid implicit/explicit flexible solvent model, nonperiodic boundary conditions, and time dependent density functional theory. As case study, the robustness of the proposed protocol has been gauged by investigating the microsolvation and electronic absorption of the anionic green fluorescent protein chromophore in methanol and aqueous solution. Satisfactory results are obtained, reproducing the microsolvation layout of the chromophore and, as a consequence, the experimental trends shown by the optical absorption in different solvents.
Collapse
Affiliation(s)
- Umberto Raucci
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
| | - Fulvio Perrella
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
| | - Greta Donati
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy.,Dipartimento di Chimica e Biologia "Adolfo Zambelli", Università di Salerno, Fisciano, Italy
| | - Maria Zoppi
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
| | - Alessio Petrone
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
| | - Nadia Rega
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy.,Center for Advanced Biomaterials for Healthcare@CRIB, Naples, Italy
| |
Collapse
|
8
|
Del Galdo S, Fusè M, Barone V. CPL Spectra of Camphor Derivatives in Solution by an Integrated QM/MD Approach. Front Chem 2020; 8:584. [PMID: 32733856 PMCID: PMC7358700 DOI: 10.3389/fchem.2020.00584] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2020] [Accepted: 06/05/2020] [Indexed: 11/16/2022] Open
Abstract
We extend a recently proposed computational strategy for the simulation of absorption spectra of semi-rigid molecular systems in condensed phases to the emission spectra of flexible chromophores. As a case study, we have chosen the CPL spectrum of camphor in methanol solution, which shows a well-defined bisignate shape. The first step of our approach is the quantum mechanical computation of reference spectra including vibrational averaging effects and taking bulk solvent effects into account by means of the polarizable continuum model. In the present case, the large amplitude inversion mode is explicitly treated by a numerical approach, whereas the other small-amplitude vibrational modes are taken into account within the harmonic approximation. Next, the snapshots of classical molecular dynamics computations are clusterized and one representative configuration from each cluster is used to compute a reference spectrum. In the present case, different clusters correspond to the two stable conformers of camphor in the S1 excited electronic state and, for each of them, to different numbers of strong solute-solvent hydrogen bonds. Finally, local fluctuation effects within each cluster are taken into account by means of the perturbed matrix model. The overall procedure leads to good agreement with experiment for absorption and emission spectra together with their chiral counterparts, thus allowing to analyze the role of different effects (stereo-electronic, vibrational, environmental) in tuning the overall experimental spectra.
Collapse
Affiliation(s)
| | - Marco Fusè
- SMART Laboratory, Scuola Normale Superiore, Pisa, Italy
| | | |
Collapse
|
9
|
Aranda D, Santoro F. The role of chlorine position in the electronic circular dichroism of chlorophenyl-ethanol investigated by vibronic calculations. Chirality 2020; 32:932-948. [PMID: 32445432 DOI: 10.1002/chir.23239] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2020] [Revised: 04/24/2020] [Accepted: 04/24/2020] [Indexed: 11/08/2022]
Abstract
(R)-1-phenyl-ethanol (PhEtOH) and the different isomers of (R)-1-(chlorophenyl)ethanol (ClPhEtOH) exhibit very interesting electronic circular dichroism (ECD) in methanol. In all cases, the spectrum shows clear vibronic features, but it is monosignated and negative for PhEtOH and meta-ClPhEtOH, positive for the ortho isomer and bisignated for the para isomer. We used computational chemistry to rationalise this behaviour adopting CAM-B3LYP/def2-TZVP, describing the bulk solvent effects with polarizable continuum models and solute-solvent specific interactions with clusters comprising the solute and two solvent molecules. We adopted harmonic vibronic models to compute the ECD spectral shapes of all stable conformers, and we obtained the room-temperature spectra by Boltzmann average. Simulated spectra are in very good agreement with experiment and allow us to rationalise their difference in terms of the relevance of Franck-Condon (FC) and Herzberg-Teller (HT) intensity-borrowing contributions, modulated by the substituent effect. The bisignated shape of the spectrum of para-ClPhEtOH arises from the competition of opposite-sign FC and HT bands, promoted by different vibrational modes. Due to the challenges we document in computing its ECD spectrum, para-ClPhEtOH represents a good test case to help the development of novel methodologies for an improved description of weak vibronic ECD spectra of flexible systems in explicit solvents.
Collapse
Affiliation(s)
- Daniel Aranda
- Consiglio Nazionale delle Ricerche-Istituto di Chimica dei Composti Organometallici (ICCOM-CNR), Pisa, Italy.,Departamento de Química Física, Facultad de Ciencias, Universidad de Málaga, Málaga, Spain
| | - Fabrizio Santoro
- Consiglio Nazionale delle Ricerche-Istituto di Chimica dei Composti Organometallici (ICCOM-CNR), Pisa, Italy
| |
Collapse
|
10
|
Cerezo J, Aranda D, Avila Ferrer FJ, Prampolini G, Santoro F. Adiabatic-Molecular Dynamics Generalized Vertical Hessian Approach: A Mixed Quantum Classical Method To Compute Electronic Spectra of Flexible Molecules in the Condensed Phase. J Chem Theory Comput 2020; 16:1215-1231. [PMID: 31855424 DOI: 10.1021/acs.jctc.9b01009] [Citation(s) in RCA: 45] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
We present a general mixed quantum classical method that couples classical molecular dynamics (MD) and vibronic models to compute the shape of electronic spectra of flexible molecules in the condensed phase without, in principle, any phenomenological broadening. It is based on a partition of the nuclear motions of the solute + solvent system in "soft" and "stiff" vibrational modes and an adiabatic hypothesis that assumes that stiff modes are much faster than soft ones. In this framework, the spectrum is rigorously expressed as a conformational integral of quantum vibronic spectra along the stiff coordinates only. Soft modes enter at the classical level through the conformational distribution that is sampled with classical MD runs. In each configuration, reduced-dimensionality quadratic Hamiltonians are built in the space of the stiff coordinates only, thanks to a generalization of the Vertical Hessian harmonic model and an iterative application of projectors in internal coordinates to remove soft modes. Quantum vibronic spectra, specific for each sampled configuration of the soft coordinates, are then computed at the desired temperature with efficient time-dependent techniques, and the global spectrum simply arises from their average. For consistency of the whole procedure, classical MD runs are performed with quantum-mechanically derived force fields, parameterized at the same level of theory selected for generating the quadratic Hamiltonians along the stiff coordinates. Application to N-methyl-6-oxyquinolinium betaine in water, dithiophene in ethanol, and cyanidine in water is presented to show the performance of the method.
Collapse
Affiliation(s)
| | - Daniel Aranda
- CNR-Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa , Area della Ricerca, via G. Moruzzi 1 , I-56124 Pisa , Italy.,Departamento de Química Física, Facultad de Ciencias, Andalucía Tech , Universidad de Málaga , E-29071 Málaga , Spain
| | - Francisco José Avila Ferrer
- Departamento de Química Física, Facultad de Ciencias, Andalucía Tech , Universidad de Málaga , E-29071 Málaga , Spain
| | - Giacomo Prampolini
- CNR-Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa , Area della Ricerca, via G. Moruzzi 1 , I-56124 Pisa , Italy
| | - Fabrizio Santoro
- CNR-Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa , Area della Ricerca, via G. Moruzzi 1 , I-56124 Pisa , Italy
| |
Collapse
|
11
|
Baiardi A, Bloino J, Barone V. Time-Dependent Formulation of Resonance Raman Optical Activity Spectroscopy. J Chem Theory Comput 2018; 14:6370-6390. [PMID: 30281300 DOI: 10.1021/acs.jctc.8b00488] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
In this work, we extend the theoretical framework recently developed for the simulation of resonance Raman (RR) spectra of medium-to-large sized systems to its chiral counterpart, namely, resonance Raman optical activity (RROA). The theory is based on a time-dependent (TD) formulation, with the transition tensors obtained as half-Fourier transforms of the appropriate cross-correlation functions. The implementation has been kept as general as possible, supporting adiabatic and vertical models for the PES representation, both in Cartesian and internal coordinates, with the possible inclusion of Herzberg-Teller (HT) effects. Thanks to the integration of this TD-RROA procedure within a general-purpose quantum-chemistry program, both solvation and leading anharmonicity effects can be included in an effective way. The implementation is validated on one of the smallest chiral molecule (methyloxirane). Practical applications are illustrated with three medium-size organic molecules (naproxen-OCD3, quinidine and 2-Br-hexahelicene), whose simulated spectra are compared to the corresponding experimental data.
Collapse
Affiliation(s)
- Alberto Baiardi
- Scuola Normale Superiore , piazza dei Cavalieri 7 , I-56126 Pisa , Italy
| | - Julien Bloino
- Scuola Normale Superiore , piazza dei Cavalieri 7 , I-56126 Pisa , Italy
| | - Vincenzo Barone
- Scuola Normale Superiore , piazza dei Cavalieri 7 , I-56126 Pisa , Italy
| |
Collapse
|
12
|
Aranda D, Cerezo J, Pescitelli G, Avila Ferrer FJ, Soto J, Santoro F. A computational study of the vibrationally-resolved electronic circular dichroism spectra of single-chain transoid and cisoid oligothiophenes in chiral conformations. Phys Chem Chem Phys 2018; 20:21864-21880. [PMID: 30105334 DOI: 10.1039/c8cp03482f] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We simulate the vibronic profile of the electronic circular dichroism (ECD) spectra of oligothiophenes in cisoid and transoid chiral arrangements. We consider oligomers of different lengths, from two to fifteen units, and investigate extensively how the ECD spectral shapes depend on the inter-ring torsions. In general, the molecular structures we consider are not stationary points of the ground state potential energy surface. Therefore, in order to perform vibronic calculations, we present a new computational protocol able to define reduced-dimensionality models where the effect of the off-equilibrium modes is removed. This is done adopting a description of the vibrational motions in curvilinear internal coordinates, and vertical harmonic models coupled with an iterative application of projectors to define energy Hessians, and therefore effective normal modes, in the space complementary to the one of the off-equilibrium coordinates. Although we consider both Franck-Condon and Herzberg-Teller contributions, the results show that transoid twisted ribbons always give rise to monosignated ECD spectra, while bi-signated and multi-signated spectra are expected for cisoid helices. These findings are explained on the basis of the different transition strengths of the lowest electronic states imparted by the different spatial arrangement, that is almost linear for transoid structures and more globular for cisoid ones. We predicted the chiroptical response of a large number of possible molecular arrangements. These data are employed to critically discuss the experimental ECD of polythiophenes in different experimental conditions, forming either aggregates or host-guest complexes. The method here proposed to perform vibronic calculations in reduced-dimensionality models is of general applicability and its potential interest goes beyond the practical application presented here.
Collapse
Affiliation(s)
- Daniel Aranda
- Department of Physical Chemistry, Faculty of Science, University of Málaga, E-29071-Málaga, Spain
| | | | | | | | | | | |
Collapse
|
13
|
Zuehlsdorff TJ, Napoli JA, Milanese JM, Markland TE, Isborn CM. Unraveling electronic absorption spectra using nuclear quantum effects: Photoactive yellow protein and green fluorescent protein chromophores in water. J Chem Phys 2018; 149:024107. [PMID: 30007372 DOI: 10.1063/1.5025517] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Many physical phenomena must be accounted for to accurately model solution-phase optical spectral line shapes, from the sampling of chromophore-solvent configurations to the electronic-vibrational transitions leading to vibronic fine structure. Here we thoroughly explore the role of nuclear quantum effects, direct and indirect solvent effects, and vibronic effects in the computation of the optical spectrum of the aqueously solvated anionic chromophores of green fluorescent protein and photoactive yellow protein. By analyzing the chromophore and solvent configurations, the distributions of vertical excitation energies, the absorption spectra computed within the ensemble approach, and the absorption spectra computed within the ensemble plus zero-temperature Franck-Condon approach, we show how solvent, nuclear quantum effects, and vibronic transitions alter the optical absorption spectra. We find that including nuclear quantum effects in the sampling of chromophore-solvent configurations using ab initio path integral molecular dynamics simulations leads to improved spectral shapes through three mechanisms. The three mechanisms that lead to line shape broadening and a better description of the high-energy tail are softening of heavy atom bonds in the chromophore that couple to the optically bright state, widening the distribution of vertical excitation energies from more diverse solvation environments, and redistributing spectral weight from the 0-0 vibronic transition to higher energy vibronic transitions when computing the Franck-Condon spectrum in a frozen solvent pocket. The absorption spectra computed using the combined ensemble plus zero-temperature Franck-Condon approach yield significant improvements in spectral shape and width compared to the spectra computed with the ensemble approach. Using the combined approach with configurations sampled from path integral molecular dynamics trajectories presents a significant step forward in accurately modeling the absorption spectra of aqueously solvated chromophores.
Collapse
Affiliation(s)
- Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Joseph A Napoli
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Joel M Milanese
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
14
|
Cerezo J, Aranda D, Avila Ferrer FJ, Prampolini G, Mazzeo G, Longhi G, Abbate S, Santoro F. Toward a general mixed quantum/classical method for the calculation of the vibronic ECD of a flexible dye molecule with different stable conformers: Revisiting the case of 2,2,2-trifluoro-anthrylethanol. Chirality 2018; 30:730-743. [DOI: 10.1002/chir.22853] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2018] [Revised: 02/19/2018] [Accepted: 02/20/2018] [Indexed: 12/30/2022]
Affiliation(s)
- Javier Cerezo
- Departamento de Química Física; Universidad de Murcia; Murcia Spain
| | - Daniel Aranda
- Consiglio Nazionale delle Ricerche-Istituto di Chimica dei Composti Organometallici (ICCOM-CNR); Pisa Italy
- Departamento de Química Física; Universidad de Málaga; Málaga Spain
| | | | - Giacomo Prampolini
- Consiglio Nazionale delle Ricerche-Istituto di Chimica dei Composti Organometallici (ICCOM-CNR); Pisa Italy
| | - Giuseppe Mazzeo
- Dipartimento di Medicina Molecolare e Traslazionale; Università di Brescia; Brescia Italy
| | - Giovanna Longhi
- Dipartimento di Medicina Molecolare e Traslazionale; Università di Brescia; Brescia Italy
| | - Sergio Abbate
- Dipartimento di Medicina Molecolare e Traslazionale; Università di Brescia; Brescia Italy
| | - Fabrizio Santoro
- Consiglio Nazionale delle Ricerche-Istituto di Chimica dei Composti Organometallici (ICCOM-CNR); Pisa Italy
| |
Collapse
|
15
|
Marazzi M, Gattuso H, Monari A, Assfeld X. Steady-State Linear and Non-linear Optical Spectroscopy of Organic Chromophores and Bio-macromolecules. Front Chem 2018; 6:86. [PMID: 29666792 PMCID: PMC5891624 DOI: 10.3389/fchem.2018.00086] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2018] [Accepted: 03/12/2018] [Indexed: 01/05/2023] Open
Abstract
Bio-macromolecules as DNA, lipid membranes and (poly)peptides are essential compounds at the core of biological systems. The development of techniques and methodologies for their characterization is therefore necessary and of utmost interest, even though difficulties can be experienced due to their intrinsic complex nature. Among these methods, spectroscopies, relying on optical properties are especially important to determine their macromolecular structures and behaviors, as well as the possible interactions and reactivity with external dyes—often drugs or pollutants—that can (photo)sensitize the bio-macromolecule leading to eventual chemical modifications, thus damages. In this review, we will focus on the theoretical simulation of electronic spectroscopies of bio-macromolecules, considering their secondary structure and including their interaction with different kind of (photo)sensitizers. Namely, absorption, emission and electronic circular dichroism (CD) spectra are calculated and compared with the available experimental data. Non-linear properties will be also taken into account by two-photon absorption, a highly promising technique (i) to enhance absorption in the red and infra-red windows and (ii) to enhance spatial resolution. Methodologically, the implications of using implicit and explicit solvent, coupled to quantum and thermal samplings of the phase space, will be addressed. Especially, hybrid quantum mechanics/molecular mechanics (QM/MM) methods are explored for a comparison with solely QM methods, in order to address the necessity to consider an accurate description of environmental effects on spectroscopic properties of biological systems.
Collapse
Affiliation(s)
- Marco Marazzi
- Laboratoire de Physique et Chimie Théoriques, Université de Lorraine-Nancy, UMR 7019, Vandoeuvre-lés-Nancy, France.,Laboratoire de Physique et Chimie Théoriques, Centre National de la Recherche Scientifique, UMR 7019, Vandoeuvre-lès-Nancy, France.,Departamento de Química, Centro de Investigacíon en Síntesis Química (CISQ), Universidad de La Rioja, Logroño, Spain
| | - Hugo Gattuso
- Laboratoire de Physique et Chimie Théoriques, Université de Lorraine-Nancy, UMR 7019, Vandoeuvre-lés-Nancy, France.,Laboratoire de Physique et Chimie Théoriques, Centre National de la Recherche Scientifique, UMR 7019, Vandoeuvre-lès-Nancy, France
| | - Antonio Monari
- Laboratoire de Physique et Chimie Théoriques, Université de Lorraine-Nancy, UMR 7019, Vandoeuvre-lés-Nancy, France.,Laboratoire de Physique et Chimie Théoriques, Centre National de la Recherche Scientifique, UMR 7019, Vandoeuvre-lès-Nancy, France
| | - Xavier Assfeld
- Laboratoire de Physique et Chimie Théoriques, Université de Lorraine-Nancy, UMR 7019, Vandoeuvre-lés-Nancy, France.,Laboratoire de Physique et Chimie Théoriques, Centre National de la Recherche Scientifique, UMR 7019, Vandoeuvre-lès-Nancy, France
| |
Collapse
|
16
|
Baiardi A, Bloino J, Barone V. Simulation of Vibronic Spectra of Flexible Systems: Hybrid DVR-Harmonic Approaches. J Chem Theory Comput 2017; 13:2804-2822. [PMID: 28467703 DOI: 10.1021/acs.jctc.7b00236] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Our general framework for the simulation of vibrational signatures in electronic spectra has been extended to treat one large-amplitude motion (LAM) at the anharmonic level, coupled to the other small-amplitude motions (SAM) treated as harmonic. The coupling between LAM and SAM is minimized thanks to the use of delocalized internal coordinates, which are built automatically from the molecular topology. General LAMs can be employed, ranging from intrinsic reaction coordinates to rigid or flexible paths based on the distinguished coordinate approach. The anharmonic model is based on a fully numerical method based on the discrete variable representation (DVR) theory, supporting different types of boundary conditions. The inclusion of this model in a general-purpose electronic structure code makes available to the user a large panel of quantum chemistry models, for both isolated and condensed phases. The flexibility and reliability of the new framework are illustrated by some case studies, covering various types of LAMs, ranging from a small test case, the photoelectron spectrum of ammonia, to larger systems, such as phenylanthracene and cyclobutanone.
Collapse
Affiliation(s)
- Alberto Baiardi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.,Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), UOS di Pisa , Area della Ricerca CNR, Via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
17
|
Padula D, Cerezo J, Pescitelli G, Santoro F. The shape of the electronic circular dichroism spectrum of (2,6-dimethylphenyl)(phenyl)methanol: interplay between conformational equilibria and vibronic effects. Phys Chem Chem Phys 2017; 19:32349-32360. [DOI: 10.1039/c7cp06369e] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Analysis of the interplay between conformational equilibria, solvent effects and vibronic contributions in the ECD spectra.
Collapse
Affiliation(s)
- Daniele Padula
- Department of Chemistry
- University of Liverpool
- Liverpool L69 7ZD
- UK
| | - Javier Cerezo
- Departamento de Química Física
- Universidad de Murcia
- 30100 Murcia
- Spain
| | - Gennaro Pescitelli
- Università di Pisa
- Dipartimento di Chimica e Chimica Industriale
- I-56124 Pisa
- Italy
| | - Fabrizio Santoro
- Consiglio Nazionale delle Ricerche – CNR
- Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR)
- Pisa
- Italy
| |
Collapse
|