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
|
Falbo E, Fusè M, Lazzari F, Mancini G, Barone V. Integration of Quantum Chemistry, Statistical Mechanics, and Artificial Intelligence for Computational Spectroscopy: The UV-Vis Spectrum of TEMPO Radical in Different Solvents. J Chem Theory Comput 2022; 18:6203-6216. [PMID: 36166322 PMCID: PMC9558374 DOI: 10.1021/acs.jctc.2c00654] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Indexed: 11/30/2022]
Abstract
The ongoing integration of quantum chemistry, statistical mechanics, and artificial intelligence is paving the route toward more effective and accurate strategies for the investigation of the spectroscopic properties of medium-to-large size chromophores in condensed phases. In this context we are developing a novel workflow aimed at improving the generality, reliability, and ease of use of the available computational tools. In this paper we report our latest developments with specific reference to unsupervised atomistic simulations employing non periodic boundary conditions (NPBC) followed by clustering of the trajectories employing optimized feature spaces. Next accurate variational computations are performed for a representative point of each cluster, whereas intracluster fluctuations are taken into account by a cheap yet reliable perturbative approach. A number of methodological improvements have been introduced including, e.g., more realistic reaction field effects at the outer boundary of the simulation sphere, automatic definition of the feature space by continuous perception of solute-solvent interactions, full account of polarization and charge transfer in the first solvation shell, and inclusion of vibronic contributions. After its validation, this new approach has been applied to the challenging case of solvatochromic effects on the UV-vis spectra of a prototypical nitroxide radical (TEMPO) in different solvents. The reliability, effectiveness, and robustness of the new platform is demonstrated by the remarkable agreement with experiment of the results obtained through an unsupervised approach characterized by a strongly reduced computational cost as compared to that of conventional quantum mechanics and molecular mechanics models without any accuracy reduction.
Collapse
Affiliation(s)
- Emanuele Falbo
- Scuola
Normale Superiore di Pisa, piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Marco Fusè
- Scuola
Normale Superiore di Pisa, piazza dei Cavalieri 7, 56126 Pisa, Italy
- Dipartimento
di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123 Brescia, Italy
| | - Federico Lazzari
- Scuola
Normale Superiore di Pisa, piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Giordano Mancini
- Scuola
Normale Superiore di Pisa, piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola
Normale Superiore di Pisa, piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
3
|
Sagresti L, Peri L, Ceccarelli G, Brancato G. Stochastic Model of Solvent Exchange in the First Coordination Shell of Aqua Ions. J Chem Theory Comput 2022; 18:3164-3173. [PMID: 35471007 PMCID: PMC9097284 DOI: 10.1021/acs.jctc.2c00181] [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
Ion microsolvation is a basic, yet fundamental, process of ionic solutions underlying many relevant phenomena in either biological or nanotechnological applications, such as solvent reorganization energy, ion transport, catalytic activity, and so on. As a consequence, it is a topic of extensive investigations by various experimental techniques, ranging from X-ray diffraction to NMR relaxation and from calorimetry to vibrational spectroscopy, and theoretical approaches, especially those based on molecular dynamics (MD) simulations. The conventional microscopic view of ion solvation is usually provided by a "static" cluster model representing the first ion-solvent coordination shell. Despite the merits of such a simple model, however, ion coordination in solution should be better regarded as a complex population of dynamically interchanging molecular configurations. Such a more comprehensive view is more subtle to characterize and often elusive to standard approaches. In this work, we report on an effective computational strategy aiming at providing a detailed picture of solvent coordination and exchange around aqua ions, thus including the main structural, thermodynamic, and dynamic properties of ion microsolvation, such as the most probable first-shell complex structures, the corresponding free energies, the interchanging energy barriers, and the solvent-exchange rates. Assuming the solvent coordination number as an effective reaction coordinate and combining MD simulations with enhanced sampling and master-equation approaches, we propose a stochastic model suitable for properly describing, at the same time, the thermodynamics and kinetics of ion-water coordination. The model is successfully tested toward various divalent ions (Ca2+, Zn2+, Hg2+, and Cd2+) in aqueous solution, considering also the case of a high ionic concentration. Results show a very good agreement with those issuing from brute-force MD simulations, when available, and support the reliable prediction of rare ion-water complexes and slow water exchange rates not easily accessible to usual computational methods.
Collapse
Affiliation(s)
- Luca Sagresti
- Scuola Normale Superiore, Piazza Dei Cavalieri 7, I-56126 Pisa, Italy.,Istituto Nazionale di Fisica Nucleare (INFN), Largo Pontecorvo 3, I-56127 Pisa, Italy
| | - Lorenzo Peri
- Scuola Normale Superiore, Piazza Dei Cavalieri 7, I-56126 Pisa, Italy
| | - Giacomo Ceccarelli
- Dipartimento di Fisica, Università di Pisa, Largo Bruno Pontecorvo 3, I-56127 Pisa, Italy
| | - Giuseppe Brancato
- Scuola Normale Superiore, Piazza Dei Cavalieri 7, I-56126 Pisa, Italy.,Istituto Nazionale di Fisica Nucleare (INFN), Largo Pontecorvo 3, I-56127 Pisa, Italy.,Consorzio Interuniversitario per Lo Sviluppo Dei Sistemi a Grande Interfase (CSGI), Via Della Lastruccia 3, I-50019 Sesto Fiorentino (FI), Italy
| |
Collapse
|
4
|
Barone V, Puzzarini C, Mancini G. Integration of theory, simulation, artificial intelligence and virtual reality: a four-pillar approach for reconciling accuracy and interpretability in computational spectroscopy. Phys Chem Chem Phys 2021; 23:17079-17096. [PMID: 34346437 DOI: 10.1039/d1cp02507d] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The established pillars of computational spectroscopy are theory and computer based simulations. Recently, artificial intelligence and virtual reality are becoming the third and fourth pillars of an integrated strategy for the investigation of complex phenomena. The main goal of the present contribution is the description of some new perspectives for computational spectroscopy, in the framework of a strategy in which computational methodologies at the state of the art, high-performance computing, artificial intelligence and virtual reality tools are integrated with the aim of improving research throughput and achieving goals otherwise not possible. Some of the key tools (e.g., continuous molecular perception model and virtual multifrequency spectrometer) and theoretical developments (e.g., non-periodic boundaries, joint variational-perturbative models) are shortly sketched and their application illustrated by means of representative case studies taken from recent work by the authors. Some of the results presented are already well beyond the state of the art in the field of computational spectroscopy, thereby also providing a proof of concept for other research fields.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| | | | | |
Collapse
|
5
|
Chiariello MG, Raucci U, Donati G, Rega N. Water-Mediated Excited State Proton Transfer of Pyranine-Acetate in Aqueous Solution: Vibrational Fingerprints from Ab Initio Molecular Dynamics. J Phys Chem A 2021; 125:3569-3578. [PMID: 33900071 PMCID: PMC8279639 DOI: 10.1021/acs.jpca.1c00692] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
In this work, we
simulate the excited state proton transfer (ESPT)
reaction involving the pyranine photoacid and an acetate molecule
as proton acceptor, connected by a bridge water molecule. We employ
ab initio molecular dynamics combined with an hybrid quantum/molecular
mechanics (QM/MM) framework. Furthermore, a time-resolved vibrational
analysis based on the wavelet-transform allows one to identify two
low frequency vibrational modes that are fingerprints of the ESPT
event: a ring wagging and ring breathing. Their composition suggests
their key role in optimizing the structure of the proton donor–acceptor
couple and promoting the ESPT event. We find that the choice of the
QM/MM partition dramatically affects the photoinduced reactivity of
the system. The QM subspace was gradually extended including the water
molecules directly interacting with the pyranine–water–acetate
system. Indeed, the ESPT reaction takes place when the hydrogen bond
network around the reactive system is taken into account at full QM
level.
Collapse
Affiliation(s)
- Maria Gabriella Chiariello
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, via Cintia, I-80126 Napoli, Italy
| | - Umberto Raucci
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, via Cintia, I-80126 Napoli, Italy
| | - Greta Donati
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, via Cintia, I-80126 Napoli, Italy
| | - Nadia Rega
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S. Angelo, via Cintia, I-80126 Napoli, Italy.,Centro Interdipartimentale di Ricerca sui Biomateriali (CRIB) Piazzale Tecchio, Largo Barsanti e Matteucci, I-80125 Napoli, Italy
| |
Collapse
|
6
|
Raucci U, Chiariello MG, Rega N. Modeling Excited-State Proton Transfer to Solvent: A Dynamics Study of a Super Photoacid with a Hybrid Implicit/Explicit Solvent Model. J Chem Theory Comput 2020; 16:7033-7043. [PMID: 33112132 PMCID: PMC8016186 DOI: 10.1021/acs.jctc.0c00782] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
![]()
The rapid growth of time-resolved
spectroscopies and the theoretical
advances in ab initio molecular dynamics (AIMD) pave the way to look
at the real-time molecular motion following the electronic excitation.
Here, we exploited the capabilities of AIMD combined with a hybrid
implicit/explicit model of solvation to investigate the ultrafast
excited-state proton transfer (ESPT) reaction of a super photoacid,
known as QCy9, in water solution. QCy9 transfers a proton to a water
solvent molecule within 100 fs upon the electronic excitation in aqueous
solution, and it is the strongest photoacid reported in the literature
so far. Because of the ultrafast kinetics, it has been experimentally
hypothesized that the ESPT escapes the solvent dynamics control (Huppert
et al., J. Photochem. Photobiol. A2014,277, 90). The sampling of the solvent configuration
space on the ground electronic state is the first key step toward
the simulation of the ESPT event. Therefore, several configurations
in the Franck–Condon region, describing an average solvation,
were chosen as starting points for the excited-state dynamics. In
all cases, the excited-state evolution spontaneously leads to the
proton transfer event, whose rate is strongly dependent on the hydrogen
bond network around the proton acceptor solvent molecule. Our study
revealed that the explicit representation at least of three solvation
shells around the proton acceptor molecule is necessary to stabilize
the excess proton. Furthermore, the analysis of the solvent molecule
motions in proximity of the reaction site suggested that even in the
case of the strongest photoacid, the ESPT is actually assisted by
the solvation dynamics of the first and second solvation shells of
the water accepting molecule.
Collapse
Affiliation(s)
- Umberto Raucci
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S.Angelo, via Cintia, I-80126 Napoli, Italy
| | - Maria Gabriella Chiariello
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S.Angelo, via Cintia, I-80126 Napoli, Italy
| | - Nadia Rega
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S.Angelo, via Cintia, I-80126 Napoli, Italy.,CRIB, Centro Interdipartimentale di Ricerca sui Biomateriali, Piazzale Tecchio, I-80125 Napoli, Italy
| |
Collapse
|
7
|
Donati G, Petrone A, Rega N. Multiresolution continuous wavelet transform for studying coupled solute-solvent vibrations via ab initio molecular dynamics. Phys Chem Chem Phys 2020; 22:22645-22661. [PMID: 33015693 DOI: 10.1039/d0cp02495c] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Vibrational analysis in solution and the theoretical determination of infrared and Raman spectra are of key importance in many fields of chemical interest. Vibrational band dynamics of molecules and their sensitivity to the environment can also be captured by these spectroscopies in their time dependent version. However, it is often difficult to provide an interpretation of the experimental data at the molecular scale, such as molecular mechanisms or the processes hidden behind them. In this work, we present a theoretical-computational protocol based on ab initio molecular dynamics simulations and a combination of normal-like (generalized) mode analysis of solute-solvent clusters with a wavelet transform, for the first time. The case study is the vibrational dynamics of N-methyl-acetamide (NMA) in water solution, a well-known model of hydration of peptides and proteins. Amide modes are typical bands of peptide and protein backbone, and their couplings with the environment are very challenging in terms of the accurate prediction of solvent induced intensity and frequency shifts. The contribution of water molecules surrounding NMA to the composition of generalized and time resolved modes is introduced in our vibrational analysis, showing unequivocally its influence on the amide mode spectra. It is also shown that such mode compositions need the inclusion of the first shell solvent molecules to be accurately described. The wavelet analysis is proven to be strongly recommended to follow the time evolution of the spectra, and to capture vibrational band couplings and frequency shifts over time, preserving at the same time a well-balanced time-frequency resolution. This peculiar feature also allows one to perform a combined structural-vibrational analysis, where the different strengths of hydrogen bond interactions can quantitatively affect the amide bands over time at finite temperature. The proposed method allows for the direct connection between vibrational modes and local structural changes, providing a link from the spectroscopic observable to the structure, in this case the peptide backbone, and its hydration layouts.
Collapse
Affiliation(s)
- Greta Donati
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M. S. Angelo, Via Cintia, I-80126 Napoli, Italy.
| | | | | |
Collapse
|
8
|
Chiariello MG, Donati G, Rega N. Time-Resolved Vibrational Analysis of Excited State Ab Initio Molecular Dynamics to Understand Photorelaxation: The Case of the Pyranine Photoacid in Aqueous Solution. J Chem Theory Comput 2020; 16:6007-6013. [PMID: 32955870 PMCID: PMC8011922 DOI: 10.1021/acs.jctc.0c00810] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
![]()
We
present a novel time-resolved vibrational analysis for studying
photoinduced nuclear relaxation. Generalized modes velocities are
defined from ab initio molecular dynamics and wavelet transformed,
providing the time localization of vibrational signals in the electronic
excited state. The photoexcited pyranine in aqueous solution is presented
as a case study. The transient and sequential activation of the simulated
vibrational signals is in good agreement with vibrational dynamics
obtained from femtosecond stimulated Raman spectroscopy data.
Collapse
Affiliation(s)
- Maria Gabriella Chiariello
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M. S. Angelo, via Cintia, I-80126 Napoli, Italy
| | - Greta Donati
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M. S. Angelo, via Cintia, I-80126 Napoli, Italy
| | - Nadia Rega
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M. S. Angelo, via Cintia, I-80126 Napoli, Italy.,CRIB Center for Advanced Biomaterials for Healthcare, Largo Barsanti e Matteucci, I-80125 Napoli, Italy
| |
Collapse
|
9
|
Mancini G, Del Galdo S, Chandramouli B, Pagliai M, Barone V. Computational Spectroscopy in Solution by Integration of Variational and Perturbative Approaches on Top of Clusterized Molecular Dynamics. J Chem Theory Comput 2020; 16:5747-5761. [PMID: 32697580 PMCID: PMC8009517 DOI: 10.1021/acs.jctc.0c00454] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
![]()
Multiscale QM/MM approaches have
become the most suitable and effective
methods for the investigation of spectroscopic properties of medium-
or large-size chromophores in condensed phases. On these grounds,
we are developing a novel workflow aimed at improving the generality,
reliability, and ease of use of the available tools. In the present
paper, we report the latest developments of such an approach with
specific reference to a general workplan starting with the addition
of acetonitrile to the panel of solvents already available in the
General Liquid Optimized Boundary (GLOB) model enforcing nonperiodic
boundary conditions (NPBC). Next, the solvatochromic shifts induced
by acetonitrile on both rigid (uracil and thymine) and flexible (thyrosine)
chromophores have been studied introducing in our software a number
of new features ranging from rigid-geometry NPBC molecular dynamics
based on the quaternion formalism to a full integration of variational
(ONIOM) and perturbative (perturbed matrix method (PMM)) approaches
for describing different solute–solvent topologies and local
fluctuations, respectively. Finally, thymine and uracil have been
studied also in methanol to point out the generality of the computational
strategy. While further developments are surely needed, the strengths
of our integrated approach even in its present version are demonstrated
by the accuracy of the results obtained by an unsupervised approach
and coupled to a computational cost strongly reduced with respect
to that of conventional QM/MM models without any appreciable accuracy
deterioration.
Collapse
Affiliation(s)
- Giordano Mancini
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.,Istituto Nazionale di Fisica Nucleare (INFN) sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
| | - Sara Del Galdo
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | | | - Marco Pagliai
- Dipartimento di Chimica "Ugo Schiff", Università degli Studi di Firenze, Via della Lastruccia 3, 50019 Sesto Fiorentino, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.,Istituto Nazionale di Fisica Nucleare (INFN) sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
| |
Collapse
|
10
|
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
|
11
|
Chiariello MG, Rega N. Exploring Nuclear Photorelaxation of Pyranine in Aqueous Solution: an Integrated Ab-Initio Molecular Dynamics and Time Resolved Vibrational Analysis Approach. J Phys Chem A 2018; 122:2884-2893. [DOI: 10.1021/acs.jpca.7b12371] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Maria Gabriella Chiariello
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S.Angelo, via Cintia, I-80126 Napoli, Italy
| | - Nadia Rega
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di M.S.Angelo, via Cintia, I-80126 Napoli, Italy
- Interdisciplinary Research Centre on Biomaterials (CRIB) Università di Napoli Federico II, Piazzale Tecchio 80, I-80125, Napoli, Italy
| |
Collapse
|
12
|
Petrone A, Williams-Young DB, Lingerfelt DB, Li X. Ab Initio Excited-State Transient Raman Analysis. J Phys Chem A 2017; 121:3958-3965. [DOI: 10.1021/acs.jpca.7b02905] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Alessio Petrone
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | | | - David B. Lingerfelt
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
13
|
Petrone A, Lingerfelt DB, Williams-Young DB, Li X. Ab Initio Transient Vibrational Spectral Analysis. J Phys Chem Lett 2016; 7:4501-4508. [PMID: 27788583 DOI: 10.1021/acs.jpclett.6b02292] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
Pump probe spectroscopy techniques have enabled the direct observation of a variety of transient molecular species in both ground and excited electronic states. Time-resolved vibrational spectroscopy is becoming an indispensable tool for investigating photoinduced nuclear dynamics of chemical systems of all kinds. On the other hand, a complete picture of the chemical dynamics encoded in these spectra cannot be achieved without a full temporal description of the structural relaxation, including the explicit time-dependence of vibrational coordinates that are substantially displaced from equilibrium by electronic excitation. Here we present a transient vibrational analysis protocol combining ab initio direct molecular dynamics and time-integrated normal modes introduced in this work, relying on the recent development of analytic time-dependent density functional theory (TDDFT) second derivatives for excited states. Prototypical molecules will be used as test cases, showing the evolution of the vibrational signatures that follow electronic excitation. This protocol provides a direct route to assigning the vibrations implicated in the (photo)dynamics of several (photoactive) systems.
Collapse
Affiliation(s)
- Alessio Petrone
- Department of Chemistry, University of Washington , Seattle, Washington 98195, United States
| | - David B Lingerfelt
- Department of Chemistry, University of Washington , Seattle, Washington 98195, United States
| | - David B Williams-Young
- Department of Chemistry, University of Washington , Seattle, Washington 98195, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington , Seattle, Washington 98195, United States
| |
Collapse
|
14
|
Šipr O, Vackář J, Kuzmin A. Effect of atomic vibrations in XANES: polarization-dependent damping of the fine structure at the Cu K-edge of (creat) 2CuCl 4. JOURNAL OF SYNCHROTRON RADIATION 2016; 23:1433-1439. [PMID: 27787249 DOI: 10.1107/s1600577516014570] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/02/2016] [Accepted: 09/13/2016] [Indexed: 06/06/2023]
Abstract
Polarization-dependent damping of the fine structure in the Cu K-edge spectrum of creatinium tetrachlorocuprate [(creat)2CuCl4] in the X-ray absorption near-edge structure (XANES) region is shown to be due to atomic vibrations. These vibrations can be separated into two groups, depending on whether the respective atoms belong to the same molecular block; individual molecular blocks can be treated as semi-rigid entities while the mutual positions of these blocks are subject to large mean relative displacements. The effect of vibrations can be efficiently included in XANES calculations by using the same formula as for static systems but with a modified free-electron propagator which accounts for fluctuations in interatomic distances.
Collapse
Affiliation(s)
- Ondřej Šipr
- Institute of Physics of the ASCR v.v.i., Cukrovarnicka 10, CZ-162 53 Prague, Czech Republic
| | - Jiří Vackář
- Institute of Physics of the ASCR v.v.i., Na Slovance 2, CZ-182 21 Prague, Czech Republic
| | - Alexei Kuzmin
- Institute of Solid State Physics, University of Latvia, Kengaraga Street 8, LV-1063 Riga, Latvia
| |
Collapse
|
15
|
|
16
|
Mancini G, Brancato G, Chandramouli B, Barone V. Organic solvent simulations under non-periodic boundary conditions: A library of effective potentials for the GLOB model. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2015.03.001] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
|
17
|
Thaunay F, Dognon JP, Ohanessian G, Clavaguéra C. Vibrational mode assignment of finite temperature infrared spectra using the AMOEBA polarizable force field. Phys Chem Chem Phys 2015. [DOI: 10.1039/c5cp02270c] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The Driven Molecular Dynamics approach has been adapted and associated with the AMOEBA polarizable force field to assign and visualize vibrational modes in infrared spectra obtained by molecular dynamics simulations.
Collapse
Affiliation(s)
- Florian Thaunay
- Laboratoire de Chimie Moléculaire
- Ecole Polytechnique
- CNRS
- 91128 Palaiseau Cedex
- France
| | - Jean-Pierre Dognon
- CEA/Saclay
- DSM/IRAMIS/NIMBE
- CNRS
- UMR 3685
- Laboratoire de Chimie Moléculaire et de Catalyse pour l'Energie
| | - Gilles Ohanessian
- Laboratoire de Chimie Moléculaire
- Ecole Polytechnique
- CNRS
- 91128 Palaiseau Cedex
- France
| | - Carine Clavaguéra
- Laboratoire de Chimie Moléculaire
- Ecole Polytechnique
- CNRS
- 91128 Palaiseau Cedex
- France
| |
Collapse
|
18
|
Petrone A, Donati G, Caruso P, Rega N. Understanding THz and IR Signals beneath Time-Resolved Fluorescence from Excited-State Ab Initio Dynamics. J Am Chem Soc 2014; 136:14866-74. [DOI: 10.1021/ja507518k] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Alessio Petrone
- Dipartimento
di Scienze Chimiche, Università di Napoli ‘Federico
II’, Complesso Universitario di M.S. Angelo, via Cintia, I-80126 Napoli, Italy
| | - Greta Donati
- Dipartimento
di Scienze Chimiche, Università di Napoli ‘Federico
II’, Complesso Universitario di M.S. Angelo, via Cintia, I-80126 Napoli, Italy
| | - Pasquale Caruso
- Italian Institute of Technology, IIT@CRIB Center for Advanced Biomaterials for Healthcare, Largo
Barsanti e Matteucci, I-80125 Napoli, Italy
| | - Nadia Rega
- Dipartimento
di Scienze Chimiche, Università di Napoli ‘Federico
II’, Complesso Universitario di M.S. Angelo, via Cintia, I-80126 Napoli, Italy
- Italian Institute of Technology, IIT@CRIB Center for Advanced Biomaterials for Healthcare, Largo
Barsanti e Matteucci, I-80125 Napoli, Italy
| |
Collapse
|
19
|
Lu X, Gaus M, Elstner M, Cui Q. Parametrization of DFTB3/3OB for magnesium and zinc for chemical and biological applications. J Phys Chem B 2014; 119:1062-82. [PMID: 25178644 PMCID: PMC4306495 DOI: 10.1021/jp506557r] [Citation(s) in RCA: 108] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
![]()
We report the parametrization of
the approximate density functional
theory, DFTB3, for magnesium and zinc for chemical and biological
applications. The parametrization strategy follows that established
in previous work that parametrized several key main group elements
(O, N, C, H, P, and S). This 3OB set of parameters can thus be used
to study many chemical and biochemical systems. The parameters are
benchmarked using both gas-phase and condensed-phase systems. The
gas-phase results are compared to DFT (mostly B3LYP), ab initio (MP2 and G3B3), and PM6, as well as to a previous DFTB parametrization
(MIO). The results indicate that DFTB3/3OB is particularly successful
at predicting structures, including rather complex dinuclear metalloenzyme
active sites, while being semiquantitative (with a typical mean absolute
deviation (MAD) of ∼3–5 kcal/mol) for energetics. Single-point
calculations with high-level quantum mechanics (QM) methods generally
lead to very satisfying (a typical MAD of ∼1 kcal/mol) energetic
properties. DFTB3/MM simulations for solution and two enzyme systems
also lead to encouraging structural and energetic properties in comparison
to available experimental data. The remaining limitations of DFTB3,
such as the treatment of interaction between metal ions and highly
charged/polarizable ligands, are also discussed.
Collapse
Affiliation(s)
- Xiya Lu
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison , 1101 University Avenue, Madison, Wisconsin 53706, United States
| | | | | | | |
Collapse
|
20
|
Barone V, Biczysko M, Bloino J. Fully anharmonic IR and Raman spectra of medium-size molecular systems: accuracy and interpretation. Phys Chem Chem Phys 2014; 16:1759-87. [PMID: 24346191 PMCID: PMC4604664 DOI: 10.1039/c3cp53413h] [Citation(s) in RCA: 295] [Impact Index Per Article: 29.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Computation of full infrared (IR) and Raman spectra (including absolute intensities and transition energies) for medium- and large-sized molecular systems beyond the harmonic approximation is one of the most interesting challenges of contemporary computational chemistry. Contrary to common beliefs, low-order perturbation theory is able to deliver results of high accuracy (actually often better than those issuing from current direct dynamics approaches) provided that anharmonic resonances are properly managed. This perspective sketches the recent developments in our research group toward the development of a robust and user-friendly virtual spectrometer rooted in second-order vibrational perturbation theory (VPT2) and usable also by non-specialists essentially as a black-box procedure. Several examples are explicitly worked out in order to illustrate the features of our computational tool together with the most important ongoing developments.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| | | | | |
Collapse
|
21
|
Riahi S, Roux B, Rowley CN. QM/MM molecular dynamics simulations of the hydration of Mg(II) and Zn(II) ions. CAN J CHEM 2013. [DOI: 10.1139/cjc-2012-0515] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
The hydration of Mg2+ and Zn2+ is examined using molecular dynamics simulations using 3 computational approaches of increasing complexity: the CHARMM nonpolarizable force field based on the TIP3P water model, the Drude polarizable force field based on the SWM4-NDP water model, and a combined QM/MM approach in which the inner coordination sphere is represented using a high-quality density functional theory (DFT) model (PBE/def2-TZVPP), and the remainder of the bulk water solvent is represented using the polarizable SWM4-NDP water model. The characteristic structural distribution functions (radial, angular, and tilt) are comparedand show very good agreement between the polarizable force field and QM/MM approaches. They predict an average Mg–O distance of 2.11 Å and an Zn–O distance of 2.13 Å, in good agreement with the available experimental neutron scattering and EXAFS data, while the Mg–O distances calculated using the nonpolarizable force field are 0.1 Å too short. Mg2+ (aq) and Zn2+ (aq) both have a coordination number of 6 and have a remarkably similar octahedral coordination mode, despite the chemical differences between these ions. Thermodynamic integration was used to calculate the relative hydration free energies of these ions (ΔΔGhydr). The nonpolarizable model is in error by 60 kcal mol– 1 and incorrectly predicts that Mg2+ has the more negative hydration energy. The Drude polarizable model predicts a ΔΔGhydr of only –13.2 kcal kcal mol– 1, an improvement over the results of the nonpolarizable force field, but still signficantly different than the experimental value of –30.1 kcal mol–1. The combined QM/MM approach performs much better, predicting a ΔΔGhydr of –34.8 kcal mol–1 in excellent agreement with experiment. These calculations support the experimental observation that Zn2+ has more favourable solvation free energy than Mg2+ despite having a very similar solvation structure.
Collapse
Affiliation(s)
- Saleh Riahi
- Department of Chemistry, Memorial University of Newfoundland, St. John's, NL A1B 3X7, Canada
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, The University of Chicago, 929 East 57th Street, Chicago, Ill., USA
| | - Christopher N. Rowley
- Department of Chemistry, Memorial University of Newfoundland, St. John's, NL A1B 3X7, Canada
| |
Collapse
|
22
|
Lipparini F, Barone V. Polarizable Force Fields and Polarizable Continuum Model: A Fluctuating Charges/PCM Approach. 1. Theory and Implementation. J Chem Theory Comput 2011; 7:3711-24. [DOI: 10.1021/ct200376z] [Citation(s) in RCA: 126] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|