1
|
Provorse Long MR, Isborn CM. Combining Explicit Quantum Solvent with a Polarizable Continuum Model. J Phys Chem B 2017; 121:10105-10117. [DOI: 10.1021/acs.jpcb.7b06693] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
| | - Christine M. Isborn
- Chemistry
and Chemical Biology, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
2
|
Husseini FS, Robinson D, Hunt NT, Parker AW, Hirst JD. Computing infrared spectra of proteins using the exciton model. J Comput Chem 2017; 38:1362-1375. [PMID: 27868210 PMCID: PMC5434914 DOI: 10.1002/jcc.24674] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2016] [Revised: 10/22/2016] [Accepted: 10/29/2016] [Indexed: 02/02/2023]
Abstract
The ability to compute from first principles the infrared spectrum of a protein in solution phase representing a biological system would provide a useful connection to atomistic models of protein structure and dynamics. Indeed, such calculations are a vital complement to 2DIR experimental measurements, allowing the observed signals to be interpreted in terms of detailed structural and dynamical information. In this article, we have studied nine structurally and spectroscopically well-characterized proteins, representing a range of structural types. We have simulated the equilibrium conformational dynamics in an explicit point charge water model. Using the resulting trajectories based on MD simulations, we have computed the one and two dimensional infrared spectra in the Amide I region, using an exciton approach, in which a local mode basis of carbonyl stretches is considered. The role of solvent in shifting the Amide I band (by 30 to 50 cm-1 ) is clearly evident. Similarly, the conformational dynamics contribute to the broadening of peaks in the spectrum. The inhomogeneous broadening in both the 1D and 2D spectra reflects the significant conformational diversity observed in the simulations. Through the computed 2D cross-peak spectra, we show how different pulse schemes can provide additional information on the coupled vibrations. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Fouad S Husseini
- School of Chemistry, University of Nottingham, Nottingham, NG7 2RD, United Kingdom
| | - David Robinson
- School of Chemistry, University of Nottingham, Nottingham, NG7 2RD, United Kingdom
| | - Neil T Hunt
- Department of Physics, University of Strathclyde, SUPA, 107 Rottenrow East, Glasgow, G4 0NG, Scotland, United Kingdom
| | - Anthony W Parker
- STFC Rutherford Appleton Laboratory, Central Laser Facility, Harwell Campus, Didcot, OX11 0QX, United Kingdom
| | - Jonathan D Hirst
- School of Chemistry, University of Nottingham, Nottingham, NG7 2RD, United Kingdom
| |
Collapse
|
3
|
Hanson-Heine MWD, Husseini FS, Hirst JD, Besley NA. Simulation of Two-Dimensional Infrared Spectroscopy of Peptides Using Localized Normal Modes. J Chem Theory Comput 2016; 12:1905-18. [DOI: 10.1021/acs.jctc.5b01198] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Fouad S. Husseini
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| | - Jonathan D. Hirst
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| | - Nicholas A. Besley
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| |
Collapse
|
4
|
Anderson BA, Literati A, Ball B, Kubelka J. Temperature dependence of C-terminal carboxylic group IR absorptions in the amide I' region. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2015; 134:473-483. [PMID: 25036456 DOI: 10.1016/j.saa.2014.06.102] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/12/2014] [Revised: 05/13/2014] [Accepted: 06/16/2014] [Indexed: 06/03/2023]
Abstract
Studies of structural changes in peptides and proteins using IR spectroscopy often rely on subtle changes in the amide I' band as a function of temperature. However, these changes can be obscured by the overlap with other absorptions, namely the side-chain and terminal carboxylic groups. The former were the subject of our previous report (Anderson et al., 2014). In this paper we investigate the IR spectra of the asymmetric stretch of α-carboxylic groups for amino acids representing all major types (Gly, Ala, Val, Leu, Ser, Thr, Asp, Glu, Lys, Asn, His, Trp, Pro) as well as the C-terminal groups of three dipeptides (Gly-Gly, Gly-Ala, Ala-Gly) in D₂O at neutral pH. Experimental temperature dependent IR spectra were analyzed by fitting of both symmetric and asymmetric pseudo-Voigt functions. Qualitatively the spectra exhibit shifts to higher frequency, loss in intensity and narrowing with increased temperature, similar to that observed previously for the side-chain carboxylic groups of Asp. The observed dependence of the band parameters (frequency, intensity, width and shape) on temperature is in all cases linear: simple linear regression is therefore used to describe the spectral changes. The spectral parameters vary between individual amino acids and show systematic differences between the free amino acids and dipeptides, particularly in the absolute peak frequencies, but the temperature variations are comparable. The relative variations between the dipeptide spectral parameters are most sensitive to the C-terminal amino acid, and follow the trends observed in the free amino acid spectra. General rules for modeling the α-carboxylic IR absorption bands in peptides and proteins as the function of temperature are proposed.
Collapse
Affiliation(s)
- Benjamin A Anderson
- University of Wyoming, 1000 E. University Ave, Laramie, WY 82071, United States
| | - Alex Literati
- University of Wyoming, 1000 E. University Ave, Laramie, WY 82071, United States
| | - Borden Ball
- University of Wyoming, 1000 E. University Ave, Laramie, WY 82071, United States
| | - Jan Kubelka
- University of Wyoming, 1000 E. University Ave, Laramie, WY 82071, United States.
| |
Collapse
|
5
|
Chi H, Welch WRW, Kubelka J, Keiderling TA. Insight into the Packing Pattern of β2 Fibrils: A Model Study of Glutamic Acid Rich Oligomers with 13C Isotopic Edited Vibrational Spectroscopy. Biomacromolecules 2013; 14:3880-91. [DOI: 10.1021/bm401015f] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Heng Chi
- Department
of Chemistry, University of Illinois at Chicago, 845 West Taylor
Street (m/c111), Chicago, Illinois 60607-7061, United States
| | - William R. W. Welch
- Department
of Chemistry, University of Wyoming, Laramie, Wyoming 82071, United States
| | - Jan Kubelka
- Department
of Chemistry, University of Wyoming, Laramie, Wyoming 82071, United States
| | - Timothy A. Keiderling
- Department
of Chemistry, University of Illinois at Chicago, 845 West Taylor
Street (m/c111), Chicago, Illinois 60607-7061, United States
| |
Collapse
|
6
|
Welch WRW, Kubelka J, Keiderling TA. Infrared, vibrational circular dichroism, and Raman spectral simulations for β-sheet structures with various isotopic labels, interstrand, and stacking arrangements using density functional theory. J Phys Chem B 2013; 117:10343-58. [PMID: 23924300 DOI: 10.1021/jp4056126] [Citation(s) in RCA: 56] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Infrared (IR), Raman, and vibrational circular dichroism (VCD) spectral variations for different β-sheet structures were studied using simulations based on density functional theory (DFT) force field and intensity computations. The DFT vibrational parameters were obtained for β-sheet fragments containing nine-amides and constrained to a variety of conformations and strand arrangements. These were subsequently transferred onto corresponding larger β-sheet models, normally consisting of five strands with ten amides each, for spectral simulations. Further extension to fibril models composed of multiple stacked β-sheets was achieved by combining the transfer of DFT parameters for each sheet with dipole coupling methods for interactions between sheets. IR spectra of the amide I show different splitting patterns for parallel and antiparallel β-sheets, and their VCD, in the absence of intersheet stacking, have distinct sign variations. Isotopic labeling by (13)C of selected residues yields spectral shifts and intensity changes uniquely sensitive to relative alignment of strands (registry) for antiparallel sheets. Stacking of multiple planar sheets maintains the qualitative spectral character of the single sheet but evidences some reduction in the exciton splitting of the amide I mode. Rotating sheets with respect to each other leads to a significant VCD enhancement, whose sign pattern and intensity is dependent on the handedness and degree of rotation. For twisted β-sheets, a significant VCD enhancement is computed even for sheets stacked with either the same or opposite alignments and the inter-sheet rotation, depending on the sense, can either further increase or weaken the enhanced VCD intensity. In twisted, stacked structures (without rotation), similar VCD amide I patterns (positive couplets) are predicted for both parallel and antiparallel sheets, but different IR intensity distributions still enable their differentiation. Our simulation results prove useful for interpreting experimental vibrational spectra in terms of β-sheet and fibril structure, as illustrated in the accompanying paper.
Collapse
Affiliation(s)
- William R W Welch
- Department of Chemistry, University of Wyoming , Laramie, Wyoming 82071, United States
| | | | | |
Collapse
|
7
|
Welch WRW, Keiderling TA, Kubelka J. Structural analyses of experimental 13C edited amide I' IR and VCD for peptide β-sheet aggregates and fibrils using DFT-based spectral simulations. J Phys Chem B 2013; 117:10359-69. [PMID: 23924239 DOI: 10.1021/jp405613r] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In the preceding paper, computational models based on density functional theory (DFT) were presented to characterize the sensitivity of vibrational spectroscopic methods (IR, VCD, and Raman) to structural features of β-sheets. Isotopically edited amide I' IR for peptides labeled with (13)C in multiple different sites provides the most structurally distinct signatures of strand alignment, while VCD is sensitive to the sheet twist and intersheet stacking. In this report, we simulate the IR and VCD spectra for models approximating structures of four β-sheet forming peptides previously experimentally studied using these methods with (13)C isotopic editing. Various register alignments are tested. Agreement with experiment is evaluated based on frequency shifts of both the (12)C and (13)C IR amide I' signals, relative intensity patterns, and VCD spectra where available. While for the simulation of IR spectra canonical planar sheets provide a sufficient model system, for VCD simulation twisted, stacked sheets are required in order to reproduce strong couplet-like amide I' VCD. Effects of the solvent (water) and amino acid side chains are also tested by using a simplified, electrostatic solvent model and atomic partial charges for the side chains. Very good agreement with experimental spectra is obtained, particularly for the relative (12)C and (13)C band frequencies. All four peptide models are shown to be antiparallel as had previously been assumed. However, in some cases our simulations are consistent with different register alignment of strands than originally proposed.
Collapse
Affiliation(s)
- William R W Welch
- Department of Chemistry, University of Wyoming , Laramie, Wyoming 82071, United States
| | | | | |
Collapse
|