1
|
Bowles J, Jähnigen S, Agostini F, Vuilleumier R, Zehnacker A, Calvo F, Clavaguéra C. Vibrational Circular Dichroism Spectroscopy with a Classical Polarizable Force Field: Alanine in the Gas and Condensed Phases. Chemphyschem 2024; 25:e202300982. [PMID: 38318765 DOI: 10.1002/cphc.202300982] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2023] [Revised: 02/05/2024] [Accepted: 02/06/2024] [Indexed: 02/07/2024]
Abstract
Polarizable force fields are an essential component for the chemically accurate modeling of complex molecular systems with a significant degree of fluxionality, beyond harmonic or perturbative approximations. In this contribution we examine the performance of such an approach for the vibrational spectroscopy of the alanine amino acid, in the gas and condensed phases, from the Fourier transform of appropriate time correlation functions generated along molecular dynamics (MD) trajectories. While the infrared (IR) spectrum only requires the electric dipole moment, the vibrational circular dichroism (VCD) spectrum further requires knowledge of the magnetic dipole moment, for which we provide relevant expressions to be used with polarizable force fields. The AMOEBA force field was employed here to model alanine in the neutral and zwitterionic isolated forms, solvated by water or nitrogen, and as a crystal. Within this framework, comparison of the electric and magnetic dipole moments to those obtained with nuclear velocity perturbation theory based on density-functional theory for the same MD trajectories are found to agree well with one another. The statistical convergence of the IR and VCD spectra is examined and found to be more demanding in the latter case. Comparisons with experimental frequencies are also provided for the condensed phases.
Collapse
Affiliation(s)
- Jessica Bowles
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405, Orsay, France
| | - Sascha Jähnigen
- PASTEUR Laboratory, Département de Chimie, Ecole Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005, Paris, France
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405, Orsay, France
| | - Rodolphe Vuilleumier
- PASTEUR Laboratory, Département de Chimie, Ecole Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005, Paris, France
| | - Anne Zehnacker
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay UMR8214, 91405, Orsay, France
| | - Florent Calvo
- Université Grenoble Alpes, CNRS, LIPhy, 38000, Grenoble, France
| | - Carine Clavaguéra
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405, Orsay, France
| |
Collapse
|
2
|
Gómez S, Ambrosetti M, Giovannini T, Cappelli C. Close-Up Look at Electronic Spectroscopic Signatures of Common Pharmaceuticals in Solution. J Phys Chem B 2024; 128:2432-2446. [PMID: 38416564 DOI: 10.1021/acs.jpcb.3c07795] [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: 02/29/2024]
Abstract
Simulating electronic properties and spectral signals requires robust computational approaches that need tuning with the system's peculiarities. In this paper, we test implicit and fully atomistic solvation models for the calculation of UV-vis and electronic circular dichroism (ECD) spectra of two pharmaceutically relevant molecules, namely, (2S)-captopril and (S)-naproxen, dissolved in aqueous solution. Room temperature molecular dynamics simulations reveal that these two drugs establish strong contacts with the surrounding solvent molecules via hydrogen bonds. Such specific interactions, which play a major role in the spectral response and are neglected in implicit approaches, are further characterized and quantified with natural bond orbital methods. Our calculations show that simulated spectra, and especially ECD, are in good agreement with experiments solely when conformational and configurational dynamics, mutual polarization, and solute-solvent repulsion effects are considered.
Collapse
Affiliation(s)
- Sara Gómez
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, Pisa 56126, Italy
| | - Matteo Ambrosetti
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, Pisa 56126, Italy
| | - Tommaso Giovannini
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, Pisa 56126, Italy
| | - Chiara Cappelli
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, Pisa 56126, Italy
| |
Collapse
|
3
|
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
|
4
|
Merten C. Modelling solute-solvent interactions in VCD spectra analysis with the micro-solvation approach. Phys Chem Chem Phys 2023; 25:29404-29414. [PMID: 37881890 DOI: 10.1039/d3cp03408a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2023]
Abstract
Vibrational circular dichroism (VCD) spectroscopy has become an important part of the (stereo-)chemists' toolbox as a reliable method for the determination of absolute configurations. Being the chiroptical version of infrared spectroscopy, it has also been recognized as being very sensitive to conformational changes and intermolecular interactions. This sensitivity originates from the fact that the VCD spectra of individual conformers are often more different than their IR spectra, so that changes in conformational distributions or band positions and intensities become more pronounced. What is an advantage for studies focussing on intermolecular interactions can, however, quickly turn into a major obstacle during AC determinations: solute-solvent interactions can have a strong influence on spectral signatures and they must be accurately treated when simulating VCD and IR spectra. In this perspective, we showcase selected examples which exhibit particularly pronounced solvent effects. It is demonstrated that it is typically sufficient to model solute-solvent interactions by placing single solvent molecules near hydrogen bonding sites of the solute and subsequently use the optimized structures for spectra simulations. This micro-solvation approach works reasonably well for medium-sized, not too conformationally flexible molecules. We thus also discuss its limitations and outline the next steps that method development needs to take in order to further improve the workflows for VCD spectra predictions.
Collapse
Affiliation(s)
- Christian Merten
- Ruhr Universität Bochum, Fakultät für Chemie und Biochemie, Organische Chemie II, Universitätsstraße 150, 44801 Bochum, Germany.
| |
Collapse
|
5
|
Eikås KDR, Beerepoot MTP, Ruud K. A Computational Protocol for Vibrational Circular Dichroism Spectra of Cyclic Oligopeptides. J Phys Chem A 2022; 126:5458-5471. [PMID: 35930395 PMCID: PMC9393892 DOI: 10.1021/acs.jpca.2c02953] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Cyclic peptides are a promising class of compounds for next-generation antibiotics as they may provide new ways of limiting antibiotic resistance development. Although their cyclic structure will introduce some rigidity, their conformational space is large and they usually have multiple chiral centers that give rise to a wide range of possible stereoisomers. Chiroptical spectroscopies such as vibrational circular dichroism (VCD) are used to assign stereochemistry and discriminate enantiomers of chiral molecules, often in combination with electronic structure methods. The reliable determination of the absolute configuration of cyclic peptides will require robust computational methods than can identify all significant conformers and their relative population and reliably assign their stereochemistry from their chiroptical spectra by comparison with ab initio calculated spectra. We here present a computational protocol for the accurate calculation of the VCD spectra of a series of flexible cyclic oligopeptides. The protocol builds on the Conformer-Rotamer Ensemble Sampling Tool (CREST) developed by Grimme and co-workers ( Phys. Chem. Chem. Phys. 2020, 22, 7169-7192 and J. Chem. Theory. Comput. 2019, 15, 2847-2862) in combination with postoptimizations using B3LYP and moderately sized basis sets. Our recommended computational protocol for the computation of VCD spectra of cyclic oligopeptides consists of three steps: (1) conformational sampling with CREST and tight-binding density functional theory (xTB); (2) energy ranking based on single-point energy calculations as well as geometry optimization and VCD calculations of conformers that are within 2.5 kcal/mol of the most stable conformer using B3LYP/6-31+G*/CPCM; and (3) VCD spectra generation based on Boltzmann weighting with Gibbs free energies. Our protocol provides a feasible basis for generating VCD spectra also for larger cyclic peptides of biological/pharmaceutical interest and can thus be used to investigate promising compounds for next-generation antibiotics.
Collapse
Affiliation(s)
- Karolina Di Remigio Eikås
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, 9037 Tromsø, Norway
| | - Maarten T P Beerepoot
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, 9037 Tromsø, Norway
| | - Kenneth Ruud
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, 9037 Tromsø, Norway.,Norwegian Defence Research Establishment, P.O. Box 25, 2027 Kjeller, Norway
| |
Collapse
|
6
|
|
7
|
Ambrosetti M, Skoko S, Giovannini T, Cappelli C. Quantum Mechanics/Fluctuating Charge Protocol to Compute Solvatochromic Shifts. J Chem Theory Comput 2021; 17:7146-7156. [PMID: 34619965 PMCID: PMC8582258 DOI: 10.1021/acs.jctc.1c00763] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
![]()
Despite the potentialities
of the quantum mechanics (QM)/fluctuating
charge (FQ) approach to model the spectral properties of solvated
systems, its extensive use has been hampered by the lack of reliable
parametrizations of solvents other than water. In this paper, we substantially
extend the applicability of QM/FQ to solvating environments of different
polarities and hydrogen-bonding capabilities. The reliability and
robustness of the approach are demonstrated by challenging the model
to simulate solvatochromic shifts of four organic chromophores, which
display large shifts when dissolved in apolar, aprotic or polar, protic
solvents.
Collapse
Affiliation(s)
| | - Sulejman Skoko
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
8
|
Vermeyen T, Brence J, Van Echelpoel R, Aerts R, Acke G, Bultinck P, Herrebout W. Exploring machine learning methods for absolute configuration determination with vibrational circular dichroism. Phys Chem Chem Phys 2021; 23:19781-19789. [PMID: 34524304 DOI: 10.1039/d1cp02428k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The added value of supervised Machine Learning (ML) methods to determine the Absolute Configuration (AC) of compounds from their Vibrational Circular Dichroism (VCD) spectra was explored. Among all ML methods considered, Random Forest (RF) and Feedforward Neural Network (FNN) yield the best performance for identification of the AC. At its best, FNN allows near-perfect AC determination, with accuracy of prediction up to 0.995, while RF combines good predictive accuracy (up to 0.940) with the ability to identify the spectral areas important for the identification of the AC. No loss in performance of either model is observed as long as the spectral sampling interval used does not exceed the spectral bandwidth. Increasing the sampling interval proves to be the best method to lower the dimensionality of the input data, thereby decreasing the computational cost associated with the training of the models.
Collapse
Affiliation(s)
- Tom Vermeyen
- Department of Chemistry, University of Antwerp, Groenenborgerlaan 171, B-2020 Antwerp, Belgium. .,Department of Chemistry, Ghent University, Krijgslaan 281, B-9000 Ghent, Belgium.
| | - Jure Brence
- Department of Knowledge Technologies, Jožef Stefan Institute, Jamova cesta 39, 1000 Ljubljana, Slovenia.,Jožef Stefan International Postgraduate School, Jamova cesta 39, 1000 Ljubljana, Slovenia
| | - Robin Van Echelpoel
- Department of Chemistry, University of Antwerp, Groenenborgerlaan 171, B-2020 Antwerp, Belgium.
| | - Roy Aerts
- Department of Chemistry, University of Antwerp, Groenenborgerlaan 171, B-2020 Antwerp, Belgium.
| | - Guillaume Acke
- Department of Chemistry, Ghent University, Krijgslaan 281, B-9000 Ghent, Belgium.
| | - Patrick Bultinck
- Department of Chemistry, Ghent University, Krijgslaan 281, B-9000 Ghent, Belgium.
| | - Wouter Herrebout
- Department of Chemistry, University of Antwerp, Groenenborgerlaan 171, B-2020 Antwerp, Belgium.
| |
Collapse
|
9
|
Loco D, Lagardère L, Adjoua O, Piquemal JP. Atomistic Polarizable Embeddings: Energy, Dynamics, Spectroscopy, and Reactivity. Acc Chem Res 2021; 54:2812-2822. [PMID: 33961401 PMCID: PMC8264944 DOI: 10.1021/acs.accounts.0c00662] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Indexed: 12/20/2022]
Abstract
The computational modeling of realistic extended systems, relevant in, e.g., Chemistry and Biophysics, is a fundamental problem of paramount importance in contemporary research. Enzymatic catalysis and photoinduced processes in pigment-protein complexes are typical problems targeted by computer-aided approaches, to complement experiments as interpretative tools at a molecular scale. The daunting complexity of this task lies in between the opposite stringent requirements of results' reliability for structural/dynamical properties and related intermolecular interactions, and a mandatory principle of realism in the modeling strategy. Therefore, in practice, a truly realistic computational model of a biologically relevant system can easily fail to meet the accuracy requirement, in order to balance the excessive computational cost necessary to reach the desired precision.To address such an "accuracy vs reality" dualistic requirement, mixed quantum mechanics/classical mechanics approaches within Atomistic (i.e., preserving the discrete particle configuration) Polarizable Embeddings (QM/APEs) methods have been proposed over the years. In this Account, we review recent developments in the design and application of general QM/APE methods, targeting situations where a local intrinsically quantum behavior is coupled to a large molecular system (i.e., an environment), often involving processes with different dynamical time scales, in order to avoid brute-force, unpractical quantum chemistry calculations on the complete system.In the first place, our interest is devoted to the available APEs models presently implemented in computational software, highlighting the quantum chemistry methods that can be used to treat the QM subsystem. We review the coupling strategy between the QM subsystem and the APE, which requires to examine the way the QM/MM mutual interactions are accounted for and how the polarization of the classical environment is considered with respect to (wrt) the quantum variables. Because of the need of reliable molecular and macromolecular structures, a pivotal aspect to address here is the handling of the system dynamics (i.e., gradients wrt nuclear positions are required), especially for large molecular assemblies composed by an overwhelming number of atoms, exploring many conformations on a complex energy landscape.Alongside, we highlight our views on the necessary steps to take toward more accurate general-purposes and transferable explicit embeddings. The main objective to achieve here is to design a more physically grounded multiscale approach. To do so, one should apply advanced new generation classical models to account for refined induction effects that are able to (i) improve the quality of QM/MM interaction energies; (ii) enhance transferability by avoiding the compulsory partial (or total) reparameterization of the classical model. Moreover, the extension of recent developments originating from the field of advanced classical molecular dynamics (MD) to the realm of QM/APE methods is a key direction to improve both speed and efficiency for the phase space exploration of systems of growing size and complexity.Lastly, we point out specific research topics where an advanced QM/APE dynamics can certainly shed some light. For example, we discuss chemical reactions in "harsh" environments and the case of spectroscopic theoretical modeling where the inclusion of refined environment effects is often mandatory.
Collapse
Affiliation(s)
- Daniele Loco
- Laboratoire
de Chimie Théorique, Sorbonne Université,
UMR 7616 CNRS, 75005 Paris, France
| | - Louis Lagardère
- Laboratoire
de Chimie Théorique, Sorbonne Université,
UMR 7616 CNRS, 75005 Paris, France
- Intitut
Parisien de Chimie Physique et Théorique, Sorbonne Université, FR 2622 CNRS, 75005 Paris, France
| | - Olivier Adjoua
- Laboratoire
de Chimie Théorique, Sorbonne Université,
UMR 7616 CNRS, 75005 Paris, France
| | - Jean-Philip Piquemal
- Laboratoire
de Chimie Théorique, Sorbonne Université,
UMR 7616 CNRS, 75005 Paris, France
- Institut
Universitaire de France, F-75005 Paris, France
- Department
of Biomedical Engineering, The University
of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
10
|
Paolino M, Giovannini T, Manathunga M, Latterini L, Zampini G, Pierron R, Léonard J, Fusi S, Giorgi G, Giuliani G, Cappelli A, Cappelli C, Olivucci M. On the Transition from a Biomimetic Molecular Switch to a Rotary Molecular Motor. J Phys Chem Lett 2021; 12:3875-3884. [PMID: 33856801 DOI: 10.1021/acs.jpclett.1c00526] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The experimental investigation of the unidirectional motion characterizing the photoisomerization of single-molecule rotary motors requires accessible lab prototypes featuring an electronic circular dichroism (ECD) signal that is sensitive to the geometrical and electronic changes occurring during an ultrafast reactive process. Here we report a combined experimental/computational study of a candidate obtained via the asymmetrization of a light-driven biomimetic molecular switch. We show that the achieved motor has an ECD band that is remarkably sensitive to the isomerization motion, and it is therefore suitable for time-resolved ECD studies. However, we also find that, unexpectedly, the synthesized motor isomerizes on a time scale longer than the subpicosecond time measured for the achiral parent, a result that points to alternative candidates conserving a high reaction speed.
Collapse
Affiliation(s)
- Marco Paolino
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, Via A. Moro 2, 53100 Siena, Italy
| | | | - Madushanka Manathunga
- Chemistry Department, Bowling Green State University, Bowling Green, Ohio 43403-0001, United States
| | - Loredana Latterini
- Dipartimento di Chimica, Biologia e Biotecnologie, Università di Perugia, Via Elce di Sotto, 8, 06123 Perugia, Italy
| | - Giulia Zampini
- Dipartimento di Chimica, Biologia e Biotecnologie, Università di Perugia, Via Elce di Sotto, 8, 06123 Perugia, Italy
| | - Robin Pierron
- Université de Strasbourg, CNRS, Institut de Physique et Chimie des Matériaux de Strasbourg, UMR 7504, F-67000 Strasbourg, France
| | - Jérémie Léonard
- Université de Strasbourg, CNRS, Institut de Physique et Chimie des Matériaux de Strasbourg, UMR 7504, F-67000 Strasbourg, France
| | - Stefania Fusi
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Gianluca Giorgi
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Germano Giuliani
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Andrea Cappelli
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Massimo Olivucci
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, Via A. Moro 2, 53100 Siena, Italy
- Chemistry Department, Bowling Green State University, Bowling Green, Ohio 43403-0001, United States
| |
Collapse
|
11
|
Polavarapu PL, Santoro E, Covington CL, Johnson JL, Puente AR, Schley ND, Kallingathodi Z, Prakasan PC, Haleema S, Thomas AA, Ibnusaud I. How important are the intermolecular hydrogen bonding interactions in methanol solvent for interpreting the chiroptical properties? SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2021; 247:119094. [PMID: 33142265 DOI: 10.1016/j.saa.2020.119094] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/07/2020] [Revised: 10/10/2020] [Accepted: 10/12/2020] [Indexed: 06/11/2023]
Abstract
Two crispine A analogs and tetrahydrofuro[2,3-b]furan-3,3a(6aH)-diol, endowed with hydroxyl groups that can participate in intramolecular hydrogen bonding, have been synthesized and experimental vibrational circular dichroism (VCD) spectra and optical rotatory dispersion (ORD) data have been measured in CD3OD/CH3OH solvents. The absolute configurations (ACs) of these compounds have been determined using their synthetic schemes, supplemented wherever possible with X-ray diffraction data. The ACs are also analyzed with quantum chemical (QC) calculations of VCD and ORD utilizing implicit solvation as well as explicit solvation models, with the later employing classical molecular dynamics (MD) simulations. It is found that VCD calculations with implicit solvation model are adequate for determining the ACs, despite propensity of studied compounds for intermolecular hydrogen bonding between solute and solvent molecules. This observation is important because time-consuming MD simulations may not be necessary in the type of situations studied here. Additionally, it is found that the QC predicted VCD spectra provided enough diastereomer discrimination for determining the correct AC of studied compounds independently. The same observation did not apply to ORD.
Collapse
Affiliation(s)
| | - Ernesto Santoro
- Department of Chemistry, Vanderbilt University, Nashville, TN 37235, USA
| | - Cody L Covington
- Department of Chemistry, Vanderbilt University, Nashville, TN 37235, USA; Department of Chemistry, Austin Peay State University, Clarksville, TN 37044, USA
| | - Jordan L Johnson
- Department of Chemistry, Vanderbilt University, Nashville, TN 37235, USA
| | - Andrew R Puente
- Department of Chemistry, Vanderbilt University, Nashville, TN 37235, USA
| | - Nathan D Schley
- Department of Chemistry, Vanderbilt University, Nashville, TN 37235, USA.
| | - Zabeera Kallingathodi
- Institute for Integrated Programmes and Research in Basic Sciences, Mahatma Gandhi University, Kottayam 686560, India
| | - Prasanth C Prakasan
- Institute for Integrated Programmes and Research in Basic Sciences, Mahatma Gandhi University, Kottayam 686560, India
| | - Simimole Haleema
- Institute for Integrated Programmes and Research in Basic Sciences, Mahatma Gandhi University, Kottayam 686560, India
| | - Annu Anna Thomas
- Institute for Integrated Programmes and Research in Basic Sciences, Mahatma Gandhi University, Kottayam 686560, India
| | - Ibrahim Ibnusaud
- Institute for Integrated Programmes and Research in Basic Sciences, Mahatma Gandhi University, Kottayam 686560, India.
| |
Collapse
|
12
|
Marrazzini G, Giovannini T, Scavino M, Egidi F, Cappelli C, Koch H. Multilevel Density Functional Theory. J Chem Theory Comput 2021; 17:791-803. [PMID: 33449681 PMCID: PMC7880574 DOI: 10.1021/acs.jctc.0c00940] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
![]()
Following recent
developments in multilevel embedding methods,
we introduce a novel density matrix-based multilevel approach within
the framework of density functional theory (DFT). In this multilevel
DFT, the system is partitioned in an active and an inactive fragment,
and all interactions are retained between the two parts. The decomposition
of the total system is performed upon the density matrix. The orthogonality
between the two parts is maintained by solving the Kohn–Sham
equations in the MO basis for the active part only, while keeping
the inactive density matrix frozen. This results in the reduction
of computational cost. We outline the theory and implementation and
discuss the differences and similarities with state-of-the-art DFT
embedding methods. We present applications to aqueous solutions of
methyloxirane and glycidol.
Collapse
Affiliation(s)
- Gioia Marrazzini
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Tommaso Giovannini
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - Marco Scavino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Franco Egidi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Henrik Koch
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
13
|
Kirchner B, Blasius J, Esser L, Reckien W. Predicting Vibrational Spectroscopy for Flexible Molecules and Molecules with Non‐Idle Environments. ADVANCED THEORY AND SIMULATIONS 2020. [DOI: 10.1002/adts.202000223] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Barbara Kirchner
- Mulliken Center for Theoretical Chemistry Rheinische Friedrich‐Wilhelms‐Universität Bonn Beringstr. 4+6 D‐53115 Bonn Germany
| | - Jan Blasius
- Mulliken Center for Theoretical Chemistry Rheinische Friedrich‐Wilhelms‐Universität Bonn Beringstr. 4+6 D‐53115 Bonn Germany
| | - Lars Esser
- Mulliken Center for Theoretical Chemistry Rheinische Friedrich‐Wilhelms‐Universität Bonn Beringstr. 4+6 D‐53115 Bonn Germany
| | - Werner Reckien
- Mulliken Center for Theoretical Chemistry Rheinische Friedrich‐Wilhelms‐Universität Bonn Beringstr. 4+6 D‐53115 Bonn Germany
| |
Collapse
|
14
|
Skoko S, Ambrosetti M, Giovannini T, Cappelli C. Simulating Absorption Spectra of Flavonoids in Aqueous Solution: A Polarizable QM/MM Study. Molecules 2020; 25:E5853. [PMID: 33322361 PMCID: PMC7764712 DOI: 10.3390/molecules25245853] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2020] [Revised: 12/04/2020] [Accepted: 12/07/2020] [Indexed: 11/16/2022] Open
Abstract
We present a detailed computational study of the UV/Vis spectra of four relevant flavonoids in aqueous solution, namely luteolin, kaempferol, quercetin, and myricetin. The absorption spectra are simulated by exploiting a fully polarizable quantum mechanical (QM)/molecular mechanics (MM) model, based on the fluctuating charge (FQ) force field. Such a model is coupled with configurational sampling obtained by performing classical molecular dynamics (MD) simulations. The calculated QM/FQ spectra are compared with the experiments. We show that an accurate reproduction of the UV/Vis spectra of the selected flavonoids can be obtained by appropriately taking into account the role of configurational sampling, polarization, and hydrogen bonding interactions.
Collapse
Affiliation(s)
- Sulejman Skoko
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy; (S.S.); (M.A.)
| | - Matteo Ambrosetti
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy; (S.S.); (M.A.)
| | - Tommaso Giovannini
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), 7491 Trondheim, Norway;
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy; (S.S.); (M.A.)
| |
Collapse
|
15
|
Marrazzini G, Giovannini T, Egidi F, Cappelli C. Calculation of Linear and Non-linear Electric Response Properties of Systems in Aqueous Solution: A Polarizable Quantum/Classical Approach with Quantum Repulsion Effects. J Chem Theory Comput 2020; 16:6993-7004. [PMID: 33058671 PMCID: PMC8015238 DOI: 10.1021/acs.jctc.0c00674] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2020] [Indexed: 01/10/2023]
Abstract
We present a computational study of polarizabilities and hyperpolarizabilities of organic molecules in aqueous solutions, focusing on solute-water interactions and the way they affect a molecule's linear and non-linear electric response properties. We employ a polarizable quantum mechanics/molecular mechanics (QM/MM) computational model that treats the solute at the QM level while the solvent is treated classically using a force field that includes polarizable charges and dipoles, which dynamically respond to the solute's quantum-mechanical electron density. Quantum confinement effects are also treated by means of a recently implemented method that endows solvent molecules with a parametric electron density, which exerts Pauli repulsion forces upon the solute. By applying the method to a set of aromatic molecules in solution we show that, for both polarizabilities and first hyperpolarizabilities, observed solution values are the result of a delicate balance between electrostatics, hydrogen-bonding, and non-electrostatic solute solvent interactions.
Collapse
Affiliation(s)
- Gioia Marrazzini
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
| | - Tommaso Giovannini
- Department
of Chemistry, Norwegian University of Science
and Technology, Trondheim 7491, Norway
| | - Franco Egidi
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
| | - Chiara Cappelli
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
| |
Collapse
|
16
|
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
|
17
|
Giovannini T, Egidi F, Cappelli C. Molecular spectroscopy of aqueous solutions: a theoretical perspective. Chem Soc Rev 2020; 49:5664-5677. [PMID: 32744278 DOI: 10.1039/c9cs00464e] [Citation(s) in RCA: 50] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/28/2024]
Abstract
Computational spectroscopy is an invaluable tool to both accurately reproduce the spectra of molecular systems and provide a rationalization for the underlying physics. However, the inherent difficulty to accurately model systems in aqueous solutions, owing to water's high polarity and ability to form hydrogen bonds, has severely hampered the development of the field. In this tutorial review we present a technique developed and tested in recent years based on a fully atomistic and polarizable classical modeling of water coupled with a quantum mechanical description of the solute. Thanks to its unparalleled accuracy and versatility, this method can change the perspective of computational and experimental chemists alike.
Collapse
Affiliation(s)
| | - Franco Egidi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy.
| |
Collapse
|
18
|
Gómez S, Giovannini T, Cappelli C. Absorption spectra of xanthines in aqueous solution: a computational study. Phys Chem Chem Phys 2020; 22:5929-5941. [PMID: 32115599 DOI: 10.1039/c9cp05420k] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We present a detailed computational analysis of the UV/Vis spectra of caffeine, paraxanthine and theophylline in aqueous solution. A hierarchy of solvation approaches for modeling the aqueous environment have been tested, ranging from the continuum model to the non-polarizable and polarizable quantum mechanical (QM)/molecular mechanics (MM) models, with and without the explicit inclusion of water molecules in the QM portion. The computed results are directly compared with the experimental data, thus highlighting the role of electrostatic, polarization and hydrogen boding solute-solvent interactions.
Collapse
Affiliation(s)
- Sara Gómez
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy.
| | | | | |
Collapse
|
19
|
Krupová M, Kessler J, Bouř P. Recent Trends in Chiroptical Spectroscopy: Theory and Applications of Vibrational Circular Dichroism and Raman Optical Activity. Chempluschem 2020; 85:561-575. [DOI: 10.1002/cplu.202000014] [Citation(s) in RCA: 39] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2020] [Revised: 02/18/2020] [Indexed: 12/13/2022]
Affiliation(s)
- Monika Krupová
- Institute of Organic Chemistry and Biochemistry Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
- Faculty of Mathematics and PhysicsCharles University Ke Karlovu 3 12116 Prague 2 Czech Republic
| | - Jiří Kessler
- Institute of Organic Chemistry and Biochemistry Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| | - Petr Bouř
- Institute of Organic Chemistry and Biochemistry Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| |
Collapse
|
20
|
Polavarapu PL, Santoro E. Vibrational optical activity for structural characterization of natural products. Nat Prod Rep 2020; 37:1661-1699. [DOI: 10.1039/d0np00025f] [Citation(s) in RCA: 44] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
This review presents the recent progress towards elucidating the structures of chiral natural products and applications using vibrational optical activity (VOA) spectroscopy.
Collapse
|
21
|
Weirich L, Magalhães de Oliveira J, Merten C. How many solvent molecules are required to solvate chiral 1,2-diols with hydrogen bonding solvents? A VCD spectroscopic study. Phys Chem Chem Phys 2020; 22:1525-1533. [DOI: 10.1039/c9cp06030h] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A VCD spectroscopic analysis of selected model systems for solute–solvent interactions of chiral diols with hydrogen bonding solvents DMSO and ACN.
Collapse
Affiliation(s)
- Luisa Weirich
- Ruhr-Universität Bochum
- Fakultät für Chemie und Biochemie
- Organische Chemie II
- 44801 Bochum
- Germany
| | | | - Christian Merten
- Ruhr-Universität Bochum
- Fakultät für Chemie und Biochemie
- Organische Chemie II
- 44801 Bochum
- Germany
| |
Collapse
|
22
|
Le Barbu-Debus K, Bowles J, Jähnigen S, Clavaguéra C, Calvo F, Vuilleumier R, Zehnacker A. Assessing cluster models of solvation for the description of vibrational circular dichroism spectra: synergy between static and dynamic approaches. Phys Chem Chem Phys 2020; 22:26047-26068. [DOI: 10.1039/d0cp03869e] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Solvation effects are essential for defining the shape of vibrational circular dichroism (VCD) spectra.
Collapse
Affiliation(s)
- Katia Le Barbu-Debus
- Institut des Sciences Moléculaires d’Orsay (ISMO)
- CNRS
- Université Paris-Saclay
- F-91405 Orsay
- France
| | - Jessica Bowles
- Université Paris-Saclay
- CNRS
- Institut de Chimie Physique
- UMR8000
- 91405 Orsay
| | - Sascha Jähnigen
- PASTEUR
- Département de Chimie
- Ecole Normale Supérieure
- PSL University
- Sorbonne Université
| | - Carine Clavaguéra
- Université Paris-Saclay
- CNRS
- Institut de Chimie Physique
- UMR8000
- 91405 Orsay
| | - Florent Calvo
- Université Grenoble Alpes
- CNRS
- LiPhy
- F-38000 Grenoble
- France
| | - Rodolphe Vuilleumier
- PASTEUR
- Département de Chimie
- Ecole Normale Supérieure
- PSL University
- Sorbonne Université
| | - Anne Zehnacker
- Institut des Sciences Moléculaires d’Orsay (ISMO)
- CNRS
- Université Paris-Saclay
- F-91405 Orsay
- France
| |
Collapse
|
23
|
Giovannini T, Riso RR, Ambrosetti M, Puglisi A, Cappelli C. Electronic transitions for a fully polarizable QM/MM approach based on fluctuating charges and fluctuating dipoles: Linear and corrected linear response regimes. J Chem Phys 2019; 151:174104. [PMID: 31703497 DOI: 10.1063/1.5121396] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
The fully polarizable Quantum Mechanics/Molecular Mechanics (QM/MM) approach based on fluctuating charges and fluctuating dipoles, named QM/FQFμ [T. Giovannini et al., J. Chem. Theory Comput. 15, 2233 (2019)], is extended to the calculation of vertical excitation energies of solvated molecular systems. Excitation energies are defined within two different solvation regimes, i.e., linear response (LR), where the response of the MM portion is adjusted to the QM transition density, and corrected-Linear Response (cLR) in which the MM response is adjusted to the relaxed QM density, thus being able to account for charge equilibration in the excited state. The model, which is specified in terms of three physical parameters (electronegativity, chemical hardness, and polarizability) is applied to vacuo-to-water solvatochromic shifts of aqueous solutions of para-nitroaniline, pyridine, and pyrimidine. The results show a good agreement with their experimental counterparts, thus highlighting the potentialities of this approach.
Collapse
Affiliation(s)
- Tommaso Giovannini
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | | | | | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
24
|
Giovannini T, Grazioli L, Ambrosetti M, Cappelli C. Calculation of IR Spectra with a Fully Polarizable QM/MM Approach Based on Fluctuating Charges and Fluctuating Dipoles. J Chem Theory Comput 2019; 15:5495-5507. [DOI: 10.1021/acs.jctc.9b00574] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Tommaso Giovannini
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim 7491, Norway
| | - Laura Grazioli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
| | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa 56126, Italy
| |
Collapse
|
25
|
Schwinn K, Ferré N, Huix-Rotllant M. Analytic QM/MM atomic charge derivatives avoiding the scaling of coupled perturbed equations with the MM subsystem size. J Chem Phys 2019; 151:041102. [DOI: 10.1063/1.5115125] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
|
26
|
Di Remigio R, Giovannini T, Ambrosetti M, Cappelli C, Frediani L. Fully Polarizable QM/Fluctuating Charge Approach to Two-Photon Absorption of Aqueous Solutions. J Chem Theory Comput 2019; 15:4056-4068. [DOI: 10.1021/acs.jctc.9b00305] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Roberto Di Remigio
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø - The Arctic University of Norway, N-9037 Tromsø, Norway
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Tommaso Giovannini
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Luca Frediani
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø - The Arctic University of Norway, N-9037 Tromsø, Norway
| |
Collapse
|
27
|
Giovannini T, Lafiosca P, Chandramouli B, Barone V, Cappelli C. Effective yet reliable computation of hyperfine coupling constants in solution by a QM/MM approach: Interplay between electrostatics and non-electrostatic effects. J Chem Phys 2019; 150:124102. [DOI: 10.1063/1.5080810] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Piero Lafiosca
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Balasubramanian Chandramouli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
- Compunet, Istituto Italiano di Tecnologia (IIT), Via Morego 30, 16163 Genova, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
28
|
Giovannini T, Puglisi A, Ambrosetti M, Cappelli C. Polarizable QM/MM Approach with Fluctuating Charges and Fluctuating Dipoles: The QM/FQFμ Model. J Chem Theory Comput 2019; 15:2233-2245. [DOI: 10.1021/acs.jctc.8b01149] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
| | | | | | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
29
|
Sidler D, Bleiziffer P, Riniker S. Beyond the Rosenfeld Equation: Computation of Vibrational Circular Dichroism Spectra for Anisotropic Solutions. J Chem Theory Comput 2019; 15:2492-2503. [PMID: 30802403 DOI: 10.1021/acs.jctc.8b01156] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The difference in absorption of left and right circularly polarized light by chiral molecules can be described by the Rosenfeld equation for isotropic samples. It allows the assignment of the absolute stereochemistry by comparing experimental and computationally derived spectra. Despite the simple form of the Rosenfeld equation, its evaluation in the infrared regime remained challenging, as the contribution from the magnetic dipole operator is zero within the Born-Oppenheimer (BO) approximation. In order to resolve this issue, "beyond BO" theories had to be developed, among which Stephen's magnetic field perturbation (MFP) approach offers a computationally easily accessible form. In this work, optical activity is discussed for cylindrically symmetric solutions, which cannot be described anymore by Rosenfeld's equation due to broken spherical symmetry. Mathematical properties of natural and electric-field induced anisotropies are discussed on the basis of the gauge-independent theoretical framework of Buckingham and Dunn. The issue of achiral noise arising from external field perturbations is considered, and potential remedies are introduced. Natural anisotropic vibrational circular dichroism (VCD) equations are solved numerically by applying the MFP approach within the Hartree-Fock (HF) formalism. Properties of anisotropic VCD spectra are discussed for R-(+)-methyloxirane and (1 S,2 S)-cyclopropane-1,2-dicarbonitrile. In particular, by using a group theoretical argument, a gauge-independent lower bound for the quadrupole contribution of C2-symmetric molecules can be identified, which allows the importance of additional quadrupole terms in anisotropic VCD spectra calculation to be assessed.
Collapse
Affiliation(s)
- Dominik Sidler
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| | - Patrick Bleiziffer
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| | - Sereina Riniker
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| |
Collapse
|
30
|
Egidi F, Giovannini T, Del Frate G, Lemler PM, Vaccaro PH, Cappelli C. A combined experimental and theoretical study of optical rotatory dispersion for (R)-glycidyl methyl ether in aqueous solution. Phys Chem Chem Phys 2019; 21:3644-3655. [PMID: 30383044 DOI: 10.1039/c8cp04445g] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The dispersive optical activity for aqueous solutions of non-rigid (R)-glycidyl methyl ether (R-GME) has been explored synergistically from experimental and theoretical perspectives. Density functional theory analyses performed with the polarizable continuum model for implicit solvation identified nine low-lying stable conformers that are interconverted by rotation about two large-amplitude torsional coordinates. The antagonistic chiroptical signatures predicted for these structural isomers were averaged under a Boltzmann-weighting ansatz to estimate the behavior expected for a thermally equilibrated ensemble. This led to optical rotatory dispersion profiles that reproduced the overall shape of observations but failed to achieve uniform agreement with measured specific-rotation values even when anharmonic vibrational corrections were applied. A mixed QM/FQ paradigm, whereby quantum-mechanical (QM) calculations of optical activity were combined with classical molecular dynamics simulations of explicit solvation that included mutual-polarization effects by means of fluctuating charges (FQ), was enlisted to elucidate the microsolvation environment and gauge its impact upon conformer distributions and response properties. Although quantitative accord with experiments remained elusive, this approach revealed strong variations in the magnitude and sign of rotatory powers for R-GME as the configuration of surrounding water molecules evolved, thereby highlighting the inherently dynamical nature of the solvated chiroptical response, calling into question the validity of "static" descriptions based on the presumption of distinct energy minima, and giving insight into the inherent complexity posed by the modeling of such properties for solvated systems.
Collapse
Affiliation(s)
- Franco Egidi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa, Italy
| | | | | | | | | | | |
Collapse
|
31
|
Shahriar I, Islam MKB, Iqfath M, Rahman A, Halim MA. Solvent effect on vibrational circular dichroism of chiral amino acids. Theor Chem Acc 2019. [DOI: 10.1007/s00214-019-2419-0] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
32
|
Puglisi A, Giovannini T, Antonov L, Cappelli C. Interplay between conformational and solvent effects in UV-visible absorption spectra: curcumin tautomers as a case study. Phys Chem Chem Phys 2019; 21:15504-15514. [PMID: 31259324 DOI: 10.1039/c9cp00907h] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
We present a combined theoretical and experimental study on the UV-vis spectra of enol-keto (EK) and keto-keto (KK) tautomeric forms of curcumin dissolved in aqueous solution. Solvent effects have been investigated by resorting to the implicit polarizable continuum model (QM/PCM) and non-polarizable and fully polarizable QM/MM approaches, the latter based on the fluctuating charges (FQ) force-field. In particular, all methods are challenged to rationalize the contribution of conformational, electrostatic and polarization effects in the calculation of the vertical excitation spectra of curcumin tautomers. The obtained results highlight that for both tautomers specific solute-solvent hydrogen-bond interactions play a minor role with respect to conformational and electrostatic effects.
Collapse
Affiliation(s)
| | | | - Liudmil Antonov
- Bulgarian Academy of Sciences, Institute of Organic Chemistry with Centre of Phytochemistry, Acad. G. Bonchev str., Bldg. 9, Sofia 1113, Bulgaria
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy.
| |
Collapse
|
33
|
Ghidinelli S, Abbate S, Boiadjiev SE, Lightner DA, Longhi G. l-Stercobilin-HCl and d-Urobilin-HCl. Analysis of Their Chiroptical and Conformational Properties by VCD, ECD, and CPL Experiments and MD and DFT Calculations. J Phys Chem B 2018; 122:12351-12362. [DOI: 10.1021/acs.jpcb.8b07954] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Simone Ghidinelli
- Dipartimento di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123 Brescia, Italy
| | - Sergio Abbate
- Dipartimento di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123 Brescia, Italy
| | - Stefan E. Boiadjiev
- Department of Chemistry and Biochemistry, Medical University-Pleven, 1 St. Kl. Ohridski Str., 5800 Pleven, Bulgaria
| | - David A. Lightner
- Department of Chemistry, University of Nevada, Reno, 1664 North Virginia Street, Reno, 89557-0020 Nevada, United States
| | - Giovanna Longhi
- Dipartimento di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123 Brescia, Italy
| |
Collapse
|
34
|
Jähnigen S, Scherrer A, Vuilleumier R, Sebastiani D. VCD‐Verstärkung durch chirale Packungseffekte in molekularen Kristallen. Angew Chem Int Ed Engl 2018. [DOI: 10.1002/ange.201805671] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Sascha Jähnigen
- Martin-Luther-Universität Halle-Wittenberg Institut für Chemie von-Danckelmann-Platz 4 Halle 06120 Deutschland
| | - Arne Scherrer
- Martin-Luther-Universität Halle-Wittenberg Institut für Chemie von-Danckelmann-Platz 4 Halle 06120 Deutschland
| | - Rodolphe Vuilleumier
- PASTEUR, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS 75005 Paris Frankreich
| | - Daniel Sebastiani
- Martin-Luther-Universität Halle-Wittenberg Institut für Chemie von-Danckelmann-Platz 4 Halle 06120 Deutschland
| |
Collapse
|
35
|
Jähnigen S, Scherrer A, Vuilleumier R, Sebastiani D. Chiral Crystal Packing Induces Enhancement of Vibrational Circular Dichroism. Angew Chem Int Ed Engl 2018; 57:13344-13348. [DOI: 10.1002/anie.201805671] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2018] [Revised: 06/25/2018] [Indexed: 01/15/2023]
Affiliation(s)
- Sascha Jähnigen
- Martin-Luther-Universität Halle-WittenbergInstitut für Chemie von-Danckelmann-Platz 4 Halle 06120 Germany
| | - Arne Scherrer
- Martin-Luther-Universität Halle-WittenbergInstitut für Chemie von-Danckelmann-Platz 4 Halle 06120 Germany
| | - Rodolphe Vuilleumier
- PASTEUR, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS 75005 Paris France
| | - Daniel Sebastiani
- Martin-Luther-Universität Halle-WittenbergInstitut für Chemie von-Danckelmann-Platz 4 Halle 06120 Germany
| |
Collapse
|
36
|
Rode JE, Górecki M, Witkowski S, Frelek J. Solvation of 2-(hydroxymethyl)-2,5,7,8-tetramethyl-chroman-6-ol revealed by circular dichroism: a case of chromane helicity rule breaking. Phys Chem Chem Phys 2018; 20:22525-22536. [PMID: 30140796 DOI: 10.1039/c8cp02491j] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
The primary goal of this work is to clarify why 2-(hydroxymethyl)-2,5,7,8-tetramethyl-chroman-6-ol {(S)-TMChM} deviates from the chromane helicity rule under solvent change. The rule, applicable to determining the absolute configuration of molecules containing the chromane chromophore, binds the sign of the 1Lb Cotton effect (CE) with the helicity of the dihydropyran ring. In case of TMChM, however, this CE exhibits extreme solvent dependence: it is negative in non-coordinating solvents and positive in coordinating ones, irrespective of the helicity of the heterocyclic ring. TD-DFT calculations using PCM and hybrid solvation models were performed to explain origin of this phenomenon. It turned out that the 1Lb CE sign directly depends on the position of the phenolic OH group at carbon atom C6 (OHC6). In the absence of interactions with solvents (as in CCl4 or nC6H14) or when a solvent plays proton donor role (as in CHCl3), the OHC6 lies in the phenyl plane and the 1Lb CE sign follows the P/M helicity rule. In contrast, in proton acceptor solvents, like DMSO, CH3OH or CH3CN, the OHC6 group is deflected from the phenyl plane, and the 1Lb CE sign of individual (S)-TMChM conformers depends on the sector in which the OHC6 is located. Thus, in solution, the 1Lb CE sign is an average over different orientations of the OHC6 group and can be positive (as in DMSO and CH3OH) or negative (as in CH3CN) which means that it does not follow the chromane helicity rule. The impact of OHC6 on the 1Lb CE sign and thus the conclusions for the stereochemistry of chromans are demonstrated here for the first time. Additionally, a comparison of experimental and simulated ECD spectra, supported by VCD data, allowed to determine the geometry of intermolecular clusters formed in different solvents.
Collapse
Affiliation(s)
- Joanna E Rode
- Institute of Organic Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland.
| | | | | | | |
Collapse
|
37
|
On the nature of charge-transfer excitations for molecules in aqueous solution: a polarizable QM/MM study. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2259-3] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
38
|
A polarizable embedding approach to second harmonic generation (SHG) of molecular systems in aqueous solutions. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2247-7] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
39
|
Morzan UN, Alonso de Armiño DJ, Foglia NO, Ramírez F, González Lebrero MC, Scherlis DA, Estrin DA. Spectroscopy in Complex Environments from QM–MM Simulations. Chem Rev 2018; 118:4071-4113. [DOI: 10.1021/acs.chemrev.8b00026] [Citation(s) in RCA: 97] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Uriel N. Morzan
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Diego J. Alonso de Armiño
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Nicolás O. Foglia
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Francisco Ramírez
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Mariano C. González Lebrero
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Damián A. Scherlis
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Darío A. Estrin
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| |
Collapse
|
40
|
Shen J, Magesh S, Chen L, Hu L, He Y. Enantiomeric characterization and structure elucidation of LH601A using vibrational circular dichroism spectroscopy. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2018; 192:312-317. [PMID: 29172127 PMCID: PMC10544735 DOI: 10.1016/j.saa.2017.11.033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/15/2017] [Revised: 10/23/2017] [Accepted: 11/14/2017] [Indexed: 06/07/2023]
Abstract
LH601A is a novel non-reactive chiral molecule inhibiting Keap1-Nrf2 protein-protein interaction. The absolute configuration (AC) was independently determined in this study using vibrational circular dichroism (VCD) spectroscopy. Because of band overlapping and broadening in the IR spectrum, a direct VCD spectrum comparison method is devised without the conventional IR band alignment. Being an unbiased AC inquiry, all possible chiralities are evaluated based on the statistical analysis of VCD similarity, Sv. The AC of three-center stereoisomer LH601A is unambiguously assigned to (S,R,S). A comparative study was also carried out to investigate the structural and energy differences of calculated conformers using the polarized continuum model of dimethyl sulfoxide.
Collapse
Affiliation(s)
- Jian Shen
- Department of Medicinal Chemistry, Ernest Mario School of Pharmacy, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; S&M Systems, P.O.Box 6583, Bridgewater, NJ 08807, USA.
| | - Sadagopan Magesh
- Department of Medicinal Chemistry, Ernest Mario School of Pharmacy, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Lin Chen
- Department of Medicinal Chemistry, Ernest Mario School of Pharmacy, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Longqin Hu
- Department of Medicinal Chemistry, Ernest Mario School of Pharmacy, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; The Cancer Institute of New Jersey, New Brunswick, NJ 08901, USA
| | - Yanan He
- Biotools, Inc., 17546 Beeline Hwy, Jupiter, FL 33458, USA
| |
Collapse
|
41
|
Giovannini T, Del Frate G, Lafiosca P, Cappelli C. Effective computational route towards vibrational optical activity spectra of chiral molecules in aqueous solution. Phys Chem Chem Phys 2018; 20:9181-9197. [DOI: 10.1039/c8cp00487k] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
A polarizable QM/MM approach to accurately compute the Vibrational Optical Activity (VOA) spectra of chiral systems is proposed and applied to aqueous solutions of (l)-methyl lactate and (S)-glycidol.
Collapse
|
42
|
Gatto P, Lipparini F, Stamm B. Computation of forces arising from the polarizable continuum model within the domain-decomposition paradigm. J Chem Phys 2017; 147:224108. [DOI: 10.1063/1.5008329] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Paolo Gatto
- Mathematics Division, Center for Computational Engineering Science, RWTH Aachen University, Aachen, Germany
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Benjamin Stamm
- Mathematics Division, Center for Computational Engineering Science, RWTH Aachen University, Aachen, Germany
| |
Collapse
|
43
|
Sæther S, Kjærgaard T, Koch H, Høyvik IM. Density-Based Multilevel Hartree–Fock Model. J Chem Theory Comput 2017; 13:5282-5290. [DOI: 10.1021/acs.jctc.7b00689] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Sandra Sæther
- Department
of Chemistry, The Norwegian University of Science and Technology, Høgskoleringen 5, 7491 Trondheim, Norway
| | - Thomas Kjærgaard
- qLEAP
Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Henrik Koch
- Department
of Chemistry, The Norwegian University of Science and Technology, Høgskoleringen 5, 7491 Trondheim, Norway
| | - Ida-Marie Høyvik
- Department
of Chemistry, The Norwegian University of Science and Technology, Høgskoleringen 5, 7491 Trondheim, Norway
| |
Collapse
|
44
|
Giovannini T, Lafiosca P, Cappelli C. A General Route to Include Pauli Repulsion and Quantum Dispersion Effects in QM/MM Approaches. J Chem Theory Comput 2017; 13:4854-4870. [PMID: 28898079 DOI: 10.1021/acs.jctc.7b00776] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
A methodology to account for nonelectrostatic interactions in Quantum Mechanical (QM)/Molecular Mechanics (MM) approaches is developed. Formulations for Pauli repulsion and dispersion energy, explicitly depending on the QM density, are derived. Such expressions are based on the definition of an auxiliary density on the MM portion and the Tkatchenko-Scheffler (TS) approach, respectively. The developed method is general enough to be applied to any QM/MM method and partition, provided an accurate tuning of a small number of parameters is obtained. The coupling of the method with both nonpolarizable and fully polarizable QM/fluctuating charge (FQ) approaches is reported and applied. A suitable parametrization for the aqueous solution, so that its most representative features are well reproduced, is outlined. Then, the obtained parametrization and method are applied to calculate the nonelectrostatic (repulsion and dispersion) interaction energy of nicotine in aqueous solution.
Collapse
Affiliation(s)
| | - Piero Lafiosca
- Scuola Normale Superiore , Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Chiara Cappelli
- Scuola Normale Superiore , Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
45
|
Giovannini T, Olszówka M, Egidi F, Cheeseman JR, Scalmani G, Cappelli C. Polarizable Embedding Approach for the Analytical Calculation of Raman and Raman Optical Activity Spectra of Solvated Systems. J Chem Theory Comput 2017; 13:4421-4435. [DOI: 10.1021/acs.jctc.7b00628] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
| | - Marta Olszówka
- Scuola Normale
Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Franco Egidi
- Scuola Normale
Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - James R. Cheeseman
- Gaussian
Inc., 340 Quinnipiac Street Building
40, Wallingford, Connecticut 06492, United States
| | - Giovanni Scalmani
- Gaussian
Inc., 340 Quinnipiac Street Building
40, Wallingford, Connecticut 06492, United States
| | - Chiara Cappelli
- Scuola Normale
Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
46
|
Reiter K, Kühn M, Weigend F. Vibrational circular dichroism spectra for large molecules and molecules with heavy elements. J Chem Phys 2017; 146:054102. [DOI: 10.1063/1.4974897] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
|