Dissecting amide-I vibration in β-peptide helices

Vibrational Energy Transfer in Energetic Ionic Liquid 4-Amino-1H-1,2,4-triazolium Nitrate: Ab Initio Molecular Dynamics Simulations

Energetic ionic liquids (EILs) represent a distinctive class of energetic materials with substantial research significance and promising energetic applications. In this work, we delved into the vibrational energy transfer mechanism within the EILs, specifically focusing on 4-amino-1H-1,2,4-triazolium nitrate (ATN), utilizing ab initio molecular dynamics simulations. Our work illustrates distinct energy transfer patterns for different vibrational modes. Upon exciting the stretching vibration of the NH group in the cationic group, vibrational energy preferentially migrates to the neighboring CH bond within the aromatic ring on the femtosecond to picosecond time scales and notably in an in-phase coherent energy transfer fashion. In contrast, exciting the stretching vibration of the N9H11 bond triggers the transfer of vibrational energy to its neighboring N9H10 bond in an out-of-phase coherent fashion. Conversely, exciting the stretching vibration of the N9H10 bond leads to energy transfer predominantly through intermolecular pathways due to the hydrogen-bonding interaction between this bond and the anion. The vibrational energy of the excited N9H10 stretch is shown to dissipate very rapidly, displaying a fast component (with a time constant as short as ca. 7 fs) and a slow component (ca. 230 fs).

Modeling amide-I vibrations of alanine dipeptide in solution by using neural network protocol

Infrared spectroscopy is a powerful tool for the understanding of molecular structure and function of polypeptides. Theoretical interpretation of IR spectra relies on ab initio calculations may be very costly in computational resources. Herein, we developed a neural network (NN) modeling protocol to evaluate a model dipeptide's backbone amide-I spectra. DFT calculations were performed for the amide-I vibrational motions and structural parameters of alanine dipeptide (ALAD) conformers in different micro-environments ranging from polar to non-polar ones. The obtained backbone dihedrals, C = O bond lengths and amide-I frequencies of ALAD were gather together for NN architecture. The applications of built NN protocols for the prediction of amide-I frequencies of ALAD in other solvation conditions are quite satisfactory with much less computational cost comparing with electronic structure calculations. The results show that this cost-effective way enables us to decipher the polypeptide's dynamic secondary structures and biological functions with their backbone vibrational probes.

Structural dynamics and vibrational feature of N-Acetyl-d-glucosamine in aqueous solution

Molecular dynamics simulations and DFT calculations were performed for the demonstration of the structural dynamics and vibrational feature of N-Acetyl-d-glucosamine (NAG) in solution phase. The interactions between NAG and solvent molecules were evaluated through spatial distribution function and radial distribution function, and the preferred conformations of NAG in aqueous solution were revealed by cluster analysis. Results from normal mode analysis show that the solvent induced structural fluctuation of NAG could be reflected in the vibrational feature of specific chromophores, thus we can evaluate the molecular structure with the help of its vibrational signature based on the built correlation between molecular structure and vibrational frequencies of specific groups.

Deciphering the structural preference encoded in amide-I vibrations of lysine dipeptide in gas phase and in aqueous solution

Protein's biological function is critically associated with its structural feature, which is encoded in its amino acid sequence. For evaluation of conformational fluctuation and folding mechanism, DFT calculations were performed on the model compound - lysine dipeptide (LYSD) in gas phase to demonstrate the correlation between amide-I vibrations and secondary structure. Molecular dynamics simulations were carried out for the structural dynamics of LYSD in aqueous solution. The results show that LYSD tends form C_7eq, C₅, β, PP_II and α conformations in the gas phase and primarily presented PP_II and α conformations in aqueous solution. The obtained amide-I vibrational frequencies of LYSD conformers were assigned, thus build the correlations between amide-I probes and secondary structure of LYSD. These results provide theoretical insights into the structural feature of LYSD through amide-I vibrations, and would shed light on site specific structural prediction of polypeptides.

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.

Application of a transparent window vibrational probe (azido probe) to the structural dynamics of model dipeptides and amyloid β-peptide

The azido asymmetric stretching motion is widely used for the elucidation of the intrinsic conformational preference and folding mechanism of protein since it has strong vibrational absorbance in the spectral transparent windows. However, the possible secondary structural disturbance induced by the insertion of azido group in the side chain of polypeptides should be carefully evaluated. Here, DFT calculation and enhanced sampling method were employed for model dipeptides with or without azido substitution, and the outcome results show that the lower potential energy basins of isolated model dipeptides are consistent with the preferred structural distributions of model dipeptides in aqueous solution. The azido asymmetric stretching frequency shows its sensitivity to the backbone configurations just like amide-I vibration does, and the azido vibration exhibits great potential as a structural reporter in the transparent window. For the evaluation of the application of azido group in biologically related system, the structural dynamics of Aβ_37-42 and N₃-Aβ_37-42 fragments and the self-assemble process of their protofiliments in aqueous solution were demonstrated. The outcome results show that the structural fluctuations of Aβ_37-42 and its protofilament in aqueous solution are quite similar with or without azido substitution, and the dewetting transitions of Aβ_37-42 and N₃-Aβ_37-42 β-sheet layers are both complete within 30 ns and assemble into stable protofilaments. Therefore, the azido asymmetric vibrational motion is a minimally invasive structural probe and would not introduce much disturbance to the structural dynamics of polypeptides.

Mapping the amide-I vibrations of model dipeptides with secondary structure sensitivity and amino acid residue specificity, and its application to amyloid β-peptide in aqueous solution

Vibrational spectroscopy has been known as particularly well-suited for deciphering the polypeptide's structure. To decode structural information encoded in IR spectra, we developed amide-I frequency maps on the basis of model dipeptides to correlate the amide-I frequency of interest to the combination of the calculated secondary structure dependent amide-I frequency by using DFT method and the electrostatic potentials that projected onto the amide unit from the micro-environment within molecular mechanics force field. The constructed maps were applied to model dipeptides and amyloid β-peptide fragment (Aβ_25-35). The dipeptide specified map (DS map) and the hybrid map (HYB map) predicted amide-I bands of Aβ_25-35 in solution satisfactorily reproduce experimental observation, and indicate the preference of forming β-sheet and random coil structure for Aβ_25-35 in D₂O just as the results of cluster analysis suggested. These maps with secondary structural sensitivity and amino acid residue specificity open up a way for the interpretation of amide-I vibrations and show their potentials in the understanding of molecular structure of polypeptides in solution.

Central-metal effect on intramolecular vibrational energy transfer of M(CO)5Br (M = Mn, Re) probed by two-dimensional infrared spectroscopy

Vibrational energy transfer in transition metal complexes with flexible structures in condensed phases is of central importance to catalytical chemistry processes. In this work, two molecules with different metal atoms, M(CO)₅Br (where M = Mn, Re), were used as model systems, and their axial and radial carbonyl stretching modes as infrared probes. The central-metal effect on intramolecular vibrational energy redistribution (IVR) in M(CO)₅Br was investigated in polar and nonpolar solvents. The linear infrared (IR) peak splitting between carbonyl vibrations increases as the metal atom changes from Mn to Re. The waiting-time dependent two-dimensional infrared diagonal- and off-diagonal peak amplitudes reveal a faster IVR process in Re(CO)₅Br than in Mn(CO)₅Br. With the aid of density functional theory (DFT) calculations, the central-metal effect on IVR time linearly correlates with the vibrational coupling strength between the two involved modes. In addition, the polar solvent is found to accelerate the IVR process by affecting the anharmonic vibrational potentials of a solute vibration mode.

Electrostatic frequency maps for amide-I mode of β-peptide: Comparison of molecular mechanics force field and DFT calculations

The spectroscopy of amide-I vibrations has been widely utilized for the understanding of dynamical structure of polypeptides. For the modeling of amide-I spectra, two frequency maps were built for β-peptide analogue (N-ethylpropionamide, NEPA) in a number of solvents within different schemes (molecular mechanics force field based, GM map; DFT calculation based, GD map), respectively. The electrostatic potentials on the amide unit that originated from solvents and peptide backbone were correlated to the amide-I frequency shift from gas phase to solution phase during map parameterization. GM map is easier to construct with negligible computational cost since the frequency calculations for the samples are purely based on force field, while GD map utilizes sophisticated DFT calculations on the representative solute-solvent clusters and brings insight into the electronic structures of solvated NEPA and its chemical environments. The results show that the maps' predicted amide-I frequencies present solvation environmental sensitivities and exhibit their specific characters with respect to the map protocols, and the obtained vibrational parameters are in satisfactory agreement with experimental amide-I spectra of NEPA in solution phase. Although different theoretical schemes based maps have their advantages and disadvantages, the present maps show their potentials in interpreting the amide-I spectra for β-peptides, respectively.

General Applicable Frequency Map for the Amide-I Mode in β-Peptides

In this work, a general applicable amide-I vibrational frequency map (GA map) for β-peptides in a number of common solvents was constructed, based on a peptide derivative, N-ethylpropionamide (NEPA). The map utilizes force fields at the ab initio computational level to accurately describe molecular structure and solute-solvent interactions, and also force fields at the molecular mechanics level to take into account long-range solute-solvent interactions. The results indicate that the GA map works reasonably for mapping the vibrational frequencies of the amide-I local-modes for β-peptides, holding promises for understanding the complicated infrared spectra of the amide-I mode in β-polypeptides.