2D IR Spectra Research Articles

Signatures of coherent vibrational dynamics in ethylene carbonate.

Despite having practical applications in battery technology and serving as a model system for Fermi resonance coupling, ethylene carbonate (EC) receives little direct attention as a vibrational probe in nonlinear vibrational spectroscopy experiments. EC contains a Fermi resonance that is well-characterized in the linear spectrum, and the environmental sensitivity of its Fermi resonance peaks could make it a good molecular probe for two-dimensional infrared spectroscopy (2DIR) experiments. As a model system, we investigate the linear and 2DIR vibrational spectrum of the carbonyl stretching region of ethylene carbonate in tetrahydrofuran. The 2DIR spectrum reveals peak dynamics that evolve coherently. We characterize these dynamics in the context of Redfield theory and find evidence that EC dynamics proceed through coherent pathways, including singular coherence transfer pathways that have not been widely observed in other studies. We find that coherent contributions play a significant role in the observed dynamics of cross-peaks in the 2DIR spectrum, which must be accounted for to extract accurate measurements of early waiting time dynamics.

Co-deposition of tannic acid and bottle-brush polymer to construct stable antifouling surface via hydrogen-bonding

The convenient preparation and excellent stability are two specific demands for antifouling surfaces. We designed a hydrogen-bonded network composed of tannic acid (TA) and bottle-brush polymer (BBP) to form an antifouling coating on diverse substrates. The BBP contained polyurethane backbone and polyethylene glycol (PEG) side chains which provided abundant hydrogen bond acceptors. Molecular dynamics (MD) simulation, 2D-IR and 1H NMR spectra confirmed that TA could interact with both the backbone and side chains of BBP. Owing to the strong affinity of TA to substrates and the hierarchical structure of bottle-brush polymer, TA/BBP coating possessed more intricate non-covalent network than TA/linear polymer coating and was more stable at physiologic pH. The excellent antifouling properties against serum albumin, red blood cells, fibroblasts and bacteria were demonstrated using different techniques. This work provided a strategy to prepare a stable antifouling coating relying solely on the hydrogen-bonded network of TA and BBP. Additionally, TA/BBP coating could be completely detached by strongly basic solution due to the dissociation of TA, which provided a possibility of cleaning and recycling substrates if necessary.

Slowdown of solvent structural dynamics in aqueous DMF solutions

This study presents a comprehensive investigation into the molecular dynamics of solvation environments through an integrated approach combining Fourier-transform infrared (FTIR) spectroscopy, molecular dynamics (MD) simulations, and two-dimensional infrared (2D IR) spectroscopy. We explore the solvation of an ionic solute (ammonium thiocyanate) in various solvent systems, including N,N-dimethylformamide (DMF), water, and a 0.5 mole fraction of DMF in water, aiming to unravel the intricate interplay between solute-solvent interactions and solvent dynamics across diverse solvation environments. By integrating FTIR spectral analysis with radial distribution functions and coordination numbers obtained from MD simulations, we decipher the solvent composition around the solute molecule. Analysis of 2D IR spectra and hydrogen bond, as well as dipolar autocorrelation function from MD simulations, further elucidates the nuances of solute-solvent interactions, highlighting the impact of solvent dynamics on solvation structures. Our results reveal a significant slowdown of the solvent structural dynamics in the equimolar binary solvent mixture compared to the neat solvents. This slowdown underscores the complex relationship between solute-solvent interactions and solvent dynamics. The integration of FTIR, MD simulations, and 2D IR spectroscopy provides a unified framework for obtaining a holistic understanding of solvation dynamics, offering valuable insights into the underlying molecular mechanisms governing solute-solvent interactions in complex systems. These results pave the way for future studies to delve deeper into the molecular intricacies of solvation phenomena.

Nanoclustering in non-ideal ethanol/heptane solutions alters solvation dynamics.

Alcohol/alkane solutions widely used in chemical synthesis and as transportation fuels are highly non-ideal due to the nanoscale clustering of the amphiphilic alcohol molecules within the nonpolar alkanes. Besides impacting reactivity, such as combustion, non-ideal solutions are likely to exhibit unusual solvation dynamics on ultrafast time scales arising from the structurally heterogeneous nature of molecular-scale association. Using a convenient transition metal carbonyl vibrational probe [(C5H5)Mn(CO)3, CMT], linear absorption and nonlinear two-dimensional infrared (2D-IR) spectroscopy reveal composition-dependent solvation dynamics as reported by the frequency fluctuation correlation function in a series of ethanol/heptane solutions. Slow spectral diffusion with dilute ethanol indicates preferential solvation of the polar solute by the alcohol with a mechanism largely dominated by solvent exchange. Comparison with an ethanol/acetonitrile solution series yields no substantial preferential solvation or solvent exchange signatures in the linear or 2D-IR spectra. In ethanol/heptane solutions, increasing the ethanol concentration speeds up the solvation dynamics, which is largely consistent with a model that includes solvent exchange and single-solvent spectral diffusion. Detailed analysis of the deviation from the experimental time constants from the model's optimal parameters yields a remarkable resemblance of the concentration-weighted Kirkwood-Buff integrals for ethanol/heptane solutions. This trend indicates that solution non-ideality alters the spectral diffusion dynamics of the probe solute. Given that nanoscale clustering drives the non-ideality, these experiments reveal a dynamical consequence of nanoscale heterogeneity on the ultrafast dynamics of the solution. Refined understanding of the structural and dynamical aspects of mixed solvents will be necessary for predictive solution strategies in chemistry.

Uncovering the binding nature of thiocyanate in contact ion pairs with lithium ions.

Ion pair formation is a fundamental molecular process that occurs in a wide variety of systems, including electrolytes, biological systems, and materials. In solution, the thiocyanate (SCN-) anion interacts with cations to form contact ion pairs (CIPs). Due to its ambidentate nature, thiocyanate can bind through either its sulfur or nitrogen atoms, depending on the solvent. This study focuses on the binding nature of thiocyanate with lithium ions as a function of the solvents using FTIR, 2D infrared spectroscopy (2DIR) spectroscopies, and theoretical calculations. The study reveals that the SCN- binding mode (S or N end) in CIPs can be identified through 2DIR spectroscopy but not by linear IR spectroscopy. Linear IR spectroscopy shows that the CN stretch frequencies are too close to one another to separate N- and S-bound CIPs. Moreover, the IR spectrum shows that the S-C stretch presents different frequencies for the salt in different solvents, but it is related to the anion speciation rather than to its binding mode. A similar trend is observed for the anion bend. 2DIR spectra show different dynamics for N-bound and S-bound thiocyanate. In particular, the frequency-frequency correlation function (FFCF) dynamics extracted from the 2DIR spectra have a single picosecond exponential decay for N-bound thiocyanate and a biexponential decay for S-bound thiocyanate, consistent with the binding mode of the anion. Finally, it is also observed that the binding mode also affects the line shape parameters, probably due to the different molecular mechanisms of the FFCF for N- and S-bound CIPs.

Understanding the Vibrational Structure and Ultrafast Dynamics of the Metal Carbonyl Precatalyst [Mn(ppy)(CO)4

The solution phase structure, vibrational spectroscopy, and ultrafast relaxation dynamics of the precatalyst species [Mn(ppy)(CO)4] (1) in solution have been investigated using ultrafast two-dimensional infrared (2D-IR) spectroscopy. By comparing 2D-IR data with the results of anharmonic density functional theory (DFT) calculations, we establish an excellent agreement between measured and predicted inter-mode couplings of the carbonyl stretching vibrational modes of 1 that relates to the atomic displacements of axial and equatorial ligands in the modes and the nature of the molecular orbitals involved in M-CO bonding. Measurements of IR pump-probe spectra and 2D-IR spectra as a function of waiting time reveal the presence of ultrafast (few ps) intramolecular vibrational energy redistribution between carbonyl stretching modes prior to vibrational relaxation. The vibrational relaxation times of the CO-stretching modes of 1 are found to be relatively solvent-insensitive, suggestive of limited solvent-solute interactions in the ground electronic state. Overall, these data provide a detailed picture of the complex potential energy surface, bonding and vibrational dynamics of 1, establishing a fundamental basis for the next steps in understanding and modulating precatalyst behavior.

Unraveling dynamic protein structures by two-dimensional infrared spectra with a pretrained machine learning model

Dynamic protein structures are crucial for deciphering their diverse biological functions. Two-dimensional infrared (2DIR) spectroscopy stands as an ideal tool for tracing rapid conformational evolutions in proteins. However, linking spectral characteristics to dynamic structures poses a formidable challenge. Here, we present a pretrained machine learning model based on 2DIR spectra analysis. This model has learned signal features from approximately 204,300 spectra to establish a "spectrum-structure" correlation, thereby tracing the dynamic conformations of proteins. It excels in accurately predicting the dynamic content changes of various secondary structures and demonstrates universal transferability on real folding trajectories spanning timescales from microseconds to milliseconds. Beyond exceptional predictive performance, the model offers attention-based spectral explanations of dynamic conformational changes. Our 2DIR-based pretrained model is anticipated to provide unique insights into the dynamic structural information of proteins in their native environments.

Study of the Myosin Relay Helix Peptide by Molecular Dynamics Simulations, Pump-Probe and 2D Infrared Spectroscopy.

The Davydov model was conjectured to describe how an amide I excitation created during ATP hydrolysis in myosin might be significant in providing energy to drive myosin's chemomechanical cycle. The free energy surfaces of the myosin relay helix peptide dissolved in 2,2,2-trifluoroethanol (TFE), determined by metadynamics simulations, demonstrate local minima differing in free energy by only ~2 kT, corresponding to broken and stabilized hydrogen bonds, respectively. Experimental pump-probe and 2D infrared spectroscopy were performed on the peptide dissolved in TFE. The relative heights of two peaks seen in the pump-probe data and the corresponding relative volumes of diagonal peaks seen in the 2D-IR spectra at time delays between 0.5 ps and 1 ps differ noticeably from what is seen at earlier or later time delays or in the linear spectrum, indicating that a vibrational excitation may influence the conformational state of this helix. Thus, it is possible that the presence of an amide I excitation may be a direct factor in the conformational state taken on by the myosin relay helix following ATP hydrolysis in myosin.

Probing Backbone Coupling within Hydrated Proteins with Two-Color 2D Infrared Spectroscopy.

The vibrational coupling between protein backbone modes and the role of water interactions are important topics in biomolecular spectroscopy. Our work reports the first study of the coupling between amide I and amide A modes within peptides and proteins with secondary structure and water contacts. We use two-color two-dimensional infrared (2D IR) spectroscopy and observe cross peaks between amide I and amide A modes. In experiments with peptides with different secondary structures and side chains, we observe that the spectra are sensitive to secondary structure. Water interactions affect the cross peaks, which may be useful as probes for the accessibility of protein sites to hydration water. Moving to two-color 2D IR spectra of proteins, the data demonstrate that the cross peaks integrate the sensitivities of both amide I and amide A spectra and that a two-color detection scheme may be a promising tool for probing secondary structures in proteins.

Experimental two-dimensional infrared spectra of methyl thiocyanate in water and organic solvents.

Thiocyanates, nitriles, and azides represent a versatile set of vibrational probes to measure the structure and dynamics in biological systems. The probes are minimally perturbative, the nitrile stretching mode appears in an otherwise uncongested spectral region, and the spectra report on the local environment around the probe. Nitrile frequencies and lineshapes, however, are difficult to interpret, and theoretical models that connect local environments with vibrational frequencies are often necessary. However, the development of both more accurate and intuitive models remains a challenge for the community. The present work provides an experimentally consistent collection of experimental measurements, including IR absorption and ultrafast two-dimensional infrared (2D IR) spectra, to serve as a benchmark in the development of future models. Specifically, we catalog spectra of the nitrile stretching mode of methyl thiocyanate (MeSCN) in fourteen different solvents, including non-polar, polar, and protic solvents. Absorption spectra indicate that π-interactions may be responsible for the line shape differences observed between aromatic and aliphatic alcohols. We also demonstrate that a recent Kamlet-Taft formulation describes the center frequency MeSCN. Furthermore, we report cryogenic infrared spectra that may lead to insights into the peak asymmetry in aprotic solvents. 2D IR spectra measured in protic solvents serve to connect hydrogen bonding with static inhomogeneity. We expect that these insights, along with the publicly available dataset, will be useful to continue advancing future models capable of quantitatively describing the relation between local environments, line shapes, and dynamics in nitrile probes.

The 2D-IR spectrum of hydrogen-bonded silanol groups in pyrogenic silica.

Pyrogenic silica is a form of amorphous silica with a high surface area and a heterogeneous distribution of silanol hydroxyl terminations and defects. In this work, the interesting and unusual form of the hydroxyl-stretch 2D-IR spectrum of pyrogenic silica is presented and explored in the deuterated (deuteroxyl) form. Transition dipole couplings between hydrogen-bonded and non-hydrogen-bonded silanol groups give a distinct cross-peak in the 2D-IR spectrum, displaying interstate coherence oscillations during the 2D-IR experimental waiting time. The strong asymmetry about the diagonal is proposed to be the result of both the relatively small transition dipole coupling strength and the extreme differences in the width of the hydrogen-bonded and non-hydrogen-bonded silanol bands. The resulting interference of negative and positive cross-peaks has minimal intensity in the below-diagonal ω3 < ω1 region of the spectrum. An additional strong positive cross-peak is observed at a position in the 2D-IR spectrum inconsistent with transition dipole coupling. An assignment as a fifth order effect is proposed.

The secondary structure of diatom silaffin peptide R5 determined by two-dimensional infrared spectroscopy.

Diatoms, unicellular marine organisms, harness short peptide repeats of the protein silaffin to transform silicic acid into biosilica nanoparticles. This process has been a white whale for material scientists due to its potential in biomimetic applications, ranging from medical to microelectronic fields. Replicating diatom biosilicification will depend on a thorough understanding of the silaffin peptide structure during the reaction, yet existing models in the literature offer conflicting views on peptide folding during silicification. In our study, we employed two-dimensional infrared spectroscopy (2DIR) within the amide I region to determine the secondary structure of the silaffin repeat unit 5 (R5), both pre- and post-interaction with silica. The 2DIR experiments are complemented by molecular dynamics (MD) simulations of pure R5 reacting with silicate. Subsequently, theoretical 2DIR spectra calculated from these MD trajectories allowed us to compare calculated spectra with experimental data, and to determine the diverse structural poses of R5. Our findings indicate that unbound R5 predominantly forms β-strand structures alongside various atypical secondary structures. Post-silicification, there's a noticeable shift: a decrease in β-strands coupled with an increase in turn-type and bend-type configurations. We theorize that this structural transformation stems from silicate embedding within R5's hydrogen-bond network, prompting the peptide backbone to contract and adapt around the biosilica precursors.

Simultaneously Measured Kinetics of Two Amyloid Polymorphs Using Cross Peak Specific 2D IR Spectroscopy.

The origin of in vitro amyloid fibril polymorphs is debated, in part, because few techniques can simultaneously monitor the formation kinetics of multiple amyloid polymorphs. Using a cross-peak specific polarization scheme, ⟨0°,0°,60°,-60°⟩, we resolve 22 previously unseen cross peaks in the 2D IR spectra of amyloid fibrils formed by the human islet amyloid polypeptide (hIAPP). Those cross peaks include a subset assigned to a second fibril polymorph, which forms on a slower time scale. We simulated the data with three different kinetic models for polymorph formation. Only a model based on secondary nucleation reproduces the cross peak kinetics. These experiments are evidence that fibrils formed by secondary nucleation have a different polymorphic structure than the parent fibrils and illustrate the enhanced structural resolution of this new cross peak specific polarization scheme.

New ψ-Santonin Derivatives from Crossostephium chinense and Their Anti-Proliferative Activities against Leishmania major and Human Cancer Cells A549.

Previously, we reported two cytotoxic ψ-santonin-amino acid conjugates isolated from the EtOAc layer of Crossostephium chinense. However, a further phytochemical investigation seems to be required because of the few reports of similar derivatives. In this study, we targeted the 1-BuOH layer, which resulted in the isolation of seven new ψ-santonin derivatives (1-7) together with ten known compounds (8-17). The structures of 1-7 were elucidated based on spectroscopic methods, including 1D and 2D NMR experiments (1H, 13C, DEPT, COSY, HSQC, and HMBC), IR spectrum, and high-resolution electrospray ionization-mass spectrometry (HR-ESI-MS). The stereochemistry of new compounds was confirmed by NOESY and ECD calculations. All isolated compounds were evaluated by in vitro experiments for their anti-proliferative activities against Leishmania major, human lung cancer cell line A549, and Vero cells. As a result, most of the ψ-santonin derivatives, especially 1-5, showed significant cytotoxicity against L. major with a lower IC50 than the positive control we used (miltefosine).

Probing the Hydrophobic Region of a Lipid Bilayer at Specific Depths Using Vibrational Spectroscopy.

A novel spectroscopic approach for studying the flexibility and mobility in the hydrophobic interior of lipid bilayers at specific depths is proposed. A set of test compounds featuring an azido moiety and a cyano or carboxylic acid moiety, connected by an alkyl chain of different lengths, was synthesized. FTIR data and molecular dynamics calculations indicated that the test compounds in a bilayer are oriented so that the cyano or carboxylic acid moiety is located in the lipid head-group region, while the azido group stays inside the bilayer at the depth determined by its alkyl chain length. We found that the asymmetric stretching mode of the azido group (νN3) can serve as a reporter of the membrane interior dynamics. FTIR and two-dimensional infrared (2DIR) studies were performed at different temperatures, ranging from 22 to 45 °C, covering the Lβ-Lα phase transition temperature of dipalmitoylphosphatidylcholine (∼41 °C). The width of the νN3 peak was found to be very sensitive to the phase transition and to the temperature in general. We introduced an order parameter, SN3, which characterizes restrictions to motion inside the bilayer. 2DIR spectra of νN3 showed different extents of inhomogeneity at different depths in the bilayer, with the smallest inhomogeneity in the middle of the leaflet. The spectral diffusion dynamics of the N3 peak was found to be dependent on the depth of the N3 group location in the bilayer. The obtained results enhance our understanding of the bilayer dynamics and can be extended to investigate membranes with more complex compositions.

Investigating Intramolecular H Atom Transfer Dynamics in β-Diketones with Ultrafast Infrared Spectroscopies and Theoretical Modeling.

The vibrational signatures and ultrafast dynamics of the intramolecular H-bond in a series of β-diketones are investigated with 2D IR spectroscopy and computational modeling. The chosen β-diketones exhibit a range of H atom donor-acceptor distances and asymmetry along the H atom transfer coordinate that tunes the intramolecular H-bond strength. The species with the strongest H-bonds are calculated to have very soft H atom potentials, resulting in highly red-shifted OH stretch fundamental frequencies and dislocation of the H atom upon vibrational excitation. These soft potentials lead to significant coupling to the other normal mode coordinates and give rise to the very broad vibrational signatures observed experimentally. The 2D IR spectra in both the OH and OD stretch regions of the light and deuterated isotopologues reveal broadened and long-lived ground-state bleach signatures of the vibrationally hot molecules. Polarization-sensitive transient absorption measurements in the OH and OD stretch regions reveal notable isotopic differences in orientational dynamics. Orientational relaxation was measured to occur on ∼600 fs and ∼2 ps time scales for the light and deuterated isotopologues, respectively. The orientational dynamics are interpreted in terms of activated H/D atom transfer events driven by collective intramolecular structural rearrangements.

Two-Dimensional Infrared Spectroscopy of Isolated Molecular Ions.

Two-dimensional infrared (2D IR) spectroscopy of mass-selected, cryogenically cooled molecular ions is presented. Nonlinear response pathways, encoded in the time-domain photodissociation action response of weakly bound N2 messenger tags, were isolated using pulse shaping techniques following excitation with four collinear ultrafast IR pulses. 2D IR spectra of Re(CO)3(CH3CN)3+ ions capture off-diagonal cross-peak bleach signals between the asymmetric and symmetric carbonyl stretching transitions. These cross peaks display intensity variations as a function of pump-probe delay time due to coherent coupling between the vibrational modes. Well-resolved 2D IR features in the congested fingerprint region of protonated caffeine (C8H10N4O2H+) are also reported. Importantly, intense cross-peak signals were observed at 3 ps waiting time, indicating that tag-loss dynamics are not competing with the measured nonlinear signals. These demonstrations pave the way for more precise studies of molecular interactions and dynamics that are not easily obtainable with current condensed-phase methodologies.

Structural Dependence of Extended Amide III Vibrations in Two-Dimensional Infrared Spectra.

Two-dimensional infrared (2D-IR) spectroscopy is a powerful experimental method for probing the structure and dynamics of proteins in aqueous solution. So far, most experimental studies have focused on the amide I vibrations, for which empirical vibrational exciton models provide a means of interpreting such experiments. However, such models are largely lacking for other regions of the vibrational spectrum. To close this gap, we employ an efficient quantum-chemical methodology for the calculation of 2D-IR spectra, which is based on anharmonic theoretical vibrational spectroscopy with localized modes. We apply this approach to explore the potential of 2D-IR spectroscopy in the extended amide III region. Using calculations for a dipeptide model as well as alanine polypeptides, we show that distinct 2D-IR cross-peaks in the extended amide III region can potentially be used to distinguish α-helix and β-strand structures. We propose that the extended amide III region could be a promising target for future 2D-IR experiments.

Probing Ion Configurations in the KcsA Selectivity Filter with Single-Isotope Labels and 2D IR Spectroscopy.

The potassium ion (K+) configurations of the selectivity filter of the KcsA ion channel protein are investigated with two-dimensional infrared (2D IR) spectroscopy of amide I vibrations. Single 13C-18O isotope labels are used, for the first time, to selectively probe the S1/S2 or S2/S3 binding sites in the selectivity filter. These binding sites have the largest differences in ion occupancy in two competing K+ transport mechanisms: soft-knock and hard-knock. According to the former, water molecules alternate between K+ ions in the selectivity filter while the latter assumes that K+ ions occupy the adjacent sites. Molecular dynamics simulations and computational spectroscopy are employed to interpret experimental 2D IR spectra. We find that in the closed conductive state of the KcsA channel, K+ ions do not occupy adjacent binding sites. The experimental data is consistent with simulated 2D IR spectra of soft-knock ion configurations. In contrast, the simulated spectra for the hard-knock ion configurations do not reproduce the experimental results. 2D IR spectra of the hard-knock mechanism have lower frequencies, homogeneous 2D lineshapes, and multiple peaks. In contrast, ion configurations of the soft-knock model produce 2D IR spectra with a single peak at a higher frequency and inhomogeneous lineshape. We conclude that under equilibrium conditions, in the absence of transmembrane voltage, both water and K+ ions occupy the selectivity filter of the KcsA channel in the closed conductive state. The ion configuration is central to the mechanism of ion transport through potassium channels.

Counteracting Effects of Trimethylamine N-Oxide against Urea in Aqueous Solutions: Insights from Theoretical Two-Dimensional Infrared Spectroscopy.

The study of small osmolytes in their aqueous solutions has gained significant attention because of their relevance to structure and thermodynamics of proteins in aqueous media. Special attention has been given to the binary and ternary aqueous solutions of urea and trimethylamine N-oxide (TMAO). Urea is a well-known protein denaturant, while TMAO protects proteins in their native states. Interestingly, TMAO counteracts urea's ability to denature proteins when present in solutions with approximately half of the concentration of urea. Vibrational spectroscopy can improve our understanding of the molecular origin of this counteracting effect because of its sensitivity to local structure and dynamics. We present results of theoretical linear vibrational and two-dimensional infrared (2DIR) spectroscopy of water in the binary and ternary aqueous solutions of TMAO and urea. The 2DIR spectra are calculated using the electronic structure/molecular dynamics approach. The non-Condon effects in spectral transitions are incorporated in the theoretical calculations of 2DIR spectra. It is found that TMAO disrupts the local structure of water, while urea leaves it essentially unaffected. The 2DIR results show that both TMAO and urea slow down the dynamics of spectral diffusion of water. The extent of slowing down is found to be particularly significant for both hydration and bulk water in the presence of TMAO which can be attributed to strong hydrogen bonds between the water and TMAO molecules. The water molecules present in the hydration layer of the solutes in the ternary solutions are found to relax at even slower rates compared to that in their binary solutions in water. The hydrogen bonds between TMAO and urea are found to be not stable. Thus, the counteracting effect of TMAO against urea is seen to take place mainly through water-mediated interactions. Such TMAO-induced effects giving rise to more structured and slower hydrogen-bonded network are successfully captured through 2DIR spectroscopic calculations.