1
|
Chen Y, Pios SV, Gelin MF, Chen L. Accelerating Molecular Vibrational Spectra Simulations with a Physically Informed Deep Learning Model. J Chem Theory Comput 2024; 20:4703-4710. [PMID: 38825857 DOI: 10.1021/acs.jctc.4c00173] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2024]
Abstract
In recent years, machine learning (ML) surrogate models have emerged as an indispensable tool to accelerate simulations of physical and chemical processes. However, there is still a lack of ML models that can accurately predict molecular vibrational spectra. Here, we present a highly efficient multitask ML surrogate model termed Vibrational Spectra Neural Network (VSpecNN), to accurately calculate infrared (IR) and Raman spectra based on dipole moments and polarizabilities obtained on-the-fly via ML-enhanced molecular dynamics simulations. The methodology is applied to pyrazine, a prototypical polyatomic chromophore. The VSpecNN-predicted energies are well within the chemical accuracy (1 kcal/mol), and the errors for VSpecNN-predicted forces are only half of those obtained from a popular high-performance ML model. Compared to the ab initio reference, the VSpecNN-predicted frequencies of IR and Raman spectra differ only by less than 5.87 cm-1, and the intensities of IR spectra and the depolarization ratios of Raman spectra are well reproduced. The VSpecNN model developed in this work highlights the importance of constructing highly accurate neural network potentials for predicting molecular vibrational spectra.
Collapse
Affiliation(s)
| | | | - Maxim F Gelin
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | | |
Collapse
|
2
|
Mamatkulov K, Zavatski S, Arynbek Y, Esawii HA, Burko A, Bandarenka H, Arzumanyan G. Conformational analysis of lipid membrane mimetics modified with A β42 peptide by Raman spectroscopy and computer simulations. J Biomol Struct Dyn 2024:1-14. [PMID: 38520152 DOI: 10.1080/07391102.2024.2330706] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2023] [Accepted: 03/08/2024] [Indexed: 03/25/2024]
Abstract
Peptide-lipid interactions play an important role in maintaining the integrity and function of the cell membrane. Even slight changes in these interactions can induce the development of various diseases. Specifically, peptide misfolding and aggregation in the membrane is considered to be one of the triggers of Alzheimer's disease (AD), however its exact mechanism is still unclear. To this end, an increase of amyloid-beta (Aβ) peptide concentration in the human brain is widely accepted to gradually produce cytotoxic Aβ aggregates (plaques). These plaques initiate a sequence of pathogenic events ending up in observable symptoms of dementia. Understanding the mechanism of the Aβ interaction with cells is crucial for early detection and prevention of Alzheimer's disease. Hence, in this work, a comprehensive Raman analysis of the Aβ42 conformational dynamics in water and in liposomes and lipodiscs that mimic the membrane system is presented. The obtained results show that the secondary structure of Aβ42 in liposomes is dominated by the α-helix conformation, which remains stable over time. However, it comes as a surprise to reveal that the lipodisc environment induces the transformation of the Aβ42 secondary structure to a β-turn/random coil. Our Raman spectroscopy findings are supported with molecular dynamics (MD) and density functional theory (DFT) simulations, showing their good agreement.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Kahramon Mamatkulov
- Laboratory of Neutron Physics, Sector of Raman Spectroscopy, Joint Institute for Nuclear Research, Dubna, Russia
| | - Siarhei Zavatski
- Applied Plasmonics Laboratory, Belarusian State University of Informatics and Radioelectronics, Minsk, Belarus
| | - Yersultan Arynbek
- Laboratory of Neutron Physics, Sector of Raman Spectroscopy, Joint Institute for Nuclear Research, Dubna, Russia
- Faculty of Physics and Technology, al-Farabi, Kazakh National University, Almaty, Kazakhstan
| | - Heba A Esawii
- Laboratory of Neutron Physics, Sector of Raman Spectroscopy, Joint Institute for Nuclear Research, Dubna, Russia
- Biophysics Department, Faculty of Science, Cairo University, Egypt
| | - Aliaksandr Burko
- Applied Plasmonics Laboratory, Belarusian State University of Informatics and Radioelectronics, Minsk, Belarus
| | - Hanna Bandarenka
- Applied Plasmonics Laboratory, Belarusian State University of Informatics and Radioelectronics, Minsk, Belarus
| | - Grigory Arzumanyan
- Laboratory of Neutron Physics, Sector of Raman Spectroscopy, Joint Institute for Nuclear Research, Dubna, Russia
| |
Collapse
|
3
|
Sahu N, Khire SS, Gadre SR. Combining fragmentation method and high-performance computing: Geometry optimization and vibrational spectra of proteins. J Chem Phys 2023; 159:044309. [PMID: 37522406 DOI: 10.1063/5.0149572] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Accepted: 07/12/2023] [Indexed: 08/01/2023] Open
Abstract
Exploring the structures and spectral features of proteins with advanced quantum chemical methods is an uphill task. In this work, a fragment-based molecular tailoring approach (MTA) is appraised for the CAM-B3LYP/aug-cc-pVDZ-level geometry optimization and vibrational infrared (IR) spectra calculation of ten real proteins containing up to 407 atoms and 6617 basis functions. The use of MTA and the inherently parallel nature of the fragment calculations enables a rapid and accurate calculation of the IR spectrum. The applicability of MTA to optimize the protein geometry and evaluate its IR spectrum employing a polarizable continuum model with water as a solvent is also showcased. The typical errors in the total energy and IR frequencies computed by MTA vis-à-vis their full calculation (FC) counterparts for the studied protein are 5-10 millihartrees and 5 cm-1, respectively. Moreover, due to the independent execution of the fragments, large-scale parallelization can also be achieved. With increasing size and level of theory, MTA shows an appreciable advantage in computer time as well as memory and disk space requirement over the corresponding FCs. The present study suggests that the geometry optimization and IR computations on the biomolecules containing ∼1000 atoms and/or ∼15 000 basis functions using MTA and HPC facility can be clearly envisioned in the near future.
Collapse
Affiliation(s)
- Nityananda Sahu
- Theoretische Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| | - Subodh S Khire
- RIKEN Center for Computational Science, Kobe 650-0047, Japan
| | - Shridhar R Gadre
- Departments of Scientific Computing, Modelling & Simulation and Chemistry, Savitribai Phule Pune University, Pune 411007, India
| |
Collapse
|
4
|
Yamamoto S, Ishiro S, Kessler J, Bouř P. Intense chiral signal from α-helical poly-L-alanine observed in low-frequency Raman optical activity. Phys Chem Chem Phys 2021; 23:26501-26509. [PMID: 34806737 DOI: 10.1039/d1cp04401j] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Raman optical activity (ROA) spectral features reliably indicate the structure of peptides and proteins, but the signal is often weak. However, we observed significantly enhanced low-frequency bands for α-helical poly-L-alanine (PLA) in solution. The biggest ROA signal at ∼100 cm-1 is about 10 times stronger than higher-frequency bands described previously, which facilitates the detection. The low-frequency bands of PLA were compared to those of α-helical proteins. For PLA, density functional simulations well reproduced the experimental spectra and revealed that about 12 alanine residues within two turns of the α-helix generate the strong ROA band. Averaging based on molecular dynamics (MD) provided an even more realistic spectrum compared to the static model. The low-frequency bands could be largely related to a collective motion of the α-helical backbone, partially modulated by the solvent. Helical and intermolecular vibrational coordinates have been introduced and the helical unwinding modes were assigned to the strongest ROA signal at 101-128 cm-1. Further analysis indicated that the helically arranged amide and methyl groups are important for the strong chiral signal of PLA, while the local chiral centers CαH contribute in a minor way only. The strong low-frequency ROA can thus provide precious information about the motions of the peptide backbone and facilitate future protein studies.
Collapse
Affiliation(s)
- Shigeki Yamamoto
- Department of Chemistry, Graduate School of Science, Osaka University, Osaka 560-0043, Japan.
| | - Shota Ishiro
- Department of Chemistry, Graduate School of Science, Osaka University, Osaka 560-0043, Japan.
| | - 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
|
5
|
Esch BVD, Peters LDM, Sauerland L, Ochsenfeld C. Quantitative Comparison of Experimental and Computed IR-Spectra Extracted from Ab Initio Molecular Dynamics. J Chem Theory Comput 2021; 17:985-995. [PMID: 33512155 DOI: 10.1021/acs.jctc.0c01279] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Experimentally measured infrared spectra are often compared to their computed equivalents. However, the accordance is typically characterized by visual inspection, which is prone to subjective judgment. The primary challenge for a similarity-based analysis is that the artifacts introduced by each approach are very different and, therefore, may require preprocessing steps to determine and correct impeding irregularities. To allow for automated objective assessment, we propose a practical and comprehensive workflow involving scaling factors, a novel baseline correction scheme, and peak smoothing. The resulting spectra can then easily be compared quantitatively using similarity measures, for which we found the Pearson correlation coefficient to be the most suitable. The proposed procedure is then applied to compare the agreement of the experimental infrared spectra from the NIST Chemistry Web book with the calculated spectra using standard harmonic frequency analysis and spectra extracted from ab initio molecular dynamics simulations at different levels of theory. We conclude that the direct, quantitative comparison of calculated and measured IR spectra might become a novel, sophisticated approach to benchmark quantum-chemical methods. In the present benchmark, simulated spectra based on ab initio molecular dynamics show in general better agreement with the experiment than static calculations.
Collapse
Affiliation(s)
- Beatriz von der Esch
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Laurens D M Peters
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Lena Sauerland
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| |
Collapse
|
6
|
Kurochka A, Průša J, Kessler J, Kapitán J, Bouř P. α-Synuclein conformations followed by vibrational optical activity. Simulation and understanding of the spectra. Phys Chem Chem Phys 2021; 23:16635-16645. [DOI: 10.1039/d1cp02574k] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
For α-synuclein novel structural markers were identified in vibrational optical activity spectra and supported by theoretical modeling.
Collapse
Affiliation(s)
- Andrii Kurochka
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences
- 16610 Prague
- Czech Republic
- Department of Analytical Chemistry
| | - Jiří Průša
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences
- 16610 Prague
- Czech Republic
- Department of Analytical Chemistry
| | - Jiří Kessler
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences
- 16610 Prague
- Czech Republic
| | - Josef Kapitán
- Department of Optics
- Palacký University Olomouc
- Olomouc
- Czech Republic
| | - Petr Bouř
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences
- 16610 Prague
- Czech Republic
| |
Collapse
|
7
|
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
|
8
|
Keiderling TA. Structure of Condensed Phase Peptides: Insights from Vibrational Circular Dichroism and Raman Optical Activity Techniques. Chem Rev 2020; 120:3381-3419. [DOI: 10.1021/acs.chemrev.9b00636] [Citation(s) in RCA: 50] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Timothy A. Keiderling
- Department of Chemistry, University of Illinois at Chicago 845 West Taylor Street m/c 111, Chicago, Illinois 60607-7061, United States
| |
Collapse
|
9
|
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
|
10
|
Profant V, Johannessen C, Blanch EW, Bouř P, Baumruk V. Effects of sulfation and the environment on the structure of chondroitin sulfate studied via Raman optical activity. Phys Chem Chem Phys 2019; 21:7367-7377. [DOI: 10.1039/c9cp00472f] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Raman optical activity reflects differences in the secondary structure of chondroitin caused by its sulfation.
Collapse
Affiliation(s)
- Václav Profant
- Institute of Physics
- Faculty of Mathematics and Physics
- Charles University
- 121 16 Prague 2
- Czech Republic
| | | | | | - Petr Bouř
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences
- 166 10 Prague 6
- Czech Republic
| | - Vladimír Baumruk
- Institute of Physics
- Faculty of Mathematics and Physics
- Charles University
- 121 16 Prague 2
- Czech Republic
| |
Collapse
|
11
|
Kessler J, Andrushchenko V, Kapitán J, Bouř P. Insight into vibrational circular dichroism of proteins by density functional modeling. Phys Chem Chem Phys 2018; 20:4926-4935. [PMID: 29384537 DOI: 10.1039/c7cp08016f] [Citation(s) in RCA: 39] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Vibrational circular dichroism (VCD) spectroscopy is an excellent method to determine the secondary structure of proteins in solution. Comparison of experimental spectra with quantum-chemical simulations represents a convenient and objective way to extract information on the structure. This has been difficult for such large molecules where approximate theoretical models have to be used. In the present study we applied the Cartesian-coordinate based tensor transfer (CCT) making it possible to extend the density functional theory (DFT) and model spectral intensities of large globular proteins nearly at quantum-chemical precision. Indeed, comparison with experiment provided a better understanding of the dependence of VCD spectral shapes on the geometry, their sensitivity to fine structural details and interactions with the environment. On a model set of globular proteins the simulated spectra correlated well with experimental data and revealed which structural information can (and cannot) be obtained from this kind of spectroscopy. Although the VCD technique has been regarded as being rather insensitive to side-chain variations, we found that the spectra of human and hen lysozyme differing by a few amino acids only are quite distinct. This has been explained by long-distance coupling of the amide vibrations. Likewise, the modeling reproduced some spectral changes caused by protein deuteration even when the protein structure was conserved.
Collapse
Affiliation(s)
- Jiří Kessler
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo náměstí 2, 16610 Prague, Czech Republic.
| | - Valery Andrushchenko
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo náměstí 2, 16610 Prague, Czech Republic.
| | - Josef Kapitán
- Department of Optics, Palacký University, 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Petr Bouř
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo náměstí 2, 16610 Prague, Czech Republic.
| |
Collapse
|
12
|
Brehm M, Thomas M. Computing Bulk Phase Raman Optical Activity Spectra from ab initio Molecular Dynamics Simulations. J Phys Chem Lett 2017; 8:3409-3414. [PMID: 28685571 DOI: 10.1021/acs.jpclett.7b01616] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We present our novel methodology for computing Raman optical activity (ROA) spectra of liquid systems from ab initio molecular dynamics (AIMD) simulations. The method is built upon the recent developments to obtain magnetic dipole moments from AIMD and to integrate molecular properties by using radical Voronoi tessellation. These techniques are used to calculate optical activity tensors for large and complex periodic bulk phase systems. Only AIMD simulations are required as input, and no time-consuming perturbation theory is involved. The approach relies only on the total electron density in each time step and can readily be combined with a wide range of electronic structure methods. To the best of our knowledge, these are the first computed ROA spectra for a periodic bulk phase system. As an example, the experimental ROA spectrum of liquid (R)-propylene oxide is reproduced very well.
Collapse
Affiliation(s)
- Martin Brehm
- Institut für Chemie - Theoretische Chemie, Martin-Luther-Universität Halle-Wittenberg , Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| | - Martin Thomas
- Institut für Chemie - Theoretische Chemie, Martin-Luther-Universität Halle-Wittenberg , Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| |
Collapse
|
13
|
Srebro-Hooper M, Autschbach J. Calculating Natural Optical Activity of Molecules from First Principles. Annu Rev Phys Chem 2017; 68:399-420. [DOI: 10.1146/annurev-physchem-052516-044827] [Citation(s) in RCA: 94] [Impact Index Per Article: 13.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260
| |
Collapse
|
14
|
Sahu N, Gadre SR. Vibrational infrared and Raman spectra of polypeptides: Fragments-in-fragments within molecular tailoring approach. J Chem Phys 2017; 144:114113. [PMID: 27004868 DOI: 10.1063/1.4943966] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
The present work reports the calculation of vibrational infrared (IR) and Raman spectra of large molecular systems employing molecular tailoring approach (MTA). Further, it extends the grafting procedure for the accurate evaluation of IR and Raman spectra of large molecular systems, employing a new methodology termed as Fragments-in-Fragments (FIF), within MTA. Unlike the previous MTA-based studies, the accurate estimation of the requisite molecular properties is achieved without performing any full calculations (FC). The basic idea of the grafting procedure is implemented by invoking the nearly basis-set-independent nature of the MTA-based error vis-à-vis the respective FCs. FIF has been tested out for the estimation of the above molecular properties for three isomers, viz., β-strand, 310- and α-helix of acetyl(alanine)nNH2 (n = 10, 15) polypeptides, three conformers of doubly protonated gramicidin S decapeptide and trpzip2 protein (PDB id: 1LE1), respectively, employing BP86/TZVP, M06/6-311G**, and M05-2X/6-31G** levels of theory. For most of the cases, a maximum difference of 3 cm(-1) is achieved between the grafted-MTA frequencies and the corresponding FC values. Further, a comparison of the BP86/TZVP level IR and Raman spectra of α-helical (alanine)20 and its N-deuterated derivative shows an excellent agreement with the existing experimental spectra. In view of the requirement of only MTA-based calculations and the ability of FIF to work at any level of theory, the current methodology provides a cost-effective solution for obtaining accurate spectra of large molecular systems.
Collapse
Affiliation(s)
- Nityananda Sahu
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208 016, India
| | - Shridhar R Gadre
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208 016, India
| |
Collapse
|
15
|
Krausbeck F, Autschbach J, Reiher M. Calculated Resonance Vibrational Raman Optical Activity Spectra of Naproxen and Ibuprofen. J Phys Chem A 2016; 120:9740-9748. [DOI: 10.1021/acs.jpca.6b09975] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Florian Krausbeck
- ETH Zürich, Laboratorium für Physikalische
Chemie, Vladimir-Prelog-Weg
2, CH-8093 Zürich, Switzerland
| | - Jochen Autschbach
- Department
of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| | - Markus Reiher
- ETH Zürich, Laboratorium für Physikalische
Chemie, Vladimir-Prelog-Weg
2, CH-8093 Zürich, Switzerland
| |
Collapse
|
16
|
Jose KVJ, Raghavachari K. Molecules-in-molecules fragment-based method for the calculation of chiroptical spectra of large molecules: Vibrational circular dichroism and Raman optical activity spectra of alanine polypeptides. Chirality 2016; 28:755-768. [DOI: 10.1002/chir.22651] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2016] [Revised: 08/29/2016] [Accepted: 08/30/2016] [Indexed: 12/17/2022]
Affiliation(s)
- K. V. Jovan Jose
- Department of Chemistry; Indiana University; Bloomington Indiana USA
| | | |
Collapse
|
17
|
Jose KVJ, Beckett D, Raghavachari K. Vibrational Circular Dichroism Spectra for Large Molecules through Molecules-in-Molecules Fragment-Based Approach. J Chem Theory Comput 2015; 11:4238-47. [DOI: 10.1021/acs.jctc.5b00647] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- K. V. Jovan Jose
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Daniel Beckett
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
18
|
Kessler J, Bouř P. Transfer of Frequency-Dependent Polarizabilities: A Tool To Simulate Absorption and Circular Dichroism Molecular Spectra. J Chem Theory Comput 2015; 11:2210-20. [DOI: 10.1021/acs.jctc.5b00136] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jiří Kessler
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo
náměstí 2, 166
10 Prague, Czech Republic
- Department
of Physical and Macromolecular Chemistry, Faculty of Science, Charles University, Hlavova 8, 128 40 Prague, Czech Republic
| | - Petr Bouř
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo
náměstí 2, 166
10 Prague, Czech Republic
| |
Collapse
|
19
|
Bieler NS, Hünenberger PH. On the ambiguity of conformational states: A B&S-LEUS simulation study of the helical conformations of decaalanine in water. J Chem Phys 2015; 142:165102. [DOI: 10.1063/1.4918548] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
|
20
|
Parchaňský V, Kapitán J, Bouř P. Inspecting chiral molecules by Raman optical activity spectroscopy. RSC Adv 2014. [DOI: 10.1039/c4ra10416a] [Citation(s) in RCA: 60] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
|
21
|
Panek PT, Jacob CR. Efficient calculation of anharmonic vibrational spectra of large molecules with localized modes. Chemphyschem 2014; 15:3365-77. [PMID: 25080397 DOI: 10.1002/cphc.201402251] [Citation(s) in RCA: 57] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2014] [Indexed: 11/11/2022]
Abstract
The analysis and interpretation of the vibrational spectra of complex (bio)molecular systems, such as polypeptides and proteins, requires support from quantum-chemical calculations. Such calculations are currently restricted to the harmonic approximation. Here, we show how one of the main bottlenecks in such calculations, the evaluation of the potential energy surface, can be overcome by using localized modes instead of the commonly employed normal modes. We apply such local vibrational self-consistent field (L-VSCF) and vibrational configuration interaction (L-VCI) calculations to a cyclic water tetramer and a helical hexa-alanine peptide. The results show that the use of localized modes is equivalent to the commonly used normal modes, but offers several advantages. First, a faster convergence with respect to the excitation level is observed in L-VCI calculations. Second, the localized modes provide a reduced representation of the couplings between modes that show a regular coupling pattern. This can be used to disregard a significant number of small two-mode potentials a priori. Several such reduced coupling approximations are explored, and we show that the number of single-point calculations required to evaluate the potential energy surface can be significantly reduced without introducing noticeable errors in the resulting vibrational spectra.
Collapse
Affiliation(s)
- Paweł T Panek
- Karlsruhe Institute of Technology (KIT), Center for Functional Nanostructures and Institute of Physical Chemistry, Wolfgang-Gaede-Str. 1a, 76131 Karlsruhe (Germany) http://www.christophjacob.eu
| | | |
Collapse
|
22
|
Kessler J, Keiderling TA, Bouř P. Arrangement of Fibril Side Chains Studied by Molecular Dynamics and Simulated Infrared and Vibrational Circular Dichroism Spectra. J Phys Chem B 2014; 118:6937-45. [DOI: 10.1021/jp502178d] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jiří Kessler
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo
náměstí 2, 166
10 Prague, Czech Republic
- Department
of Physical and Macromolecular Chemistry, Faculty of Science, Charles University, Hlavova 8, 128 40 Prague, Czech Republic
| | - Timothy A. Keiderling
- Department
of Chemistry, University of Illinois at Chicago, 845 West Taylor
Street, Chicago, Illinois 60607-7061, United States
| | - Petr Bouř
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo
náměstí 2, 166
10 Prague, Czech Republic
| |
Collapse
|
23
|
Yamamoto S, Furukawa T, Bouř P, Ozaki Y. Solvated States of Poly-l-alanine α-Helix Explored by Raman Optical Activity. J Phys Chem A 2014; 118:3655-62. [DOI: 10.1021/jp500794s] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Shigeki Yamamoto
- Department
of Chemistry, Graduate School of Science, Osaka University, Toyonaka, Osaka 560-0043, Japan
| | - Tatsuya Furukawa
- Department
of Chemistry, School of Science and Technology, Kwansei Gakuin University, Gakuen 2-1, Sanda, Hyogo 669-1337, Japan
| | - Petr Bouř
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences, 166 10 Prague, Czech Republic
| | - Yukihiro Ozaki
- Department
of Chemistry, School of Science and Technology, Kwansei Gakuin University, Gakuen 2-1, Sanda, Hyogo 669-1337, Japan
| |
Collapse
|
24
|
Nagy PR, Surján PR, Szabados Á. Vibrational optical activity of chiral carbon nanoclusters treated by a generalized π-electron method. J Chem Phys 2014; 140:044112. [DOI: 10.1063/1.4862682] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
25
|
Thomas M, Brehm M, Fligg R, Vöhringer P, Kirchner B. Computing vibrational spectra from ab initio molecular dynamics. Phys Chem Chem Phys 2013; 15:6608-22. [DOI: 10.1039/c3cp44302g] [Citation(s) in RCA: 319] [Impact Index Per Article: 29.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
|
26
|
Thorvaldsen AJ, Gao B, Ruud K, Fedorovsky M, Zuber G, Hug W. Efficient Calculation of ROA Tensors with Analytical Gradients and Fragmentation. Chirality 2012; 24:1018-30. [DOI: 10.1002/chir.22090] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2012] [Accepted: 03/28/2012] [Indexed: 11/10/2022]
Affiliation(s)
- Andreas J. Thorvaldsen
- Centre for Theoretical and Computational Chemistry, Department of Chemistry; University of Tromsø; Tromsø Norway
| | - Bin Gao
- Centre for Theoretical and Computational Chemistry, Department of Chemistry; University of Tromsø; Tromsø Norway
| | - Kenneth Ruud
- Centre for Theoretical and Computational Chemistry, Department of Chemistry; University of Tromsø; Tromsø Norway
| | - Maxim Fedorovsky
- Department of Chemistry; University of Fribourg; Fribourg Switzerland
| | - Gérard Zuber
- Department of Chemistry; University of Fribourg; Fribourg Switzerland
| | - Werner Hug
- Department of Chemistry; University of Fribourg; Fribourg Switzerland
| |
Collapse
|
27
|
Guo W, Wu A, Zhang IY, Xu X. XO: An extended ONIOM method for accurate and efficient modeling of large systems. J Comput Chem 2012; 33:2142-60. [DOI: 10.1002/jcc.23051] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2012] [Revised: 06/04/2012] [Accepted: 06/05/2012] [Indexed: 11/10/2022]
|
28
|
Borowski P, Pasieczna-Patkowska S, Barczak M, Pilorz K. Theoretical determination of the infrared spectra of amorphous polymers. J Phys Chem A 2012; 116:7424-35. [PMID: 22686640 DOI: 10.1021/jp303209v] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The simple procedure of calculating the infrared spectra of polymers is presented. It is based on selecting the relevant, medium-size representative fragments of a polymer, for which the vibrational frequencies are computed within the harmonic approximation, in conjunction with the multiparameter scaling techniques. Scaling is necessary to predict the reliable fundamentals, which, along with the calculated intensities and properly chosen band widths, reproduce the observed band shapes with high accuracy. Applications to the three polymers: poly(methyl methacrylate), poly(vinyl acetate), and poly(isopropenyl acetate) are presented. The simulated spectra are in good agreement with the experiment. The assignment of bands is reported. The obtained results indicate strong delocalization of the vibrational modes within polymers, which is in accord with the most recent experimental finding [Macromolecules2008, 41, 2494-2501]. Good agreement between the observed and the calculated spectra of deuterated PMMA confirms the correctness of our approach. The preliminary results obtained for the highly irregular macromolecular compound (vinyl-functionalized silica) are also shown.
Collapse
Affiliation(s)
- Piotr Borowski
- Faculty of Chemistry, Maria Curie-Skłodowska University, pl. Marii Curie-Skłodowskiej 3, 20-031 Lublin, Poland.
| | | | | | | |
Collapse
|
29
|
Weymuth T, Haag MP, Kiewisch K, Luber S, Schenk S, Jacob CR, Herrmann C, Neugebauer J, Reiher M. MOVIPAC: Vibrational spectroscopy with a robust meta-program for massively parallel standard and inverse calculations. J Comput Chem 2012; 33:2186-98. [DOI: 10.1002/jcc.23036] [Citation(s) in RCA: 53] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2012] [Revised: 05/16/2012] [Accepted: 05/19/2012] [Indexed: 11/08/2022]
|
30
|
Simmen B, Weymuth T, Reiher M. How Many Chiral Centers Can Raman Optical Activity Spectroscopy Distinguish in a Molecule? J Phys Chem A 2012; 116:5410-9. [DOI: 10.1021/jp303428f] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Benjamin Simmen
- Laboratorium
für Physikalische Chemie, ETH Zurich, Wolfgang-Pauli-Strasse 10, 8093 Zurich,
Switzerland
| | - Thomas Weymuth
- Laboratorium
für Physikalische Chemie, ETH Zurich, Wolfgang-Pauli-Strasse 10, 8093 Zurich,
Switzerland
| | - Markus Reiher
- Laboratorium
für Physikalische Chemie, ETH Zurich, Wolfgang-Pauli-Strasse 10, 8093 Zurich,
Switzerland
| |
Collapse
|
31
|
Conformational analyses of peptides and proteins by vibrational Raman optical activity. Anal Bioanal Chem 2012; 403:2203-12. [DOI: 10.1007/s00216-012-5891-0] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2011] [Revised: 02/09/2012] [Accepted: 02/20/2012] [Indexed: 10/28/2022]
|
32
|
Dračínský M, Bouř P. Vibrational averaging of the chemical shift in crystalline α-glycine. J Comput Chem 2012; 33:1080-9. [PMID: 22410968 DOI: 10.1002/jcc.22940] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2011] [Revised: 01/04/2012] [Accepted: 01/04/2012] [Indexed: 11/06/2022]
Abstract
Averaging of the chemical shift over the molecular motion improves the simulated data and provides additional information about the temperature dependence and system dynamics. However, crystal modeling is difficult due to the limited precision of the plane-wave density functional theory (DFT) methods and approximate vibrational schemes. On the glycine example, we investigate how the averaging can be achieved within the periodic boundary conditions at the DFT level. The nuclear motion is modeled with the vibrational configuration interaction, with other simplified quantum anharmonic schemes, and the classical Born-Oppenheimer molecular dynamics (BOMD). The results confirm a large vibrational contribution to the isotropic shielding values. Both the first and second derivatives of the shielding were found important for the quantum averaging. The first derivatives influence the shielding mostly due to the anharmonic character of the CH and NH stretching modes, whereas second derivatives produce most vibrational corrections associated with the lower-frequency vibrational modes. Temperature excitations of the lowest-frequency vibrational states and the expansion of the crystal cell both determine the temperature dependence of nuclear magnetic resonance parameters. The vibrational quantum approach as well as classical BOMD schemes provided temperature dependencies of the chemical shifts that are consistent with the previous experimental data.
Collapse
Affiliation(s)
- Martin Dračínský
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo náměstí 2, Prague 166 10, Czech Republic.
| | | |
Collapse
|
33
|
Yamamoto S, Li X, Ruud K, Bouř P. Transferability of Various Molecular Property Tensors in Vibrational Spectroscopy. J Chem Theory Comput 2012; 8:977-85. [DOI: 10.1021/ct200714h] [Citation(s) in RCA: 57] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Shigeki Yamamoto
- Institute of Organic Chemistry
and Biochemistry, Academy of Sciences,
166 10 Prague, Czech Republic
| | - Xiaojun Li
- Centre
for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø, N-9037 Tromsø, Norway
| | - Kenneth Ruud
- Centre
for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø, N-9037 Tromsø, Norway
| | - Petr Bouř
- Institute of Organic Chemistry
and Biochemistry, Academy of Sciences,
166 10 Prague, Czech Republic
| |
Collapse
|
34
|
Yamamoto S, Kaminský J, Bouř P. Structure and Vibrational Motion of Insulin from Raman Optical Activity Spectra. Anal Chem 2012; 84:2440-51. [DOI: 10.1021/ac2032436] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Shigeki Yamamoto
- Institute of Organic Chemistry and
Biochemistry, Academy of Sciences, Flemingovo
nám. 2, 166 10, Prague 6, Czech Republic
| | - Jakub Kaminský
- Institute of Organic Chemistry and
Biochemistry, Academy of Sciences, Flemingovo
nám. 2, 166 10, Prague 6, Czech Republic
| | - Petr Bouř
- Institute of Organic Chemistry and
Biochemistry, Academy of Sciences, Flemingovo
nám. 2, 166 10, Prague 6, Czech Republic
| |
Collapse
|
35
|
Profant V, Baumruk V, Li X, Šafařík M, Bouř P. Tracking of the Polyproline Folding by Density Functional Computations and Raman Optical Activity Spectra. J Phys Chem B 2011; 115:15079-89. [DOI: 10.1021/jp207706p] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Václav Profant
- Faculty of Mathematics and Physics, Charles University, Ke Karlovu 5, 12116, Prague 2, Czech Republic
| | - Vladimír Baumruk
- Faculty of Mathematics and Physics, Charles University, Ke Karlovu 5, 12116, Prague 2, Czech Republic
| | - Xiaojun Li
- Department of Chemistry, Centre for Theoretical and Computational Chemistry, University of Tromsø, N-9037 Tromsø, Norway
| | - Martin Šafařík
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo nám. 2, 166 10 Prague, Czech Republic
| | - Petr Bouř
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo nám. 2, 166 10 Prague, Czech Republic
| |
Collapse
|
36
|
Jacob CR. Theoretical Study of the Raman Optical Activity Spectra of 310-Helical Polypeptides. Chemphyschem 2011; 12:3291-306. [DOI: 10.1002/cphc.201100593] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2011] [Indexed: 11/06/2022]
|