Calculations of intermode coupling constants and simulations of amide I, II, and III vibrational spectra of dipeptides

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.

Effect of radiotherapy on the chemical composition of root dentin

BACKGROUND

The radiotherapy can directly affect the bond strength of the adhesive materials, interfering in the prognosis of restorative treatments, which may be caused by chemical changes in dentin structure.

METHODS

Twenty inferior homologues premolars were distributed in 2 groups (in vitro study) (n = 10): nonirradiated and irradiated. The specimens were submitted to the analysis of phosphate (ν₁ PO₄ ^3- ;ν₂ PO₄ ^3- ;ν₄ PO₄ ^3- ), carbonate (ν₃ CO₃ ^2- ), amide I, CH₂ , amide III, and amide I/III ratio by confocal Raman spectroscopy. Data were submitted to statistical analysis (T test, P < .05).

RESULTS

In intracanal dentin, the irradiated group had lower ν₄ PO₄ ^3- values (1.23 ± 0.06) compared to nonirradiated group (1.40 ± 0.18) (P < .05), with no difference for ν₁ PO₄ ^3- and ν₂ PO₄ ³ peaks (P > .05). The irradiated (1.56 ± 0.06) had lower carbonate, amide III (1.05 ± 0.19), and amide I/III ratio values (0.19 ± 0.06) compared to nonirradiated group (1.42 ± 0.10, 1.28 ± 0.24, and 0.31 ± 0.10, respectively) (P < .05). For medium dentin irradiated group (1.30 ± 0.12) had lower phosphate values compared to nonirradiated group (1.48 ± 0.22) (P < .05). In cementum, there was no statistical difference between the groups.

CONCLUSION

The radiotherapy was able to cause changes in ν₄ PO₄ ^3- , carbonate, and amide III peaks of root dentin.

In vivo study of dermal collagen of striae distensae by confocal Raman spectroscopy

This research work mainly deals with studying qualitatively the changes in the dermal collagen of two forms of striae distensae (SD) namely striae rubrae (SR) and striae albae (SA) when compared to normal skin (NS) using confocal Raman spectroscopy. The methodology includes an in vivo human skin study for the comparison of confocal Raman spectra of dermis region of SR, SA, and NS by supervised multivariate analysis using partial least squares discriminant analysis (PLS-DA) to determine qualitatively the changes in dermal collagen. These groups are further analyzed for the extent of hydration of dermal collagen by studying the changes in the water content bound to it. PLS-DA score plot showed good separation of the confocal Raman spectra of dermis region into SR, SA, and NS data groups. Further analysis using loading plot and S-plot indicated the participation of various components of dermal collagen in the separation of these groups. Bound water content analysis showed that the extent of hydration of collagen is more in SD when compared to NS. Based on the results obtained, this study confirms the active involvement of dermal collagen in the formation of SD. It also emphasizes the need to study quantitatively the role of these various biochemical changes in the dermal collagen responsible for the variance between SR, SA, and NS.

Quantifying Carboxylic Acid Concentration in Model Polyamide Desalination Membranes via Fourier Transform Infrared Spectroscopy

Carboxylic acid groups impart hydrophilicity and ionizable moieties to polyamide membranes for desalination, hence influencing water and ion transport through the material. Model polyamide films were synthesized via molecular layer-by-layer deposition on planar substrates to study the formation process of these materials and overcome the chemical and topological inhomogeneity inherent to conventional interfacially polymerized polyamide membranes. The carboxylic acid content in these model films was characterized using Fourier transform infrared (FTIR) spectroscopy by quantifying the C=O band at 1718 cm-1. The concentration of carboxylic acid groups decreased as the thickness of the membrane increased, suggestive of an increase in crosslink density as the polyamide network develops. For the thinnest molecular layer-by-layer (mLbL) samples, the carboxylic acid concentration for films on gold was 0.35 mmol g-1, whereas analogous films on silicon had an acid content of 0.56 mmol g-1, indicating a clear influence of the substrate on the initial network formation. As the thickness of the membrane increased, the influence of the substrate and initial layer growth became less significant as the carboxylic acid concentration on both substrates reached a value of 0.12 mmol g-1. We demonstrate that FTIR spectroscopy is a practical and accessible way to quantify the carboxylic acid content in these types of extremely thin polyamide membranes to help quantify network formation in these materials.

Effective representation of amide III, II, I, and A modes on local vibrational modes: Analysis of ab initio quantum calculation results

The Hamiltonian matrix for the first excited vibrational states of a protein can be effectively represented by local vibrational modes constituting amide III, II, I, and A modes to simulate various vibrational spectra. Methods for obtaining the Hamiltonian matrix from ab initio quantum calculation results are discussed, where the methods consist of three steps: selection of local vibrational mode coordinates, calculation of a reduced Hessian matrix, and extraction of the Hamiltonian matrix from the Hessian matrix. We introduce several methods for each step. The methods were assessed based on the density functional theory calculation results of 24 oligopeptides with four different peptide lengths and six different secondary structures. The completeness of a Hamiltonian matrix represented in the reduced local mode space is improved by adopting a specific atom group for each amide mode and reducing the effect of ignored local modes. The calculation results are also compared to previous models using C=O stretching vibration and transition dipole couplings. We found that local electric transition dipole moments of the amide modes are mainly bound on the local peptide planes. Their direction and magnitude are well conserved except amide A modes, which show large variation. Contrary to amide I modes, the vibrational coupling constants of amide III, II, and A modes obtained by analysis of a dipeptide are not transferable to oligopeptides with the same secondary conformation because coupling constants are affected by the surrounding atomic environment.

Assessment of amide I spectroscopic maps for a gas-phase peptide using IR-UV double-resonance spectroscopy and density functional theory calculations

The spectroscopy of amide I vibrations has become a powerful tool for exploring protein structure and dynamics. To help with spectral interpretation, it is often useful to perform molecular dynamics (MD) simulations. To connect spectroscopic experiments to simulations in an efficient manner, several researchers have proposed "maps," which relate observables in classical MD simulations to quantum spectroscopic variables. It can be difficult to discern whether errors in the theoretical results (compared to experiment) arise from inaccuracies in the MD trajectories or in the maps themselves. In this work, we evaluate spectroscopic maps independently from MD simulations by comparing experimental and theoretical spectra for a single conformation of the α-helical model peptide Ac-Phe-(Ala)5-Lys-H(+) in the gas phase. Conformation-specific experimental spectra are obtained for the unlabeled peptide and for several singly and doubly (13)C-labeled variants using infrared-ultraviolet double-resonance spectroscopy, and these spectra are found to be well-modeled by density functional theory (DFT) calculations at the B3LYP/6-31G** level. We then compare DFT results for the deuterated and (13)C(18)O-labeled peptide with those from spectroscopic maps developed and used previously by the Skinner group. We find that the maps are typically accurate to within a few cm(-1) for both frequencies and couplings, having larger errors only for the frequencies of terminal amides.

The intrinsic conformational features of amino acids from a protein coil library and their applications in force field development

The local conformational (φ, ψ, χ) preferences of amino acid residues remain an active research area, which are important for the development of protein force fields. In this perspective article, we first summarize spectroscopic studies of alanine-based short peptides in aqueous solution. While most studies indicate a preference for the P(II) conformation in the unfolded state over α and β conformations, significant variations are also observed. A statistical analysis from various coil libraries of high-resolution protein structures is then summarized, which gives a more coherent view of the local conformational features. The φ, ψ, χ distributions of the 20 amino acids have been obtained from a protein coil library, considering both backbone and side-chain conformational preferences. The intrinsic side-chain χ(1) rotamer preference and χ(1)-dependent Ramachandran plot can be generally understood by combining the interaction of the side-chain Cγ/Oγ atom with two neighboring backbone peptide groups. Current all-atom force fields such as AMBER ff99sb-ILDN, ff03 and OPLS-AA/L do not reproduce these distributions well. A method has been developed by combining the φ, ψ plot of alanine with the influence of side-chain χ(1) rotamers to derive the local conformational features of various amino acids. It has been further applied to improve the OPLS-AA force field. The modified force field (OPLS-AA/C) reproduces experimental (3)J coupling constants for various short peptides quite well. It also better reproduces the temperature-dependence of the helix-coil transition for alanine-based peptides. The new force field can fold a series of peptides and proteins with various secondary structures to their experimental structures. MD simulations of several globular proteins using the improved force field give significantly less deviation (RMSD) to experimental structures. The results indicate that the local conformational features from coil libraries are valuable for the development of balanced protein force fields.

Parallel β-sheet vibrational couplings revealed by 2D IR spectroscopy of an isotopically labeled macrocycle: quantitative benchmark for the interpretation of amyloid and protein infrared spectra

Infrared spectroscopy is playing an important role in the elucidation of amyloid fiber formation, but the coupling models that link spectra to structure are not well tested for parallel β-sheets. Using a synthetic macrocycle that enforces a two stranded parallel β-sheet conformation, we measured the lifetimes and frequency for six combinations of doubly (13)C═(18)O labeled amide I modes using 2D IR spectroscopy. The average vibrational lifetime of the isotope labeled residues was 550 fs. The frequencies of the labels ranged from 1585 to 1595 cm(-1), with the largest frequency shift occurring for in-register amino acids. The 2D IR spectra of the coupled isotope labels were calculated from molecular dynamics simulations of a series of macrocycle structures generated from replica exchange dynamics to fully sample the conformational distribution. The models used to simulate the spectra include through-space coupling, through-bond coupling, and local frequency shifts caused by environment electrostatics and hydrogen bonding. The calculated spectra predict the line widths and frequencies nearly quantitatively. Historically, the characteristic features of β-sheet infrared spectra have been attributed to through-space couplings such as transition dipole coupling. We find that frequency shifts of the local carbonyl groups due to nearest neighbor couplings and environmental factors are more important, while the through-space couplings dictate the spectral intensities. As a result, the characteristic absorption spectra empirically used for decades to assign parallel β-sheet secondary structure arises because of a redistribution of oscillator strength, but the through-space couplings do not themselves dramatically alter the frequency distribution of eigenstates much more than already exists in random coil structures. Moreover, solvent exposed residues have amide I bands with >20 cm(-1) line width. Narrower line widths indicate that the amide I backbone is solvent protected inside the macrocycle. This work provides calculated and experimentally verified couplings for parallel β-sheets that can be used in structure-based models to simulate and interpret the infrared spectra of β-sheet containing proteins and protein assemblies, such as amyloid fibers.