1
|
Liasi Z, Jensen L, Mikkelsen KV. A Combined Quantum Mechanics and Molecular Mechanics Approach for Simulating the Optical Properties of DNA-Stabilized Silver Nanoclusters. J Chem Theory Comput 2024; 20:937-945. [PMID: 38164716 DOI: 10.1021/acs.jctc.3c01022] [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: 01/03/2024]
Abstract
DNA-stabilized silver nanoclusters have emerged as an intriguing type of nanomaterial due to their unique optical and electronic properties, with potential applications in areas such as biosensing and imaging. The development of efficient methods for modeling these properties is paramount for furthering the understanding and utilization of these clusters. In this study, a hybrid quantum mechanical and molecular mechanical approach for modeling the optical properties of a DNA-templated silver nanocluster is evaluated. The influence of different parameters, including ligand fragmentation, damping, embedding potential, basis set, and density functional, is investigated. The results demonstrate that the most important parameter is the type of atomic properties used to represent the ligands, with isotropic dipole-dipole polarizabilities outperforming the rest. This underscores the importance of an appropriate representation of the ligands, particularly through the selection of the properties used to represent them. Moreover, the results are compared to experimental data, showing that the applied methodology is reliable and effective for the modeling of DNA-stabilized silver nanoclusters. These findings offer valuable insights that may guide future computational efforts to explore and harness the potential of these novel systems.
Collapse
Affiliation(s)
- Zacharias Liasi
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, Copenhagen Ø 2100, Denmark
| | - Lasse Jensen
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Kurt V Mikkelsen
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, Copenhagen Ø 2100, Denmark
| |
Collapse
|
2
|
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
|
3
|
Jørgensen FK, Reinholdt P, Hedegård ED, Kongsted J. Nuclear Magnetic Shielding Constants with the Polarizable Density Embedding Model. J Chem Theory Comput 2022; 18:7384-7393. [PMID: 36332108 DOI: 10.1021/acs.jctc.2c00829] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
We extend the polarizable density embedding (PDE) model to support the calculation of nuclear magnetic resonance (NMR) shielding constants using gauge-including atomic orbitals (GIAOs) within a density functional theory (DFT) framework. The PDE model divides the total system into fragments, describing some by quantum mechanics (QM) and the others through an embedding model. The PDE model uses anisotropic polarizabilities, inter-fragment two-electron Coulomb integrals, and a non-local repulsion operator to emulate the QM effects. The terms involving Coulomb integrals are straightforwardly extended with GIAOs. In contrast, we consider two approaches to handle the gauge dependency of the non-local operator, employing either simple symmetrization or a gauge transformation. We find the latter approach to be most stable with respect to increasing the basis set size of the QM region. We examine the accuracy of the PDE model for calculating NMR shielding constants on several solutes in a water solution. The performance is compared with the classical polarizable embedding (PE) model in addition to supermolecular reference calculations. Based on these systems, we address the basis set convergence characteristics and the QM region size requirements. Furthermore, we investigate the performance of the PDE model for a system with significant electron spill-out. In many cases, we find that the PDE model outperforms the PE model, especially regarding the accuracy of nuclear shielding constants when using small QM region sizes and in systems with significant electron spill-out.
Collapse
Affiliation(s)
- Frederik Kamper Jørgensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230Odense M, Denmark
| | - Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230Odense M, Denmark
| | - Erik Donovan Hedegård
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230Odense M, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230Odense M, Denmark
| |
Collapse
|
4
|
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
|
5
|
Fu M, Wesolowski TA. The Challenge of Accurate Computation of Two-Photon Absorption Properties of Organic Chromophores in the Condensed Phase. J Chem Theory Comput 2021; 17:3652-3665. [PMID: 33944563 DOI: 10.1021/acs.jctc.1c00204] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Two strategies are applied to evaluate the effect of the environment on the two-photon absorption (TPA) cross sections for two characteristic excited states of C2H4 upon complexation with H2O. The supermolecular strategy provides the reference complexation-induced shifts and uses either the EOM-CCSD or ADC(2) method. The embedding strategy is based on frozen-density-embedding theory (FDET) and uses only fundamental constants. The TPA cross sections from high-level supermolecular calculations are extremely basis-set-sensitive. Literature data and the present study indicate that accuracy of the absolute TPA cross sections below 100 atomic units and their shifts below 10 atomic units remains a challenge. The obtained FDET results show a similar basis-set behavior. For the largest basis set (d-aug-cc-pVQZ), TPA cross sections obtained from these two strategies are in excellent agreement. The complexation-induced shifts have the correct sign of the effect and a small (12-33%) relative error in magnitude. The deviations of the FDET-derived shifts from the reference are of similar magnitude as the reliability threshold of the reference shifts.
Collapse
Affiliation(s)
- Mingxue Fu
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - Tomasz A Wesolowski
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| |
Collapse
|
6
|
Grabarek D, Andruniów T. Removing artifacts in polarizable embedding calculations of one- and two-photon absorption spectra of fluorescent proteins. J Chem Phys 2021; 153:215102. [PMID: 33291919 DOI: 10.1063/5.0023434] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
The multiscale calculations involving excited states may suffer from the electron spill-out (ESO) problem. This seems to be especially the case when the environment of the core region, described with the electronic structure method, is approximated by a polarizable force field. The ESO effect often leads to incorrect physical character of electronic excitations, spreading outside the quantum region, which, in turn, results in erroneous absorption spectra. In this work, we investigate means to remove the artifacts in one-photon absorption (OPA) and two-photon absorption (TPA) spectra of green and yellow fluorescent protein representatives. This includes (i) using different basis sets, (ii) extending the core subsystem beyond the chromophore, (iii) modification of polarization interaction between the core region and its environment, and (iv) including the Pauli repulsion through effective core potentials (ECPs). Our results clearly show that ESO is observed when diffuse functions are used to assemble the multielectron wave function regardless of the exchange-correlation functional used. Furthermore, extending the core region, thus accounting for exchange interactions between the chromophore and its environment, leads to even more spurious excited states. Also, damping the interactions between the core subsystem and the polarizable force field is hardly helpful. In contrast, placing ECPs in the position of sites creating the embedding potential leads to the removal of artificious excited states that presumably should not be observed in the OPA and TPA spectra. We prove that it is a reliable and cost-effective approach for systems where the covalent bond(s) between the core region and its environment must be cut.
Collapse
Affiliation(s)
- Dawid Grabarek
- Advanced Materials Engineering and Modelling Group, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, 50-370 Wroclaw, Poland
| | - Tadeusz Andruniów
- Advanced Materials Engineering and Modelling Group, Wroclaw University of Science and Technology, Wyb. Wyspianskiego 27, 50-370 Wroclaw, Poland
| |
Collapse
|
7
|
Attia AAA, Reinholdt P, Kongsted J. Modeling One‐ and Two‐Photon Excitation of 4′‐(Hydroxymethyl)‐4,5′,8‐trimethylpsoralen in Complex with DNA: Solving Electron Spill‐Out Problems in Polarizable QM/MM Calculations. ADVANCED THEORY AND SIMULATIONS 2021. [DOI: 10.1002/adts.202000294] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Amr A. A. Attia
- Faculty of Chemistry and Chemical Engineering Babeş‐Bolyai University Cluj‐Napoca Romania
| | - Peter Reinholdt
- Department of Physics Chemistry and Pharmacy University of Southern Denmark Campusvej 55 Odense M DK–5230 Denmark
| | - Jacob Kongsted
- Department of Physics Chemistry and Pharmacy University of Southern Denmark Campusvej 55 Odense M DK–5230 Denmark
| |
Collapse
|