Applications of 2D IR spectroscopy to peptides, proteins, and hydrogen-bond dynamics

Unveiling the Structural and Dynamic Characteristics of Concentrated LiNO3 Aqueous Solutions through Ultrafast Infrared Spectroscopy and Molecular Dynamics Simulations

Highly concentrated aqueous electrolytes have attracted a significant amount of attention for their potential applications in lithium-ion batteries. Nevertheless, a comprehensive understanding of the Li+ solvation structure and its migration within electrolyte solutions remains elusive. This study employs linear vibrational spectroscopy, ultrafast infrared spectroscopy, and molecular dynamics (MD) simulations to elucidate the structural dynamics in LiNO3 solutions by using intrinsic and extrinsic vibrational probes. The N-O stretching vibrations of NO3- exhibit a distinct spectral splitting, attributed to its asymmetric interaction with the surrounding solvation structure. Analysis of the vibrational relaxation dynamics of intrinsic and extrinsic probes, in combination with MD simulations, reveals cage-like networks formed through electrostatic interactions between Li+ and NO3-. This microscopic heterogeneity is reflected in the intertwined arrangement of ions and water molecules. Furthermore, both vehicular transport and structural diffusion assisted by solvent rearrangement for Li+ were analyzed, which are closely linked with the bulk concentration.

Viscosity effects on the dynamics of diols and diol-based deep eutectic solvents

Diols, characterized by the presence of two hydroxyl groups, form extended hydrogen-bonded networks. Increasing hydrocarbon chain length is known to elevate the viscosity of diols. Given the established influence of viscosity on solvent dynamics, it becomes imperative to comprehend the impact of viscosity on the fluctuation dynamics within diols and establish connections with hydrogen bond formation and breaking dynamics. In this study, we employ two-dimensional infrared spectroscopy to investigate the viscosity dependence of the structural evolution dynamics in three diols with varying chain lengths. Complementing our experimental approach, molecular dynamics simulations are conducted to extract hydrogen bond lifetimes. Our findings reveal a linear correlation between bulk viscosity, solvent fluctuation timescales, and hydrogen bond lifetimes. Notably, the selected diols exhibit the capability to form deep eutectic solvents upon mixing with choline chloride at specific molar ratios. In contrast to molecular solvents like diols, deep eutectic solvents are characterized by the formation of heterogeneous nanodomains, comprising various intercomponent hydrogen-bonded networks. Interestingly, our observations indicate that while the fluctuation dynamics decelerate with increasing bulk viscosity in diol-based deep eutectic solvents, the relationship between viscosity and dynamics is not linear, in contrast to the observed linearity in diols. This nuanced understanding contributes to the broader comprehension of the interplay between viscosity and dynamics in both molecular and deep eutectic solvents.

Chasing weakly-bound biological water in aqueous environment near the peptide backbone by ultrafast 2D infrared spectroscopy

There has been a long-standing debate as to how many hydrogen bonds a peptide backbone amide can form in aqueous solution. Hydrogen-bonding structural dynamics of N-ethylpropionamide (a β-peptide model) in water was examined using infrared (IR) spectroscopy. Two amide-I sub bands arise mainly from amide C=O group that forms strong H-bonds with solvent water molecules (SHB state), and minorly from that involving one weak H-bond with water (WHB state). This picture is supported by molecular dynamics simulations and ab-initio calculations. Further, thermodynamics and kinetics of the SHB and WHB species were examined mainly by chemical-exchange two-dimensional IR spectroscopy, yielding an activation energy for the SHB-to-WHB exchange of 13.25 ± 0.52 kJ mol‒1, which occurs in half picosecond at room temperature. Our results provided experimental evidence of an unstable water molecule near peptide backbone, allowing us to gain more insights into the dynamics of the protein backbone hydration.

Intermolecular Protein-Water Coupling Impedes the Coupling Between the Amide A and Amide I Mode in Interfacial Proteins

The coupling between different vibrational modes in proteins is essential for chemical dynamics and biological functions and is linked to the propagation of conformational changes and pathways of allosteric communication. However, little is known about the influence of intermolecular protein-H2O coupling on the vibrational coupling between amide A (NH) and amide I (C═O) bands. Here, we investigate the NH/CO coupling strength in various peptides with different secondary structures at the lipid cell membrane/H2O interface using femtosecond time-resolved sum frequency generation vibrational spectroscopy (SFG-VS) in which a femtosecond infrared pump is used to excite the amide A band, and SFG-VS is used to probe transient spectral evolution in the amide A and amide I bands. Our results reveal that the NH/CO coupling strength strongly depends on the bandwidth of the amide I mode and the coupling of proteins with water molecules. A large extent of protein-water coupling significantly reduces the delocalization of the amide I mode along the peptide chain and impedes the NH/CO coupling strength. A large NH/CO coupling strength is found to show a strong correlation with the high energy transfer rate found in the light-harvesting proteins of green sulfur bacteria, which may understand the mechanism of energy transfer through a molecular system and assist in controlling vibrational energy transfer by engineering the molecular structures to achieve high energy transfer efficiency.

Critical Evaluation of Polarizable and Nonpolarizable Force Fields for Proteins Using Experimentally Derived Nitrile Electric Fields

Molecular dynamics (MD) simulations are frequently carried out for proteins to investigate the role of electrostatics in their biological function. The choice of force field (FF) can significantly alter the MD results, as the simulated local electrostatic interactions lack benchmarking in the absence of appropriate experimental methods. We recently reported that the transition dipole moment (TDM) of the popular nitrile vibrational probe varies linearly with the environmental electric field, overcoming well-known hydrogen bonding (H-bonding) issues for the nitrile frequency and, thus, enabling the unambiguous measurement of electric fields in proteins (J. Am. Chem. Soc. 2022, 144 (17), 7562-7567). Herein, we utilize this new strategy to enable comparisons of experimental and simulated electric fields in protein environments. Specifically, previously determined TDM electric fields exerted onto nitrile-containing o-cyanophenylalanine residues in photoactive yellow protein are compared with MD electric fields from the fixed-charge AMBER FF and the polarizable AMOEBA FF. We observe that the electric field distributions for H-bonding nitriles are substantially affected by the choice of FF. As such, AMBER underestimates electric fields for nitriles experiencing moderate field strengths; in contrast, AMOEBA robustly recapitulates the TDM electric fields. The FF dependence of the electric fields can be partly explained by the presence of additional negative charge density along the nitrile bond axis in AMOEBA, which is due to the inclusion of higher-order multipole parameters; this, in turn, begets more head-on nitrile H-bonds. We conclude by discussing the implications of the FF dependence for the simulation of nitriles and proteins in general.

Using 2D-IR Spectroscopy to Measure the Structure, Dynamics, and Intermolecular Interactions of Proteins in H2O

Infrared (IR) spectroscopy probes molecular structure at the level of the chemical bond or functional group. In the case of proteins, the most informative band in the IR spectrum is the amide I band, which arises predominantly from the C═O stretching vibration of the peptide link. The folding of proteins into secondary and tertiary structures leads to vibrational coupling between peptide units, generating specific amide I spectral signatures that provide a fingerprint of the macromolecular conformation. Ultrafast two-dimensional IR (2D-IR) spectroscopy allows the amide I band of a protein to be spread over a second frequency dimension in a way that mirrors 2D-NMR methods. This means that amide I 2D-IR spectroscopy produces a spectral map that is exquisitely sensitive to protein structure and dynamics and so provides detailed insights that cannot be matched by IR absorption spectroscopy. As a result, 2D-IR spectroscopy has emerged as a powerful tool for probing protein structure and dynamics over a broad range of time and length scales in the solution phase at room temperature. However, the protein amide I band coincides with an IR absorption from the bending vibration of water (δHOH), the natural biological solvent. To circumvent this problem, protein IR studies are routinely performed in D2O solutions because H/D substitution shifts the solvent bending mode (δDOD) to a lower frequency, revealing the amide I band. While effective, this method raises fundamental questions regarding the impact of the change in solvent mass on the structural or solvation dynamics of the protein and the removal of the energetic resonance between solvent and solute.In this Account, a series of studies applying 2D-IR to study the spectroscopy and dynamics of proteins in H2O-rich solvents is reviewed. A comparison of IR absorption spectroscopy and 2D-IR spectroscopy of protein-containing fluids is used to demonstrate the basis of the approach before a series of applications is presented. These range from measurements of fundamental protein biophysics to recent applications of machine learning to gain insight into protein-drug binding in complex mixtures. An outlook is presented, considering the potential for 2D-IR measurements to contribute to our understanding of protein behavior under near-physiological conditions, along with an evaluation of the obstacles that still need to be overcome.

Dynamics of Formamide-Water Mixtures Investigated by Linear and Nonlinear Infrared Spectroscopy

Formamide (FA) exhibits complete miscibility with water, offering a simplified model for exploring the solvation dynamics of peptide linkages in biophysical processes. Its liquid state demonstrates a three-dimensional hydrogen bonding network akin to water, reflecting solvent-like behavior. Analyzing the microscopic structure and dynamics of FA-water mixtures is expected to provide crucial insights into hydrogen bonding dynamics─a key aspect of various biophysical phenomena. This study is focused on the dynamics of FA-water mixtures using linear and femtosecond infrared spectroscopies. By using the intrinsic OD stretch and extrinsic probe SCN-, the local vibrational behaviors across various FA-water compositions were systematically investigated. The vibrational relaxation of OD stretch revealed a negligible impact of FA addition on the vibrational lifetime of water molecules, underscoring the mixture's water-like behavior. However, the reorientational dynamics of OD stretch slowed with increasing FA mole fraction (XFA), plateauing beyond XFA > 0.5. This suggests a correlation between OD's reorientational time and the strength of the hydrogen bond network, likely tied to the solution's changing dielectric constant. Conversely, the vibrational relaxation dynamics of SCN- was strongly correlated with XFA, highlighting a competition between water and FA molecules in solvating SCN-. Moreover, a linear relationship between rising viscosity and the prolonged correlation time of SCN-'s slow dynamics indicates that the solution's macroscopic viscosity is dictated by the extended structures formed between FA and water molecules. The relation between the reorientation dynamics of the SCN- and the macroscopic viscosity in aqueous FA-water mixture solutions was analyzed by using the Stokes-Einstein-Debye equations. The direct viscosity-diffusion coupling is observed, which can be attributed to the homogeneous dynamics feature in FA-water mixture solutions. The inclusion of these intrinsic and extrinsic probes not only enhances the comprehensiveness of our analysis but also provides valuable insights into various aspects of the dynamics within the FA-water system. This investigation sheds light on the fundamental dynamics of FA-water mixtures, emphasizing their molecular-level homogeneity in this binary mixture solution.

Biomolecular infrared spectroscopy: making time for dynamics

Time resolved infrared spectroscopy of biological molecules has provided a wealth of information relating to structural dynamics, conformational changes, solvation and intermolecular interactions. Challenges still exist however arising from the wide range of timescales over which biological processes occur, stretching from picoseconds to minutes or hours. Experimental methods are often limited by vibrational lifetimes of probe groups, which are typically on the order of picoseconds, while measuring an evolving system continuously over some 18 orders of magnitude in time presents a raft of technological hurdles. In this Perspective, a series of recent advances which allow biological molecules and processes to be studied over an increasing range of timescales, while maintaining ultrafast time resolution, will be reviewed, showing that the potential for real-time observation of biomolecular function draws ever closer, while offering a new set of challenges to be overcome.

Theoretical Study of the Two-Dimensional Vibrational Sum Frequency Generation Spectroscopy of the Air-Water Interface at Varying Temperature and Its Connections to the Interfacial Structure and Dynamics

We performed a theoretical study of the temperature variation of two-dimensional vibrational sum frequency generation (2D-VSFG) spectra of the OH stretch modes at air-water interfaces in the mid-IR region. The calculations are performed at four different temperatures from 250 to 325 K by using a combination of techniques involving response function formalism of nonlinear spectroscopy, electronic structure calculations, and molecular dynamics simulations. Also, the calculations are performed for isotopically dilute solutions so that the intra- and intermolecular coupling between the vibrational modes of interest can be ignored. We have established the connections of temperature variation of various frequency- and time-dependent features of the calculated spectra to the changes in the underlying structure and dynamics of the interfaces. The results reveal that interfacial water is dynamically more heterogeneous than bulk water, with three dominant dynamical processes exhibiting their corresponding time-dependent features in the 2D-VSFG spectrum. These are the spectral diffusion of hydrogen-bonded OH groups at the interface, conversion of an initially hydrogen-bonded OH group to a dangling OH which is a stable state for surface water, unlike the bulk water, and the third one, which involves the conversion of an initially free or dangling OH group to its hydrogen-bonded state at the interface. The temporal appearance of the cross peaks corresponding to interconversion of the hydrogen-bonded state to the dangling state or vice versa of an interfacial OH group is found to take place at a slower rate than the dynamics of spectral diffusion of hydrogen-bonded molecules at the interface, which, in turn, is slower than the corresponding spectral diffusion of bulk water molecules. The temperature variation of these dynamic processes can be linked to the decay of appropriate hydrogen-bond and non-hydrogen-bond time correlation functions of interfacial water molecules for the different air-water systems studied in this work.

Ultrafast Spectroscopy Reveals Slow Water Dynamics in Biocondensates

Cells achieve high spatiotemporal control over biochemical processes through compartmentalization to membrane-bound as well as membraneless organelles that assemble by liquid-liquid phase separation. Characterizing the balance of forces within these environments is essential to understanding their stability and function, and water is an integral part of the condensate, playing an important role in mediating electrostatic and hydrogen-bonding interactions. Here, we investigate the ultrafast, picosecond hydrogen-bond dynamics of a model biocondensate consisting of a peptide poly-l-arginine (Poly-R) and the nucleic acid adenosine monophosphate (AMP) using coherent two-dimensional infrared (2D IR) spectroscopy. We investigated three vibrational modes: the arginine side-chain C═N stretches, an AMP ring mode, and the amide backbone carbonyl stretching modes. Dynamics slow considerably between the dilute phase and the condensate phase for each vibrational probe. For example, the arginine side-chain C═N modes slow from 0.38 to 2.26 ps due to strong electrostatic interactions. All-atom molecular dynamics simulations provide an atomistic interpretation of the H-bond network disruption resulting from electrostatic contributions as well as collapse within the condensate. Simulations predict that a fraction of water molecules are highly constrained within the condensate, explaining the observed slowdown in the H-bond dynamics.

Liquid-Liquid Phase Separation of the Intrinsically Disordered Domain of the Fused in Sarcoma Protein Results in Substantial Slowing of Hydration Dynamics

Formation of liquid condensates plays a critical role in biology via localization of different components or via altered hydrodynamic transport, yet the hydrogen-bonding environment within condensates, pivotal for solvation, has remained elusive. We explore the hydrogen-bond dynamics within condensates formed by the low-complexity domain of the fused in sarcoma protein. Probing the hydrogen-bond dynamics sensed by condensate proteins using two-dimensional infrared spectroscopy of the protein amide I vibrations, we find that frequency-frequency correlations of the amide I vibration decay on a picosecond time scale. Interestingly, these dynamics are markedly slower for proteins in the condensate than in a homogeneous protein solution, indicative of different hydration dynamics. All-atom molecular dynamics simulations confirm that lifetimes of hydrogen-bonds between water and the protein are longer in the condensates than in the protein in solution. Altered hydrogen-bonding dynamics may contribute to unique solvation and reaction dynamics in such condensates.

Crowding alters F-actin secondary structure and hydration

Actin, an important component of eukaryotic cell cytoskeleton, regulates cell shape and transport. The morphology and biochemical properties of actin filaments are determined by their structure and protein-protein contacts. Crowded environments can organize filaments into bundles, but less is known about how they affect F-actin structure. This study used 2D IR spectroscopy and spectral calculations to examine how crowding and bundling impact the secondary structure and local environments in filaments and weakly or strongly bundled networks. The results reveal that bundling induces changes in actin's secondary structure, leading to a decrease in β-sheet and an increase in loop conformations. Strongly bundled networks exhibit a decrease in backbone solvent exposure, with less perturbed α-helices and nearly "locked" β-sheets. Similarly, the loops become less hydrated but maintain a dynamic environment. These findings highlight the role of loop structure in actin network morphology and stability under morphology control by PEG.

Mulberry Leaves (Morus Rubra)-Derived Blue-Emissive Carbon Dots Fed to Silkworms to Produce Augmented Silk Applicable for the Ratiometric Detection of Dopamine

Silk fibers (SF) reeled from silkworms are constituted by natural proteins, and their characteristic structural features render them applicable as materials for textiles and packaging. Modification of SF with functional materials can facilitate their applications in additional areas. In this work, the preparation of functional SF embedded with carbon dots (CD) is reported through the direct feeding of a CD-modified diet to silkworms. Fluorescent and mechanically robust SF are obtained from silkworms (Bombyx mori) that are fed on CDs synthesized from the Morus rubra variant of mulberry leaves (MB-CDs). MB-CDs are introduced to silkworms from the third instar by spraying them on the silkworm feed, the mulberry leaves. MB-CDs are synthesized hydrothermally without adding surface passivating agents and are observed to have a quantum yield of 22%. With sizes of ≈4 nm, MB-CDs exhibited blue fluorescence, and they can be used as efficient fluorophores to detect Dopamine (DA) up to the limit of 4.39 nM. The nanostructures and physical characteristics of SF weren't altered when the SF are infused with MB-CDs. Also, a novel DA sensing application based on fluorescence with the MB-CD incorporated SF is demonstrated.

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.

Chain-like Structures Facilitate Li+ Transport in Concentrated Aqueous Electrolytes: Insights from Ultrafast Infrared Spectroscopy and Molecular Dynamics Simulations

Highly concentrated aqueous electrolytes have attracted attention due to their unique applications in lithium ion batteries (LIBs). However, the solvation structure and transport mechanism of Li+ cations at concentrated concentrations remain largely unexplored. To address this gap in knowledge, we employ ultrafast infrared spectroscopy and molecular dynamics (MD) simulations to reveal the dynamic and spatial structural heterogeneity in aqueous lithium chloride (LiCl) solutions. The coupling between the reorientation dynamics of the extrinsic probe and the macroscopic viscosity in aqueous LiCl solutions was analyzed using the Stokes-Einstein-Debye (SED) equations. MD simulations reveal that the Cl- and Li+ form chain-like structures through electrostatic interactions, supporting the vehicular migration of Li+ through the chain-like structure. The concentration dependent conductivity of the LiCl solution is well reproduced, where Li(H2O)2+ and Li(H2O)3+ are the dominant species that contribute to the conduction of Li+. This study is expected to establish correlations between ion pair structures and macroscopic properties.

Theoretical Two Dimensional Infrared Spectroscopy of Aqueous Solutions of tert-Butyl Alcohol: Variation of the Dynamics of Spectral Diffusion along the Percolation Transition

Binary mixtures of water and tert-butyl alcohol (TBA) are known to exhibit the so-called percolation transition where small clusters of TBA molecules span into large aggregates beyond a threshold concentration of the alcohol. In the present study, we have investigated the linear and two-dimensional infrared spectral features of aqueous solutions of TBA for varying concentration of the alcohol along the percolation transition. The percolation transition is characterized through calculations of intermolecular radial distribution functions and average size of the largest cluster of TBA molecules. It is found that, with variation of alcohol concentration, the radial distribution functions of the central carbon atoms of TBA molecules show a nonmonotonic change in the height of the first peak and also the size of the largest cluster of TBA molecules show a jump in the increase of its size for TBA mole fraction between 0.04 and 0.06 corresponding to a transition from smaller clusters to larger spanning aggregates. However, it is found that the linear infrared spectrum of water does not exhibit any noticeable changes on variation of TBA concentration along the percolation transition. Subsequently, two-dimensional infrared (2DIR) spectra and vibrational frequency time correlation function of water are calculated for all the TBA-water solutions considered in this study. The spectral diffusion of water calculated from 2DIR is found to slow down with increase of the TBA concentration. The time scales of spectral diffusion of water, as characterized by the relaxation of frequency time correlation function, 2DIR metric of central line slope, and also the hydrogen bond time correlation functions, are found to exhibit a noticeable jump along the percolation transition. The hydrophilic group of TBA is found to retard the water dynamics more effectively than the hydrophobic groups. Also, the jump in the dynamical slowdown along the percolation transition is found to be more significant for water molecules at the hydrophilic sites.

Probing local changes to α-helical structures with 2D IR spectroscopy and isotope labeling

α-Helical secondary structures impart specific mechanical and physiochemical properties to peptides and proteins, enabling them to perform a vast array of molecular tasks ranging from membrane insertion to molecular allostery. Loss of α-helical content in specific regions can inhibit native protein function or induce new, potentially toxic, biological activities. Thus, identifying specific residues that exhibit loss or gain of helicity is critical for understanding the molecular basis of function. Two-dimensional infrared (2D IR) spectroscopy coupled with isotope labeling is capable of capturing detailed structural changes in polypeptides. Yet, questions remain regarding the inherent sensitivity of isotope-labeled modes to local changes in α-helicity, such as terminal fraying; the origin of spectral shifts (hydrogen-bonding versus vibrational coupling); and the ability to definitively detect coupled isotopic signals in the presence of overlapping side chains. Here, we address each of these points individually by characterizing a short, model α-helix (DPAEAAKAAAGR-NH2) with 2D IR and isotope labeling. These results demonstrate that pairs of 13C18O probes placed three residues apart can detect subtle structural changes and variations along the length of the model peptide as the α-helicity is systematically tuned. Comparison of singly and doubly labeled peptides affirm that frequency shifts arise primarily from hydrogen-bonding, while vibrational coupling between paired isotopes leads to increased peak areas that can be clearly differentiated from underlying side-chain modes or uncoupled isotope labels not participating in helical structures. These results demonstrate that 2D IR in tandem with i,i+3 isotope-labeling schemes can capture residue-specific molecular interactions within a single turn of an α-helix.

Forced Interactions: Ionic Polymers at Charged Surfactant Interfaces

Characterizing electrostatic interactions at heterogeneous interfaces is critical for developing a fundamental description of the dynamic processes at charged interfaces. Water-in-oil reverse micelles (RMs) offer a high degree of tunability across composition, polarity, and temperature, making them ideal systems for studying interactions at heterogeneous liquid-liquid interfaces. In the present study, we use a combination of ultrafast two-dimensional infrared spectroscopy and molecular dynamics (MD) simulations to determine the picosecond interfacial dynamics in RMs containing binary compositions of sorbitan monostearate and anionic or cationic cosurfactants, which are used to tune the ratio of charged to nonionic surfactants at the interface. The positively charged polyethylenimine (PEI) polymer is encapsulated within the RMs, and the carbonyl stretching mode of sorbitan monostearate reports on the interfacial hydrogen-bond populations and dynamics. The results show that hydrogen-bond populations are altered through the inclusion of both negatively and positively charged cosurfactants. Charged surfactants increase interfacial water penetration into the surfactant layer, and the surface localization of polymers decreases water penetration. Local hydrogen-bond dynamics undergo a slowdown with the inclusion of charged surfactants, and the encapsulation of polymers results in similar effects, irrespective of the charge.

Determining the impact of gold nanoparticles on amyloid aggregation with 2D IR spectroscopy

As nanomaterials become more prevalent in both industry and medicine, it is crucial to fully understand their health risks. One area of concern is the interaction of nanoparticles with proteins, including their ability to modulate the uncontrolled aggregation of amyloid proteins associated with diseases, such as Alzheimer's disease and type II diabetes, and potentially extend the lifetime of cytotoxic soluble oligomers. This work demonstrates that two-dimensional infrared spectroscopy and ¹³C¹⁸O isotope labeling can be used to follow the aggregation of human islet amyloid polypeptide (hIAPP) in the presence of gold nanoparticles (AuNPs) with single-residue structural resolution. 60 nm AuNPs were found to inhibit hIAPP, tripling the aggregation time. Furthermore, calculating the actual transition dipole strength of the backbone amide I' mode reveals that hIAPP forms a more ordered aggregate structure in the presence of AuNPs. Ultimately, such studies can provide insight into how mechanisms of amyloid aggregation are altered in the presence of nanoparticles, furthering our understanding of protein-nanoparticle interactions.

Dicyanamide Anion Reports on Water Induced Local Structural and Dynamic Heterogeneity in Ionic Liquid Mixtures

The effects of limited amounts (under 21.6% χ_Water) of water on 1-butyl-3-methylimidazolium tetrafluoroborate (BmimBF₄) and 1-butyl-3-methylimidazolium dicyanamide (BmimDCA) room-temperature ionic liquid (RTIL) mixtures were characterized by tracking changes in the linear and two-dimensional infrared (2D IR) vibrational features of the dicyanamide anion (DCA). Peak shifts with increasing water suggest the formation of water-associated and nonwater-associated DCA populations. Further results showed clear differences in the dynamic behavior of these different populations of DCA at low (defined here as below 2.5% χ_Water), mid (defined here as between 2.5% χ_Water and 9.6% χ_Water), and high (defined here as between 11.6% χ_Water and 21.6% χ_Water) range water concentrations. Vibrational relaxation is accelerated with increasing water content for water-associated populations of DCA, indicating water facilitates population relaxation, possibly through the provision of additional bath modes. Conversely, spectral diffusion of water-associated populations slowed dramatically with increasing water, suggesting that water drives the formation of distinct and noninterchangeable or very slowly interchangeable local solvent environments.

Ultrafast Two-Dimensional Infrared Spectroscopy Resolved a Structured Lysine 159 on the Cytoplasmic Surface of the Microbial Photoreceptor Bacteriorhodopsin

Bacteriorhodopsin (bR) is a light-driven microbial receptor, and lysine 159 (K159) is a charged residue on the cytoplasmic (CP) side of its E-F loop. However, its conformation and function remain unknown due to fast surface dynamics. By utilizing a ¹³C, ¹⁵N-labeled lysine (K) as an isotope probe, we created a network of site-specific amide-I vibrational signatures (backbone carbonyl stretch) to identify the frequency contribution of the labeled residues to the amide-I excitonic band structure. Thus, the red-shifted amide-I frequency in the ¹³C, ¹⁵N-lysine-labeled bR (uK-bR) to the unlabeled bR (WT-bR) could be differentiated and examined by ultrafast two-dimensional vibrational echo infrared (2D IR) spectroscopy. Our results showed that the backbone carbonyl of K159 is located at a high frequency of ca. 1693 cm^-1 and has a vibrational excited-state relaxation time shorter than the bulk helical amide-I mode at the same frequency, suggesting that K159 may possess a hydrogen-bonded γ-turn structure with E161, one of the carboxylate residues on the CP surface of bR. The 2D solid-state NMR study of uK-bR also revealed conformational dependent lysine residues, from which K159 was found to involve the turn motif. This γ-turn structure maintained by K159 may help to stabilize the E-F loop and support E161 in attracting protons from the bulk during the late stage of the bR photocycle. The combined spectroscopic approach illustrated in this work may be applied to map residue-specific local structures and dynamics of other receptors and large proteins.

Distinct Hydrogen Bonding Dynamics Underlies the Microheterogeneity in DMF-Water Mixtures

The hydrogen bonding interaction between the amide functional group and water is fundamental to understanding the liquid-liquid heterogeneity in biological systems. Herein, the structure and dynamics of the N,N-dimethylformamide (DMF)-water mixtures have been investigated by linear and nonlinear IR spectroscopies, using the hydroxyl stretch and extrinsic probe of thiocyanate as local vibrational reporters. According to vibrational relaxation dynamics measurements, the orientational dynamics of water is not directly tied to those of DMF molecules. Wobbling-in-a-cone analysis demonstrates that the water molecules have varying degrees of angular restriction depending on their composition due to the formation of specific water-DMF networks. Because of the preferential solvation by DMF molecules, the rotational dynamics of the extrinsic probe is slowed significantly, and its rotational time constants are correlated to the change of solution viscosity. The unique structural dynamics observed in the DMF-water mixtures is expected to provide important insights into the underlying mechanism of microscopic heterogeneity in binary mixtures.

Effect of the N-Alkyl Side Chain on the Amide-Water Interactions

Amide-water interactions influence the structure and functions of amide-based systems, such as proteins and homopolymers. In particular, the N-alkylation of the amide unit appears to play a critical role in defining the interactions of the amide group. Previous studies have linked the thermal behavior of amide-based polymers to the nature of their N-alkyl side chain. However, the connection between the chemical structure of the N-alkyl and the hydration of the amide remains elusive. In this study, the solvation structure and dynamics of amides, having differing N-alkyl groups, are investigated using a combination of linear and nonlinear infrared spectroscopies and computational methods. Interestingly, the dynamics of the amide local environment do not slow down as the N-alkyl side chain becomes bulkier, but rather speeds up. Computational calculations confirm the hydration dynamics and assign the effect to smaller amplitude and faster rotations of the bulkier group. It is also observed experimentally that the hydrogen-bond making and breaking between water and the amide carbonyl do not directly relate to the size of the N-alkyl side chain. The bulkier N-isopropyl substituent presents significantly slower chemical exchange dynamics than smaller chains (ethyl and methyl), but the two small groups do not present a major difference. The hydrogen-bond making and breaking disparities and similarities among groups are well modeled by the theory demonstrating that the N-alkyl group affects the amide hydration structure and dynamics via a steric effect. In summary, the results presented here show that the size of the N-substituted alkyl group significantly influences the hydration dynamics of amides and stress the importance of considering this effect on much larger systems, such as polymers.

Atomistic Simulations for Reactions and Vibrational Spectroscopy in the Era of Machine Learning─Quo Vadis?

Atomistic simulations using accurate energy functions can provide molecular-level insight into functional motions of molecules in the gas and in the condensed phase. This Perspective delineates the present status of the field from the efforts of others and some of our own work and discusses open questions and future prospects. The combination of physics-based long-range representations using multipolar charge distributions and kernel representations for the bonded interactions is shown to provide realistic models for the exploration of the infrared spectroscopy of molecules in solution. For reactions, empirical models connecting dedicated energy functions for the reactant and product states allow statistically meaningful sampling of conformational space whereas machine-learned energy functions are superior in accuracy. The future combination of physics-based models with machine-learning techniques and integration into all-purpose molecular simulation software provides a unique opportunity to bring such dynamics simulations closer to reality.

Investigating the effects of N-terminal acetylation on KFE8 self-assembly with 2D IR spectroscopy

Peptide self-assembly is an exciting and robust approach to create novel nanoscale materials for biomedical applications. However, the complex interplay between intra- and intermolecular interactions in peptide aggregation means that minor changes in peptide sequence can yield dramatic changes in supramolecular structure. Here, we use two-dimensional infrared (2D IR) spectroscopy to study a model amphiphilic peptide, KFE8, and its N-terminal acetylated counterpart, AcKFE8. 2D IR spectra of isotope-labeled peptides reveal that AcKFE8 aggregates comprise two distinct β-sheet structures while KFE8 aggregates comprise only one of these structures. Using an excitonic Hamiltonian to simulate the vibrational spectra of model β-sheets, we determine that the spectra are consistent with antiparallel β-sheets with different strand alignments, specifically a two-residue shift in the register of the β-strands. These findings bring forth new insights into how N-terminal acetylation may subtly impact secondary structure, leading to larger effects on overall aggregate morphology. Additionally, these results highlight the importance of understanding the residue-level structural differences that result from changes in peptide sequence in order to facilitate the rational design of peptide materials.

The polarization dependence of 2D IR cross-peaks distinguishes parallel-stranded and antiparallel-stranded DNA G-quadruplexes

Guanine-rich nucleic acid sequences have a tendency to form four-stranded non-canonical motifs known as G-quadruplexes. These motifs may adopt a wide range of structures characterized by size, strand orientation, guanine base conformation, and fold topology. Using three K⁺-bound model systems, we show that vibrational coupling between guanine C6 = O and ring modes varies between parallel-stranded and antiparallel-stranded G-quadruplexes, and that such structures can be distinguished by comparison of the polarization dependences of cross-peaks in their two-dimensional infrared (2D IR) spectra. Combined with previously defined vibrational frequency trends, this analysis reveals key features of a 30-nucleotide unimolecular variant of the Bcl-2 proximal promoter that are consistent with its reported structure. This study shows that 2D IR spectroscopy is a convenient method for analyzing G-quadruplex structures that can be applied to complex sequences where traditional high-resolution methods are limited by solubility and disorder.

Cross-α/β polymorphism of PSMα3 fibrils

Phenol soluble modulins (PSMs) are an important class of peptides secreted by Staphylococcus aureus bacteria. The toxicity to human cells and unique ability of one such peptide, PSMα3, to aggregate into an α-helical amyloid-like structure may hold a key to a better understanding of the virulence of dangerous pathogens such as methicillin resistant S. aureus. In reporting a detailed two-dimensional infrared (2DIR) analysis of PSMα3, we found direct evidence of multiple aggregate architectures existing in equilibrium with one another. We also discovered a unique and characteristic 2DIR spectroscopic signature that unambiguously reports on the presence of the unusual and highly cytotoxic cross-α amyloid structure.

The formation of ordered cross-β amyloid protein aggregates is associated with a variety of human disorders. While conventional infrared methods serve as sensitive reporters of the presence of these amyloids, the recently discovered amyloid secondary structure of cross-α fibrils presents new questions and challenges. Herein, we report results using Fourier transform infrared spectroscopy and two-dimensional infrared spectroscopy to monitor the aggregation of one such cross-α–forming peptide, phenol soluble modulin alpha 3 (PSMα3). Phenol soluble modulins (PSMs) are involved in the formation and stabilization of Staphylococcus aureus biofilms, making sensitive methods of detecting and characterizing these fibrils a pressing need. Our experimental data coupled with spectroscopic simulations reveals the simultaneous presence of cross-α and cross-β polymorphs within samples of PSMα3 fibrils. We also report a new spectroscopic feature indicative of cross-α fibrils.

Exploring the Structure and Complexation Dynamics of Azide Anion Recognition by Calix[4]pyrroles in Solution

The structure and anion recognition dynamics between calix[4]pyrroles and azide (N₃^-) anions in the form of its TBA⁺ and Na⁺ salts were investigated in dimethyl sulfoxide solutions by Fourier transform infrared (FTIR) spectroscopy and ultrafast IR spectroscopy. Vibrational energy redistribution of the N₃^- anion in the complex is accelerated through hydrogen bonding interactions with the N-H proton of the receptor. Rotational dynamics of the bound N₃^- is greatly restricted, demonstrating a distinct countercation effect. The detailed binding modes of N₃^- with the receptor were further evaluated by the density functional theoretical (DFT) calculations and nuclear magnetic resonance (NMR) spectroscopy. All of these measurements support the notion that the calix[4]pyrroles are capable of capturing the azide anion in solution. However, the calix[4]pyrroles may not necessarily undergo a conformational change to a cone-like geometry when they bind to the azide anion in the solution.

Cationic Effects on the Structural Dynamics of the Metal Ion-Crown Ether Complexes Investigated by Ultrafast Infrared Spectroscopy

It is usually believed that the binding affinity and selectivity of an alkali metal ion with crown ether are defined by the size matching model. However, the underlying mechanism of the specific host-guest interactions and the structural dynamics of the metal ions confined in the cavity of the crown ethers in the solutions are still not clear. In this report, a series of alkali thiocyanate salts (XSCN; X = Li, Na, K, and Cs) complexed with 18-crown-6 (a typical crown ether) in the chloroform solutions were studied by the polarization-selective infrared pump-probe spectroscopy and the ultrafast two-dimensional infrared (2D IR) spectroscopy. The SCN^- counteranions were employed as the local vibrational probe to reveal the specific host-guest interactions in the crown ether complexes. The rotational dynamics and spectral diffusion of SCN^- vibration were both measured by ultrafast IR spectroscopy, and it was found that the metal cations hosted by the crown ethers can have a pronounced effect on the rotational dynamics of the counteranions. The reorientational time constants of the SCN^- vibration in the complexation follow the order Li⁺ > Na⁺ > K⁺ ≃ Cs⁺. More importantly, the spectral diffusion dynamics of SCN^-, which quantifies the decay of the correlation of the frequency fluctuations in the complexation, was also affected by the metal ions but showed a different order of cationic effect. A detailed analysis of the 2D IR data showed that the spectral diffusion of SCN^- counteranion clearly decayed with two different time scales in the complex of 18-crown-6 with K⁺. The 3-4-fold slowdown in spectral diffusion indicated that the fluctuation of SCN^- vibrational transition frequency was strongly affected by the K⁺ cation due to the geometric constraint imposed by the crown ether. The results should help the researchers to unravel the specific host-guest interactions and further reveal the origination of the binding selectivity of crown ether for metal cations in the condensed phases from the perspective of structural dynamics.

Two-dimensional Infrared Spectroscopy Reveals Better Insights of Structure and Dynamics of Protein

Proteins play an important role in biological and biochemical processes taking place in the living system. To uncover these fundamental processes of the living system, it is an absolutely necessary task to understand the structure and dynamics of the protein. Vibrational spectroscopy is an established tool to explore protein structure and dynamics. In particular, two-dimensional infrared (2DIR) spectroscopy has already proven its versatility to explore the protein structure and its ultrafast dynamics, and it has essentially unprecedented time resolutions to observe the vibrational dynamics of the protein. Providing several examples from our theoretical and experimental efforts, it is established here that two-dimensional vibrational spectroscopy provides exceptionally more information than one-dimensional vibrational spectroscopy. The structural information of the protein is encoded in the position, shape, and strength of the peak in 2DIR spectra. The time evolution of the 2DIR spectra allows for the visualisation of molecular motions.

Long Vibrational Lifetime R-Selenocyanate Probes for Ultrafast Infrared Spectroscopy: Properties and Synthesis

Ultrafast infrared vibrational spectroscopy is widely used for the investigation of dynamics in systems from water to model membranes. Because the experimental observation window is limited to a few times the probe's vibrational lifetime, a frequent obstacle for the measurement of a broad time range is short molecular vibrational lifetimes (typically a few to tens of picoseconds). Five new long-lifetime aromatic selenocyanate vibrational probes have been synthesized and their vibrational properties characterized. These probes are compared to commercial phenyl selenocyanate. The vibrational lifetimes range between ∼400 and 500 ps in complex solvents, which are some of the longest room-temperature vibrational lifetimes reported to date. In contrast to vibrations that are long-lived in simple solvents such as CCl₄, but become much shorter in complex solvents, the probes discussed here have ∼400 ps lifetimes in complex solvents and even longer in simple solvents. One of them has a remarkable lifetime of 1235 ps in CCl₄. These probes have a range of molecular sizes and geometries that can make them useful for placement into different complex materials due to steric reasons, and some of them have functionalities that enable their synthetic incorporation into larger molecules, such as industrial polymers. We investigated the effect of a range of electron-donating and electron-withdrawing para-substituents on the vibrational properties of the CN stretch. The probes have a solvent-independent linear relationship to the Hammett substituent parameter when evaluated with respect to the CN vibrational frequency and the ipso ¹³C NMR chemical shift.

Transparent window 2D IR spectroscopy of proteins

Proteins are complex, heterogeneous macromolecules that exist as ensembles of interconverting states on a complex energy landscape. A complete, molecular-level understanding of their function requires experimental tools to characterize them with high spatial and temporal precision. Infrared (IR) spectroscopy has an inherently fast time scale that can capture all states and their dynamics with, in principle, bond-specific spatial resolution. Two-dimensional (2D) IR methods that provide richer information are becoming more routine but remain challenging to apply to proteins. Spectral congestion typically prevents selective investigation of native vibrations; however, the problem can be overcome by site-specific introduction of amino acid side chains that have vibrational groups with frequencies in the "transparent window" of protein spectra. This Perspective provides an overview of the history and recent progress in the development of transparent window 2D IR of proteins.

Protein Dynamics by Two-Dimensional Infrared Spectroscopy

Proteins function as ensembles of interconverting structures. The motions span from picosecond bond rotations to millisecond and longer subunit displacements. Characterization of functional dynamics on all spatial and temporal scales remains challenging experimentally. Two-dimensional infrared spectroscopy (2D IR) is maturing as a powerful approach for investigating proteins and their dynamics. We outline the advantages of IR spectroscopy, describe 2D IR and the information it provides, and introduce vibrational groups for protein analysis. We highlight example studies that illustrate the power and versatility of 2D IR for characterizing protein dynamics and conclude with a brief discussion of the outlook for biomolecular 2D IR.

Substituent Effects on the Vibrational Properties of the CN Stretch Mode of Aromatic Nitriles: IR Probes Useful for Time-resolved IR Spectroscopy

Developing ideal IR probes is essential to understand the structure and dynamics of biomolecules with time-resolved IR spectroscopies and imaging techniques. Especially, nitrile (CN) group has recently been proposed to serve as IR probes of the local environment of proteins. Herein, we investigated the effect of a substituent on the vibrational properties of the benzonitrile. The electron-donating and withdrawing character of p-substituent on benzonitrile are expected to modulate the vibrational frequency, molar extinction coefficient, and vibrational lifetime of CN probe. FT-IR revealed the positive correlation between electron-donating character and the molar extinction coefficient of CN stretch mode. Infrared pump-probe (IR-PP) measurements showed that the vibrational lifetime of CN stretch mode exhibits a relatively weak correlation with the electron-donating strength. Among the investigated samples, 4-dimethylamino benzonitrile with the strongest electron-donating strength shows enhanced absorption and extended vibrational lifetime. Utilizing substituent effects will be a practical strategy to improve the performance of the IR probe.

A computational insight into the relationship between side chain IR line shapes and local environment in fibril-like structures

Infrared spectroscopy is a widely used technique to characterize protein structures and protein mediated processes. While the amide I band provides information on proteins' secondary structure, amino acid side chains are used as infrared probes for the investigation of protein reactions and local properties. In this paper, we use a hybrid quantum mechanical/classical molecular dynamical approach based on the perturbed matrix method to compute the infrared band due to the C=O stretching mode of amide-containing side chains. We calculate, at first, the infrared band of zwitterionic glutamine in water and obtain results in very good agreement with the experimental data. Then, we compute the signal arising from glutamine side chains in a microcrystal of the yeast prion Sup35-derived peptide, GNNQQNY, with a fibrillar structure. The infrared bands obtained by selective isotopic labeling of the two glutamine residues, Q4 and Q5, of each peptide were experimentally used to investigate the local hydration in the fibrillar microcrystal. The experimental spectra of the two glutamine residues, which experience different hydration environments, feature different spectral signals that are well reproduced by the corresponding calculated spectra. In addition, the analysis of the simulated spectra clarifies the molecular origin of the experimentally observed spectroscopic differences that arise from the different local electric field experienced by the two glutamine residues, which is, in turn, determined by a different hydrogen bonding pattern.

Theory of two-dimensional spectroscopy with intense laser fields

Two-dimensional vibrational and electronic spectroscopic observables of isotropically oriented molecular samples in solution are sensitive to laser field intensities and polarization. The third-order response function formalism predicts a signal that grows linearly with the field strength of each laser pulse, thus lacking a way of accounting for non-trivial intensity-dependent effects, such as saturation and finite bleaching. An analytical expression to describe the orientational part of the molecular response, which, in the weak-field limit, becomes equivalent to a four-point correlation function, is presented. Such an expression is evaluated for Liouville-space pathways accounting for diagonal and cross peaks for all-parallel and cross-polarized pulse sequences, in both the weak- and strong-field conditions, via truncation of a Taylor series expansion at different orders. The results obtained in the strong-field conditions suggest how a careful analysis of two-dimensional spectroscopic experimental data should include laser pulse intensity considerations when determining molecular internal coordinates.

Elucidating the mechanism behind the infrared spectral features and dynamics observed in the carbonyl stretch region of organic carbonates interacting with lithium ions

Ultrafast infrared spectroscopy has become a very important tool for studying the structure and ultrafast dynamics in solution. In particular, it has been recently applied to investigate the molecular interactions and motions of lithium salts in organic carbonates. However, there has been a discrepancy in the molecular interpretation of the spectral features and dynamics derived from these spectroscopies. Hence, the mechanism behind spectral features appearing in the carbonyl stretching region was further investigated using linear and nonlinear spectroscopic tools and the co-solvent dilution strategy. Lithium perchlorate in a binary mixture of dimethyl carbonate (DMC) and tetrahydrofuran was used as part of the dilution strategy to identify the changes of the spectral features with the number of carbonates in the first solvation shell since both solvents have similar interaction energetics with the lithium ion. Experiments showed that more than one carbonate is always participating in the lithium ion solvation structures, even at the low concentration of DMC. Moreover, temperature-dependent study revealed that the exchange of the solvent molecules coordinating the lithium ion is not thermally accessible at room temperature. Furthermore, time-resolved IR experiments confirmed the presence of vibrationally coupled carbonyl stretches among coordinated DMC molecules and demonstrated that this process is significantly altered by limiting the number of carbonate molecules in the lithium ion solvation shell. Overall, the presented experimental findings strongly support the vibrational energy transfer as the mechanism behind the off-diagonal features appearing on the 2DIR spectra of solutions of lithium salt in organic carbonates.

CH Mode Mixing Determines the Band Shape of the Carboxylate Symmetric Stretch in Apo-EDTA, Ca2+-EDTA, and Mg2+-EDTA

The infrared spectra of EDTA complexed with Ca²⁺ and Mg²⁺ contain, to date, unidentified vibrational bands. This study assigns the peaks in the linear and two-dimensional infrared spectra of EDTA, with and without either Ca²⁺ or Mg²⁺ ions. Two-dimensional infrared spectroscopy and DFT calculations reveal that, in both the presence and absence of ions, the carboxylate symmetric stretch and the terminal CH bending vibrations mix. We introduce a method to calculate participation coefficients that quantify the contribution of the carboxylate symmetric stretch, CH wag, CH twist, and CH scissor in the 1400-1550 cm^-1 region. With the help of participation coefficients, we assign the 1400-1430 cm^-1 region to the carboxylate symmetric stretch, which can mix with CH modes. We assign the 1000-1380 cm^-1 region to CH twist modes, the 1380-1430 cm^-1 region to wag modes, and the 1420-1650 cm^-1 region to scissor modes. The difference in binding geometry between the carboxylate-Ca²⁺ and carboxylate-Mg²⁺ complex manifests as new diagonal and cross-peaks between the mixed modes in the two complexes. The small Mg²⁺ ion binds EDTA tighter than the Ca²⁺ ion, which causes a redshift of the COO symmetric stretches of the sagittal carboxylates. Energy decomposition analysis further characterizes the importance of electrostatics and deformation energy in the bound complexes.

Bursting the bubble: A molecular understanding of surfactant-water interfaces

Surfactant science has historically emphasized bulk, thermodynamic measurements to understand the microemulsion properties of greatest industrial significance, such as interfacial tensions, phase behavior, and thermal stability. Recently, interest in the molecular properties of surfactants has grown among the physical chemistry community. This has led to the application of cutting-edge spectroscopic methods and advanced simulations to understand the specific interactions that give rise to the previously studied bulk characteristics. In this Perspective, we catalog key findings that describe the surfactant-oil and surfactant-water interfaces in molecular detail. We emphasize the role of ultrafast spectroscopic methods, including two-dimensional infrared spectroscopy and sum-frequency-generation spectroscopy, in conjunction with molecular dynamics simulations, and the role these techniques have played in advancing our understanding of interfacial properties in surfactant microemulsions.

Lessons from combined experimental and theoretical examination of the FTIR and 2D-IR spectroelectrochemistry of the amide I region of cytochrome c

Amide I difference spectroscopy is widely used to investigate protein function and structure changes. In this article, we show that the common approach of assigning features in amide I difference signals to distinct secondary structure elements in many cases may not be justified. Evidence comes from Fourier transform infrared (FTIR) and 2D-IR spectroelectrochemistry of the protein cytochrome c in the amide I range, in combination with computational spectroscopy based on molecular dynamics (MD) simulations. This combination reveals that each secondary structure unit, such as an alpha-helix or a beta-sheet, exhibits broad overlapping contributions, usually spanning a large part of the amide I region, which in the case of difference absorption experiments (such as in FTIR spectroelectrochemistry) may lead to intensity-compensating and even sign-changing contributions. We use cytochrome c as the test case, as this small electron-transferring redox-active protein contains different kinds of secondary structure units. Upon switching its redox-state, the protein exhibits a different charge distribution while largely retaining its structural scaffold. Our theoretical analysis suggests that the change in charge distribution contributes to the spectral changes and that structural changes are small. However, in order to confidently interpret FTIR amide I difference signals in cytochrome c and proteins in general, MD simulations in combination with additional experimental approaches such as isotope labeling, the insertion of infrared labels to selectively probe local structural elements will be required. In case these data are not available, a critical assessment of previous interpretations of protein amide I 1D- and 2D-IR difference spectroscopy data is warranted.

Membrane-dependent amyloid aggregation of human BAX α9 (173-192)

Mitochondrial outer membrane permeabilization, which is a critical step in apoptosis, is initiated upon transmembrane insertion of the C-terminal α-helix (α9) of the proapoptotic Bcl-2 family protein BAX. The isolated α9 fragment (residues 173-192) is also competent to disrupt model membranes, and the structures of its membrane-associated oligomers are of interest in understanding the potential roles of this sequence in apoptosis. Here, we used ultrafast two-dimensional infrared (2D IR) spectroscopy, thioflavin T binding, and transmission electron microscopy to show that the synthetic BAX α9 peptide (α9p) forms amyloid aggregates in aqueous environments and on the surfaces of anionic small unilamellar vesicles. Its inherent amyloidogenicity was predicted by sequence analysis, and 2D IR spectra reveal that vesicles modulate the β-sheet structures of insoluble aggregates, motivating further examination of the formation or suppression of BAX amyloids in apoptosis.

Detection of Glycine as a Model Protein in Blood Serum Using 2D-IR Spectroscopy

Glycine (Gly) is used as a model system to evaluate the ability of ultrafast two-dimensional infrared (2D-IR) spectroscopy to detect and quantify the low-molecular-weight proteinaceous components of blood serum. Combining data acquisition schemes to suppress absorption bands of H₂O that overlap with the protein amide I band with analysis of peak patterns appearing in the off-diagonal region of the 2D-IR spectrum allows separation of the Gly spectral signature from that of the dominant protein fraction of serum in a transmission-mode 2D-IR measurement without any sample manipulation, e.g., filtration or drying. 2D-IR spectra of blood serum samples supplemented with varying concentrations of Gly were obtained, and a range of data analysis methods compared, leading to a detection limit of ∼3 mg/mL for Gly. The reported methodology provides a platform for a critical assessment of the sensitivity of 2D-IR for measuring the concentrations of amino acids, peptides, and low-molecular-weight proteins present in serum samples. We conclude that, in the case of several clinically relevant diagnostic molecules and their combinations, the potential exists for 2D-IR to complement IR absorption methods as the benefits of the second frequency dimension offered by 2D-IR spectroscopy outweigh the added technical complexity of the measurement.

Bottom-Up Approach to Assess the Molecular Structure of Aqueous Poly(N-Isopropylacrylamide) at Room Temperature via Infrared Spectroscopy

The structure of poly(N-isopropylacrylamide) (PNIPAM) in solution is still an unresolved topic. Here, the PNIPAM structure in water was investigated using a bottom-up approach, involving the monomer, dimer, and trimer, and a combination of infrared (IR) spectroscopies as well as molecular dynamics simulations. The experiments show that the monomer and oligomers exhibit a broad and asymmetric amide I band with two underlying transitions, while PNIPAM presents the same major transitions and a minor one. Analysis of the 2D IR spectra and theoretical modeling of the amide I band indicates that the two transitions of the monomer do not have the same molecular origin as the oligomers and the polymer. In the monomer, the two bands originate from the ultrafast rotation of its ethyl group, which leads to different solvation structures for the various rotational conformers. In the case of the oligomers, the asymmetry and splitting of the amide I band is caused by the vibrational coupling among adjacent amide side chains. Moreover, it is deduced from the simulations that the oligomers have three distinct backbone conformations for neighboring amides. In particular, two of the backbone conformations have a closed and compact structure, while in the third, the backbone is open and elongated. The bottom-up approach allowed us to infer that such backbone conformations exist in PNIPAM as well. Consequently, the two major amide I transitions of the polymer are also assigned to split amide I transitions resulting from the vibrationally coupled nearest-neighboring amides. In contrast, the additional minor transition observed in PNIPAM is assigned to unsolvated amide units of the polymer. The proposed molecular model successfully describes that PNIPAM amide I band changes with temperature in terms of its molecular structure. This new model strongly suggests that PNIPAM does not have a completely random backbone structure, but has distinct backbone conformers between neighboring amides.

Specific Host-Guest Interactions in the Crown Ether Complexes with K+ and NH4+ Revealed from the Vibrational Relaxation Dynamics of the Counteranion

The specific host-guest interactions in the corresponding complexes of K⁺ and NH₄⁺ with typical crown ethers were investigated by using FTIR and ultrafast IR spectroscopies. The counteranions, i.e., SCN^-, were employed as a local vibrational probe to report the structural dynamics of the complexation. It was found that the vibrational relaxation dynamics of the SCN^- was strongly affected by the cations confined in the cavities of the crown ethers. The time constant of the vibrational population decay of SCN^- in the complex of NH₄⁺ with the 18-crown-6 was determined to be 6 ± 2 ps, which is ∼30 times faster than that in the complex of K⁺ with the crown ethers. Control experiments showed that the vibrational population decay of SCN^- depended on the size of the cavities of the crown ethers. A theoretical calculation further indicated that the nitrogen atom of SCN^- showed preferential coordination to the K⁺ ions hosted by the crown ethers, while the NH₄⁺ can form hydrogen bonds with the oxygen atoms in the studied crown ethers. The geometric constraints formed in the complex of crown ethers can cause a specific interaction between the NH₄⁺ and SCN^-, which can facilitate the intermolecular vibrational energy redistribution of the SCN^-.

Temperature-Jump 2D IR Spectroscopy with Intensity-Modulated CW Optical Heating

Pulsed temperature-jump (T-jump) spectroscopy with infrared (IR) detection has been widely used to study biophysical processes occurring from nanoseconds to ∼1 ms with structural sensitivity. However, many systems exhibit structural dynamics on time scales longer than the millisecond barrier that is set by the time scale for thermal relaxation of the sample. We developed a linear and nonlinear infrared spectrometer coupled to an intensity-modulated continuous wave (CW) laser to probe T-jump-initiated chemical reactions from <1 ms to seconds. Time-dependent modulation of the CW laser leads to a <1 ms heating time as well as a constant final temperature (±3%) for the duration of the heating time. Temperature changes of up to 75 °C in D₂O are demonstrated, allowing for nonequilibrium measurements inaccessible to standard pulsed optical T-jump setups. T-jump linear absorption, pump-probe, and two-dimensional IR (2D IR) spectroscopy are applied to the unfolding and refolding of ubiquitin and a model intercalated motif (i-motif) DNA sequence, and analysis of the observed signals is used to demonstrate the limits and utility of each method. Overall, the ability to probe temperature-induced chemical processes from <1 ms to many seconds with 2D IR spectroscopy provides multiple new avenues for time-dependent spectroscopy in chemistry and biophysics.

Extending the vibrational lifetime of azides with heavy atoms

The development of novel vibrational reporters (VRs), aka infrared (IR) probes, to study local environments and dynamic processes in biomolecules and materials continues to be an important area of research. Azides are important VRs because of their small size and large transition dipole strengths, however, their relatively short vibrational lifetimes (<2 ps) have limited their full potential. Herein we report that the vibrational lifetimes of azides can be increased by attaching them to heavy atoms and by using heavy 15N isotopes. Three group 14 atom triphenyl azides (Ph3CN3, Ph3SiN3, Ph3SnN3), and their triple-15N isotopomers, were synthesized in good yields. Tributyltin azide and its heavy isotopomer (Bu3Sn15N3) were also prepared to probe the effect of molecular scaffolding. The extinction coefficients for the natural abundance azides were determined, ranging from 900 to 1500 M-1 cm-1. The vibrational lifetimes of all azides were measured by pump-probe IR spectroscopy and each showed a major component with a short-to-moderate vibrational lifetime and a minor component with a much longer vibrational lifetime. Based on these results, the lifetime, aka the observation window, of an azide reporter can be extended from ∼2 ps to as long as ∼300 ps by a combination of isotopic labeling and heavy atom effect. 2D IR measurements of these compounds further confirmed the ability to observe these azide transitions at much longer timescales showing their utility to capture dynamic processes from tens to hundreds of picoseconds.

Infrared Spectroscopy Reveals the Preferred Motif Size and Local Disorder in Parallel Stranded DNA G-Quadruplexes

Infrared spectroscopy detects the formation of G-quadruplexes in guanine-rich nucleic acid sequences through shifts in the guanine C=O stretch mode. Here, we use ultrafast 2D infrared (IR) spectroscopy and isotope substitution to show that these shifts arise from vibrational delocalization among stacked G-quartets. This provides a direct measure of the sizes of locally ordered motifs in heterogeneous samples with substantial disordered regions. We find that parallel-stranded, potassium-bound DNA G-quadruplexes are limited to five consecutive G-quartets and 3-4 consecutive layers are preferred for longer polyguanine tracts. The resulting potassium-dependent G-quadruplex assembly landscape reflects the polyguanine tract lengths found in genomes, the ionic conditions prevalent in healthy mammalian cells, and the onset of structural disorder in disease states. Our study describes spectral markers that can be used to probe other G-quadruplex structures and provides insight into the fundamental limits of their formation in biological and artificial systems.