1
|
Sepali C, Lafiosca P, Gómez S, Giovannini T, Cappelli C. Effective fully polarizable QM/MM approaches to compute Raman and Raman Optical Activity spectra in aqueous solution. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2024; 305:123485. [PMID: 37827000 DOI: 10.1016/j.saa.2023.123485] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/23/2023] [Revised: 09/22/2023] [Accepted: 10/01/2023] [Indexed: 10/14/2023]
Abstract
Raman and Raman Optical Activity (ROA) signals are amply affected by solvent effects, especially in the presence of strongly solute-solvent interactions such as Hydrogen Bonding (HB). In this work, we extend the fully atomistic polarizable Quantum Mechanics/Molecular Mechanics approach, based on the Fluctuating Charges and Fluctuating Dipoles force field to the calculation of Raman and ROA spectra. Such an approach is able to accurately describe specific HB interactions, by also accounting for anisotropic contributions due to the inclusion of fluctuating dipoles. To highlight the potentiality of the novel approach, Raman and ROA spectra of L-Serine and L-Cysteine dissolved in aqueous solution are computed and compared both with alternative theoretical approaches and experimental measurements.
Collapse
Affiliation(s)
- Chiara Sepali
- Scuola Normale Superiore, Piazza dei Cavalieri, 7, Pisa, 56126, Italy
| | - Piero Lafiosca
- Scuola Normale Superiore, Piazza dei Cavalieri, 7, Pisa, 56126, Italy
| | - Sara Gómez
- Scuola Normale Superiore, Piazza dei Cavalieri, 7, Pisa, 56126, Italy
| | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri, 7, Pisa, 56126, Italy.
| |
Collapse
|
2
|
Rojas-Valencia N, Gómez S, Giovannini T, Cappelli C, Restrepo A, Núñez Zarur F. Water Maintains the UV-Vis Spectral Features During the Insertion of Anionic Naproxen and Ibuprofen into Model Cell Membranes. J Phys Chem B 2023; 127:2146-2155. [PMID: 36877579 DOI: 10.1021/acs.jpcb.2c08332] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2023]
Abstract
UV-vis spectra of anionic ibuprofen and naproxen in a model lipid bilayer of the cell membrane are investigated using computational techniques in combination with a comparative analysis of drug spectra in purely aqueous environments. The simulations aim at elucidating the intricacies behind the negligible changes in the maximum absorption wavelength in the experimental spectra. A set of configurations of the systems constituted by lipid, water, and drugs or just water and drugs are obtained from classical Molecular Dynamics simulations. UV-vis spectra are computed in the framework of atomistic Quantum Mechanical/Molecular Mechanics (QM/MM) approaches together with Time-Dependent Density Functional Theory (TD-DFT). Our results suggest that the molecular orbitals involved in the electronic transitions are the same, regardless of the chemical environment. A thorough analysis of the contacts between the drug and water molecules reveals that no significant changes in UV-vis spectra are a consequence of ibuprofen and naproxen molecules being permanently microsolvated by water molecules, despite the presence of lipid molecules. Water molecules microsolvate the charged carboxylate group as expected but also microsolvate the aromatic regions of the drugs.
Collapse
Affiliation(s)
- Natalia Rojas-Valencia
- Facultad de Ciencias Básicas, Universidad de Medellín, Carrera 87 No. 30-65, 050026, Medellín, Colombia
| | - Sara Gómez
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126, Pisa, Italy
| | - Tommaso Giovannini
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126, Pisa, Italy
| | - Chiara Cappelli
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126, Pisa, Italy
| | - Albeiro Restrepo
- Instituto de Química, Universidad de Antioquia, UdeA, Calle 70 No. 52-21 050010, Medellín, Colombia
| | - Francisco Núñez Zarur
- Facultad de Ciencias Básicas, Universidad de Medellín, Carrera 87 No. 30-65, 050026, Medellín, Colombia
| |
Collapse
|
3
|
Xu Q, Liu Y, Wang M, Cerezo J, Improta R, Santoro F. The Resonance Raman Spectrum of Cytosine in Water: Analysis of the Effect of Specific Solute-Solvent Interactions and Non-Adiabatic Couplings. Molecules 2023; 28:molecules28052286. [PMID: 36903532 PMCID: PMC10005559 DOI: 10.3390/molecules28052286] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Revised: 02/21/2023] [Accepted: 02/22/2023] [Indexed: 03/06/2023] Open
Abstract
In this contribution, we report a computational study of the vibrational Resonance Raman (vRR) spectra of cytosine in water, on the grounds of potential energy surfaces (PES) computed by time-dependent density functional theory (TD-DFT) and CAM-B3LYP and PBE0 functionals. Cytosine is interesting because it is characterized by several close-lying and coupled electronic states, challenging the approach commonly used to compute the vRR for systems where the excitation frequency is in quasi-resonance with a single state. We adopt two recently developed time-dependent approaches, based either on quantum dynamical numerical propagations of vibronic wavepackets on coupled PES or on analytical correlation functions for cases in which inter-state couplings were neglected. In this way, we compute the vRR spectra, considering the quasi-resonance with the eight lowest-energy excited states, disentangling the role of their inter-state couplings from the mere interference of their different contributions to the transition polarizability. We show that these effects are only moderate in the excitation energy range explored by experiments, where the spectral patterns can be rationalized from the simple analysis of displacements of the equilibrium positions along the different states. Conversely, at higher energies, interference and inter-state couplings play a major role, and the adoption of a fully non-adiabatic approach is strongly recommended. We also investigate the effect of specific solute-solvent interactions on the vRR spectra, by considering a cluster of cytosine, hydrogen-bonded by six water molecules, and embedded in a polarizable continuum. We show that their inclusion remarkably improves the agreement with the experiments, mainly altering the composition of the normal modes, in terms of internal valence coordinates. We also document cases, mostly for low-frequency modes, in which a cluster model is not sufficient, and more elaborate mixed quantum classical approaches, in explicit solvent models, need to be applied.
Collapse
Affiliation(s)
- Qiushuang Xu
- School of Physics Engineering, Qufu Normal University, Qufu 273165, China
- School of Physics and Optoelectronics Engineering, Ludong University, Yantai 264025, China
- 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
| | - Yanli Liu
- School of Physics and Optoelectronics Engineering, Ludong University, Yantai 264025, China
| | - Meishan Wang
- School of Physics Engineering, Qufu Normal University, Qufu 273165, China
- School of Physics and Optoelectronics Engineering, Ludong University, Yantai 264025, China
| | - Javier Cerezo
- 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 and Institute for Advanced Research in Chemical Sciences (IAdChem), Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Roberto Improta
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), Via De Amicis 95, I-80145 Napoli, Italy
- Correspondence: (R.I.); (F.S.)
| | - Fabrizio Santoro
- 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
- Correspondence: (R.I.); (F.S.)
| |
Collapse
|
4
|
Gómez S, Lafiosca P, Egidi F, Giovannini T, Cappelli C. UV-Resonance Raman Spectra of Systems in Complex Environments: A Multiscale Modeling Applied to Doxorubicin Intercalated into DNA. J Chem Inf Model 2023; 63:1208-1217. [PMID: 36745496 PMCID: PMC9976284 DOI: 10.1021/acs.jcim.2c01495] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
UV-Resonance Raman (RR) spectroscopy is a valuable tool to study the binding of drugs to biomolecular receptors. The extraction of information at the molecular level from experimental RR spectra is made much easier and more complete thanks to the use of computational approaches, specifically tuned to deal with the complexity of the supramolecular system. In this paper, we propose a protocol to simulate RR spectra of complex systems at different levels of sophistication, by exploiting a quantum mechanics/molecular mechanics (QM/MM) approach. The approach is challenged to investigate RR spectra of a widely used chemotherapy drug, doxorubicin (DOX) intercalated into a DNA double strand. The computed results show good agreement with experimental data, thus confirming the reliability of the computational protocol.
Collapse
Affiliation(s)
- Sara Gómez
- Scuola
Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy,E-mail:
| | - Piero Lafiosca
- Scuola
Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Franco Egidi
- Software
for Chemistry and Materials BV, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Tommaso Giovannini
- Scuola
Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Chiara Cappelli
- Scuola
Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy,E-mail:
| |
Collapse
|
5
|
Gómez S, Giovannini T, Cappelli C. Multiple Facets of Modeling Electronic Absorption Spectra of Systems in Solution. ACS PHYSICAL CHEMISTRY AU 2022; 3:1-16. [PMID: 36718266 PMCID: PMC9881242 DOI: 10.1021/acsphyschemau.2c00050] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/30/2022] [Revised: 11/08/2022] [Accepted: 11/08/2022] [Indexed: 11/24/2022]
Abstract
In this Perspective, we outline the essential physicochemical aspects that need to be considered when building a reliable approach to describe absorption properties of solvated systems. In particular, we focus on how to properly model the complexity of the solvation phenomenon, arising from dynamical aspects and specific, strong solute-solvent interactions. To this end, conformational and configurational sampling techniques, such as Molecular Dynamics, have to be coupled to accurate fully atomistic Quantum Mechanical/Molecular Mechanics (QM/MM) methodologies. By exploiting different illustrative applications, we show that an effective reproduction of experimental spectral signals can be achieved by delicately balancing exhaustive sampling, hydrogen bonding, mutual polarization, and nonelectrostatic effects.
Collapse
|
6
|
Strunge K, Madzharova F, Jensen F, Weidner T, Nagata Y. Theoretical Sum Frequency Generation Spectra of Protein Amide with Surface-Specific Velocity-Velocity Correlation Functions. J Phys Chem B 2022; 126:8571-8578. [PMID: 36194760 DOI: 10.1021/acs.jpcb.2c04321] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Vibrational sum frequency generation (vSFG) spectroscopy is widely used to probe the protein structure at interfaces. Because protein vSFG spectra are complex, they can only provide detailed structural information if combined with computer simulations of protein molecular dynamics and spectra calculations. We show how vSFG spectra can be accurately modeled using a surface-specific velocity-velocity scheme based on ab initio normal modes. Our calculated vSFG spectra show excellent agreement with the experimental sum frequency spectrum of LTα14 peptide and provide insight into the origin of the characteristic α-helical amide I peak. Analysis indicates that the peak shape can be explained largely by two effects: (1) the uncoupled response of amide groups located on opposite sides of the α-helix will have different orientations with respect to the interface and therefore different local environments affecting the local mode vibrations and (2) vibrational splitting from nearest neighbor coupling evaluated as inter-residue vibrational correlation. The conclusion is consistent with frequency mapping techniques with an empirically based ensemble of peptide structures, thus showing how time correlation approaches and frequency mapping techniques can give independent yet complementary molecular descriptions of protein vSFG. These models reveal the sensitive relationship between protein structure and their amide I response, allowing exploitation of the complicated molecular vibrations and their interference to derive the structures of proteins under native conditions at interfaces.
Collapse
Affiliation(s)
- Kris Strunge
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Fani Madzharova
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Frank Jensen
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Tobias Weidner
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Yuki Nagata
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| |
Collapse
|
7
|
Xu Q, Aranda D, Yaghoubi Jouybari M, Liu Y, Wang M, Cerezo J, Improta R, Santoro F. Nonadiabatic Vibrational Resonance Raman Spectra from Quantum Dynamics Propagations with LVC Models. Application to Thymine. J Phys Chem A 2022; 126:7468-7479. [PMID: 36099554 PMCID: PMC9596142 DOI: 10.1021/acs.jpca.2c05271] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
We present a viable protocol to compute vibrational resonance
Raman
(vRR) spectra for systems with several close-lying and potentially
coupled electronic states. It is based on the parametrization of linear
vibronic coupling (LVC) models from time-dependent density functional
theory (TD-DFT) calculations and quantum dynamics propagations of
vibronic wavepackets with the multilayer version of the multiconfiguration
time-dependent Hartree (ML-MCTDH) method. Our approach is applied
to thymine considering seven coupled electronic states, comprising
the three lowest bright states, and all vibrational coordinates. Computed
vRR at different excitation wavelengths are in good agreement with
the available experimental data. Up to 250 nm the signal is dominated
by the lowest HOMO → LUMO transition, whereas at 233 nm, in
the valley between the two lowest energy absorption bands, the contributions
of all the three bright states, and their interferences and couplings,
are important. Inclusion of solvent (water) effects improves the agreement
with experiment, reproducing the coalescence of vibrational bands
due to CC and C=O stretchings. With our approach we disentangle
and assess the effect of interferences between the contribution of
different quasi-resonant states to the transition polarizability and
the effect of interstate couplings. Our findings strongly suggest
that in cases of close-lying and potentially coupled states a simple
inclusion of interference effects is not sufficient, and a fully nonadiabatic
computation should instead be performed. We also document that for
systems with strong couplings and quasi-degenerate states, the use
of HT perturbative approach, not designed for these cases, may lead
to large artifacts.
Collapse
Affiliation(s)
- Qiushuang Xu
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China.,School of Physics Engineering, Qufu Normal University, 2673100 Qufu, Shandong, PR China.,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
| | - Daniel Aranda
- Instituto de Ciencia Molecular (ICMol)., Universidad de Valencia, c/Catedrático José Beltrán, 2, 46980 Paterna, Spain
| | - Martha Yaghoubi Jouybari
- 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
| | - Yanli Liu
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China
| | - Meishan Wang
- School of Physics and Optoelectronics Engineering, Ludong University, 264025 Yantai, Shandong, PR China
| | - Javier Cerezo
- 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, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Roberto Improta
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), Via De Amicis 95, I-80145 Napoli, Italy.,DTU Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
| | - Fabrizio Santoro
- 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
|
8
|
Lafiosca P, Gómez S, Giovannini T, Cappelli C. Absorption Properties of Large Complex Molecular Systems: The DFTB/Fluctuating Charge Approach. J Chem Theory Comput 2022; 18:1765-1779. [PMID: 35184553 PMCID: PMC8908768 DOI: 10.1021/acs.jctc.1c01066] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
![]()
We report on the
first formulation of a novel polarizable QM/MM
approach, where the density functional tight binding (DFTB) is coupled
to the fluctuating charge (FQ) force field. The resulting method (DFTB/FQ)
is then extended to the linear response within the TD-DFTB framework
and challenged to study absorption spectra of large condensed-phase
systems.
Collapse
Affiliation(s)
- Piero Lafiosca
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Sara Gómez
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Tommaso Giovannini
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Chiara Cappelli
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
9
|
Gómez S, Rojas-Valencia N, Giovannini T, Restrepo A, Cappelli C. Ring Vibrations to Sense Anionic Ibuprofen in Aqueous Solution as Revealed by Resonance Raman. Molecules 2022; 27:molecules27020442. [PMID: 35056755 PMCID: PMC8780161 DOI: 10.3390/molecules27020442] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Revised: 12/31/2021] [Accepted: 01/04/2022] [Indexed: 12/07/2022] Open
Abstract
We unravel the potentialities of resonance Raman spectroscopy to detect ibuprofen in diluted aqueous solutions. In particular, we exploit a fully polarizable quantum mechanics/molecular mechanics (QM/MM) methodology based on fluctuating charges coupled to molecular dynamics (MD) in order to take into account the dynamical aspects of the solvation phenomenon. Our findings, which are discussed in light of a natural bond orbital (NBO) analysis, reveal that a selective enhancement of the Raman signal due to the normal mode associated with the C-C stretching in the ring, νC=C, can be achieved by properly tuning the incident wavelength, thus facilitating the recognition of ibuprofen in water samples.
Collapse
Affiliation(s)
- Sara Gómez
- Classe di Scienze, Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy;
- Correspondence: (S.G.); (C.C.)
| | - Natalia Rojas-Valencia
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellin 050010, Colombia; (N.R.-V.); (A.R.)
| | - Tommaso Giovannini
- Classe di Scienze, Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy;
| | - Albeiro Restrepo
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellin 050010, Colombia; (N.R.-V.); (A.R.)
| | - Chiara Cappelli
- Classe di Scienze, Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy;
- Correspondence: (S.G.); (C.C.)
| |
Collapse
|