Highly sensitive and label-free protein immunoassay-based biosensor comprising infrared metamaterial absorber inducing strong coupling

A mid-infrared label-free immunoassay-based biosensor is an effective device to help identify and quantify biomolecules. This biosensor employs a surface-enhanced infrared absorption spectroscopy, which is a highly potent sensing technique for detecting minute quantities of analytes. In this study, a biosensor was constructed using a metamaterial absorber, which facilitated strong coupling effects. For maximum coupling effect, it is necessary to enhance the near-field intensity and the spatial and spectral overlap between the optical cavity resonance and the vibrational mode of the analyte. Due to significant peak splitting, conventional baseline correction methods fail to adequately analyze such a coupling system. Therefore, we employed a coupled harmonic oscillation model to analyze the spectral distortion resulting from the peak splitting induced by the strong coupling effect. The proposed biosensor with a thrombin-binding aptamer-based immunoassay could achieve a limit of detection of 267.4 pM, paving the way for more efficient protein detection in clinical practice.

Stimuli-Responsive Oligourea Molecular Films

We have designed and synthesized a helical cysteamine-terminated oligourea foldamer composed of ten urea residues featuring side carboxyl and amine groups. The carboxyl group is located in proximity to the C-terminus of the oligourea and hence at the negative pole of the helix dipole. The amine group is located close to the N-terminus and hence at the positive pole of the helix dipole. Beyond the already remarkable dipole moment inherent in oligourea 2.5 helices, the incorporation of additional charges originating from the carboxylic and amine groups is supposed to impact the overall charge distribution along the molecule. These molecules were self-assembled into monolayers on a gold substrate, allowing us to investigate the influence of an electric field on these polar helices. By applying surface-enhanced infrared reflection-absorption spectroscopy, we proved that molecules within the monolayers tend to reorient themselves more vertically when a negative bias is applied to the surface. It was also found that surface-confined oligourea molecules affected by the external electric field tend to rearrange the electron density at urea groups, leading to the stabilization of the resonance structure with charge transfer character. The presence of the external electric field also affected the nanomechanical properties of the oligourea films, suggesting that molecules also tend to reorient in the ambient environment without an electrolyte solution. Under the same conditions, the helical oligourea displayed a robust piezoresponse, particularly noteworthy given the slender thickness of the monolayer, which measured approximately 1.2 nm. This observation demonstrates that thin molecular films composed of oligoureas may exhibit stimulus-responsive properties. This, in turn, may be used in nanotechnology systems as actuators or functional films, enabling precise control of their thickness in the range of even fractions of nanometers.

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.

Increased accuracy of vibrational circular dichroism calculations for isotopically labeled helical peptides

Vibrational circular dichroism (VCD) spectra have been computed with qualitatively correct sign patterns for α-helical peptides using various methods, ranging from empirical models to ab initio quantum mechanical computations. However, some details, such as deuteration effects and isotope substitution shifts and sign patterns for the resultant amide I' band shape, have remained a predictive challenge. Fully optimized computations for a 25-residue Ala-rich peptide, including implicit solvent corrections and explicit side chains that experimentally stabilize these model helical peptides in water, have been carried out using density functional theory (DFT). These fully minimized structures show minor changes in the (ϕ,ψ) torsions at the termini and yield an extra negative band to the low energy side of the characteristic amide I' couplet VCD, in agreement with experiments. Additionally, these calculations give the right sign and relative intensity patterns, as compared to experimental results, for several 13C=O substituted variants. The differences from previously reported computations that used ideal helical structures and vacuum conditions imply that inclusion of distorted termini and solvent effects can have an impact on the final detailed spectral patterns. Inclusion of side chains in these calculations had very little effect on the computed amide I' IR and VCD. Tests of constrained geometries, varying dielectric, and different functionals indicate that each can affect the band shapes, particularly for the 12C=O components, but these aspects do not fully explain the difference from previous spectral simulations. Inclusion of long-range amide coupling, as obtained from DFT computation of the full structure, or transfer of parameters from a somewhat longer peptide model, rather than shorter model, seems to be more important for the final detailed band shape under isotopic substitution. However, these corrections can also induce other changes, suggesting that previously reported, limited calculations may have been qualitatively useful due to a balance of errors. This may also explain the success of simple empirical IR models.

Lipids and proteins: Insights into the dynamics of assembly, recognition, condensate formation. What is still missing?

Lipid membranes and proteins, which are part of us throughout our lives, have been studied for decades. However, every year, new discoveries show how little we know about them. In a reader-friendly manner for people not involved in the field, this paper tries to serve as a bridge between physicists and biologists and new young researchers diving into the field to show its relevance, pointing out just some of the plethora of lines of research yet to be unraveled. It illustrates how new ways, from experimental to theoretical approaches, are needed in order to understand the structures and interactions that take place in a single lipid, protein, or multicomponent system, as we are still only scratching the surface.

Water inside the Selectivity Filter of a K+ Ion Channel: Structural Heterogeneity, Picosecond Dynamics, and Hydrogen Bonding

Water inside biological ion channels regulates the key properties of these proteins, such as selectivity, ion conductance, and gating. In this article, we measure the picosecond spectral diffusion of amide I vibrations of an isotope-labeled KcsA potassium channel using two-dimensional infrared (2D IR) spectroscopy. By combining waiting time (100-2000 fs) 2D IR measurements of the KcsA channel including 13C18O isotope-labeled Val76 and Gly77 residues with molecular dynamics simulations, we elucidated the site-specific dynamics of water and K+ ions inside the selectivity filter of KcsA. We observe inhomogeneous 2D line shapes with extremely slow spectral diffusion. Our simulations quantitatively reproduce the experiments and show that water is the only component with any appreciable dynamics, whereas K+ ions and the protein are essentially static on a picosecond timescale. By analyzing simulated and experimental vibrational frequencies, we find that water in the selectivity filter can be oriented to form hydrogen bonds with adjacent or nonadjacent carbonyl groups with the reorientation timescales being three times slower and comparable to that of water molecules in liquid, respectively. Water molecules can reside in the cavity sufficiently far from carbonyls and behave essentially like "free" gas-phase-like water with fast reorientation times. Remarkably, no interconversion between these configurations was observed on a picosecond timescale. These dynamics are in stark contrast with liquid water, which remains highly dynamic even in the presence of ions at high concentrations.

Serum-based ATR-FTIR spectroscopy combined with multivariate analysis for the diagnosis of pre-diabetes and diabetes

Diabetes mellitus (DM) is a metabolic disease with an increasing prevalence that is causing worldwide concern. The pre-diabetes stage is the only reversible stage in the patho-physiological process towards DM. Due to the limitations of traditional methods, the diagnosis and detection of DM and pre-diabetes are complicated, expensive, and time-consuming. Therefore, it would be of great benefit to develop a simple, rapid and inexpensive diagnostic test. Herein, the infrared (IR) spectra of serum samples from 111 DM patients, 111 pre-diabetes patients and 333 healthy volunteers were collected using attenuated total reflection Fourier-transform IR (ATR-FTIR) spectroscopy and this was combined with the multivariate analysis of principal component analysis linear discriminant analysis (PCA-LDA) to develop a discriminant model to verify the diagnostic potential of this approach. The study found that the accuracy of the test model established by ATR-FTIR spectroscopy combined with PCA-LDA was 97%, and the sensitivity and specificity were 100% and 100% in the control group, 94% and 98% in the pre-diabetes group, and 91% and 98% in the DM group, respectively. This indicates that this method can effectively diagnose DM and pre-diabetes, which has far-reaching clinical significance.

Refining protein amide I spectrum simulations with simple yet effective electrostatic models for local wavenumbers and dipole derivative magnitudes

Analysis of the amide I band of proteins is probably the most wide-spread application of bioanalytical infrared spectroscopy. Although highly desirable for a more detailed structural interpretation, a quantitative description of this absorption band is still difficult. This work optimized several electrostatic models with the aim to reproduce the effect of the protein environment on the intrinsic wavenumber of a local amide I oscillator. We considered the main secondary structures - α-helices, parallel and antiparallel β-sheets - with a maximum of 21 amide groups. The models were based on the electric potential and/or the electric field component along the CO bond at up to four atoms in an amide group. They were bench-marked by comparison to Hessian matrices reconstructed from density functional theory calculations at the BPW91, 6-31G** level. The performance of the electrostatic models depended on the charge set used to calculate the electric field and potential. Gromos and DSSP charge sets, used in common force fields, were not optimal for the better performing models. A good compromise between performance and the stability of model parameters was achieved by a model that considered the electric field at the positions of the oxygen, nitrogen, and hydrogen atoms of the considered amide group. The model describes also some aspects of the local conformation effect and performs similar on its own as in combination with an explicit implementation of the local conformation effect. It is better than a combination of a local hydrogen bonding model with the local conformation effect. Even though the short-range hydrogen bonding model performs worse, it captures important aspects of the local wavenumber sensitivity to the molecular surroundings. We improved also the description of the coupling between local amide I oscillators by developing an electrostatic model for the dependency of the dipole derivative magnitude on the protein environment.

Water inside the selectivity filter of a K+ ion channel: structural heterogeneity, picosecond dynamics, and hydrogen-bonding

Water inside biological ion channels regulates the key properties of these proteins such as selectivity, ion conductance, and gating. In this Article we measure the picosecond spectral diffusion of amide I vibrations of an isotope labeled KcsA potassium channel using two-dimensional infrared (2D IR) spectroscopy. By combining waiting time (100 - 2000 fs) 2D IR measurements of the KcsA channel including 13C18O isotope labeled Val76 and Gly77 residues with molecular dynamics simulations, we elucidated the site-specific dynamics of water and K+ ions inside the selectivity filter of KcsA. We observe inhomogeneous 2D lineshapes with extremely slow spectral diffusion. Our simulations quantitatively reproduce the experiments and show that water is the only component with any appreciable dynamics, whereas K+ ions and the protein are essentially static on a picosecond timescale. By analyzing simulated and experimental vibrational frequencies, we find that water in the selectivity filter can be oriented to form hydrogen bonds with adjacent, or non-adjacent carbonyl groups with the reorientation timescales being three times slower and comparable to that of water molecules in liquid, respectively. Water molecules can reside in the cavity sufficiently far from carbonyls and behave essentially like "free" gas-phase-like water with fast reorientation times. Remarkably, no interconversion between these configurations were observed on a picosecond timescale. These dynamics are in stark contrast with liquid water that remains highly dynamic even in the presence of ions at high concentrations.

Hidden Beneath the Layers: Extending the Core/Shell Model of Reverse Micelles

Reverse micelles (RMs) provide a unique and highly tunable model system to study water in confined environments. The complex properties of water within RMs arise from the disruption of extended hydrogen bond (H-bond) networks that mediate local and long-range dynamics in bulk aqueous systems. Modulating the water pool size influences its H-bond dynamics, with smaller RMs increasingly restricting the H-bond network rearrangements leading to slower dynamics; however, within small confined systems, the dynamics of the surfactants also influence the water dynamics. Using ultrafast two-dimensional infrared spectroscopy, we investigate the effects of RM size on the surfactant headgroup rotamer populations and picosecond interfacial H-bond dynamics of aerosol-OT surfactants. We find that the increased water penetration accelerates H-bond dynamics, with larger RMs showing faster dynamics. These results imply that the changes in the RM structure alter the physical structure of the RM interface and thus alter the solvation dynamics. The findings in this study can be used for developing models for structure-specific solvation dynamics that account for the surfactant packing and hydration at the interface.

Unpolarized laser method for infrared spectrum calculation of amide I CO bonds in proteins using molecular dynamics simulation

The investigation of the strong infrared (IR)-active amide I modes of peptides and proteins has received considerable attention because a wealth of detailed information on hydrogen bonding, dipole-dipole interactions, and the conformations of the peptide backbone can be derived from the amide I bands. The interpretation of experimental spectra typically requires substantial theoretical support, such as direct ab-initio molecular dynamics simulation or mixed quantum-classical description. However, considering the difficulties associated with these theoretical methods and their applications are limited in small peptides, it is highly desirable to develop a simple yet efficient approach for simulating the amide I modes of any large proteins in solution. In this work, we proposed a comprehensive computational method that extends the well-established molecular dynamics (MD) simulation method to include an unpolarized IR laser for exciting the CO bonds of proteins. We showed the amide I frequency corresponding to the frequency of the laser pulse which resonated with the CO bond vibration. At this frequency, the protein energy and the CO bond length fluctuation were maximized. Overall, the amide I bands of various single proteins and amyloids agreed well with experimental data. The method has been implemented into the AMBER simulation package, making it widely available to the scientific community. Additionally, the application of the method to simulate the transient amide I bands of amyloid fibrils during the IR laser-induced disassembly process was discussed in details.

Studies of Local DNA Backbone Conformation and Conformational Disorder Using Site-Specific Exciton-Coupled Dimer Probe Spectroscopy

The processes of genome expression, regulation, and repair require direct interactions between proteins and DNA at specific sites located at and near single-stranded-double-stranded DNA (ssDNA-dsDNA) junctions. Here, we review the application of recently developed spectroscopic methods and analyses that combine linear absorbance and circular dichroism spectroscopy with nonlinear 2D fluorescence spectroscopy to study the local conformations and conformational disorder of the sugar-phosphate backbones of ssDNA-dsDNA fork constructs that have been internally labeled with exciton-coupled cyanine (iCy3)2 dimer probes. With the application of these methods, the (iCy3)2 dimer can serve as a reliable probe of the mean local conformations and conformational distributions of the sugar-phosphate backbones of dsDNA at various critical positions. The results of our studies suggest a possible structural framework for understanding the roles of DNA breathing in driving the processes of protein-DNA complex assembly and function.

Equivalence of quantum and classical third order response for weakly anharmonic coupled oscillators

Two-dimensional (2D) infrared (IR) spectra are commonly interpreted using a quantum diagrammatic expansion that describes the changes to the density matrix of quantum systems in response to light-matter interactions. Although classical response functions (based on Newtonian dynamics) have shown promise in computational 2D IR modeling studies, a simple diagrammatic description has so far been lacking. Recently, we introduced a diagrammatic representation for the 2D IR response functions of a single, weakly anharmonic oscillator and showed that the classical and quantum 2D IR response functions for this system are identical. Here, we extend this result to systems with an arbitrary number of bilinearly coupled, weakly anharmonic oscillators. As in the single-oscillator case, quantum and classical response functions are found to be identical in the weakly anharmonic limit or, in experimental terms, when the anharmonicity is small relative to the optical linewidth. The final form of the weakly anharmonic response function is surprisingly simple and offers potential computational advantages for application to large, multi-oscillator systems.

Simulation of two-dimensional infrared Raman spectroscopy with application to proteins

Two-dimensional infrared Raman spectroscopy is a powerful technique for studying the structure and interaction in molecular and biological systems. Here, we present a new implementation of the simulation of the two-dimensional infrared Raman signals. The implementation builds on the numerical integration of the Schrödinger equation approach. It combines the prediction of dynamics from molecular dynamics with a map-based approach for obtaining Hamiltonian trajectories and response function calculations. The new implementation is tested on the amide-I region for two proteins, where one is dominated by α-helices and the other by β-sheets. We find that the predicted spectra agree well with experimental observations. We further find that the two-dimensional infrared Raman spectra at least of the studied proteins are much less sensitive to the laser polarization used compared to conventional two-dimensional infrared experiments. The present implementation and findings pave the way for future applications for the interpretation of two-dimensional infrared Raman spectra.

Direct calculation of the temperature dependence of 2D-IR spectra: Urea in water

A method for directly calculating the temperature derivative of two-dimensional infrared (2D-IR) spectra from simulations at a single temperature is presented. The approach is demonstrated by application to the OD stretching spectrum of isotopically dilute aqueous (HOD in H₂O) solutions of urea as a function of concentration. Urea is an important osmolyte because of its ability to denature proteins, which has motivated significant interest in its effect on the structure and dynamics of water. The present results show that the temperature dependence of both the linear IR and 2D-IR spectra, which report on the underlying energetic driving forces, is more sensitive to urea concentration than the spectra themselves. Additional physical insight is provided by calculation of the contributions to the temperature derivative from different interactions, e.g., water-water, water-urea, and urea-urea, present in the system. Finally, it is demonstrated how 2D-IR spectra at other temperatures can be obtained from only room temperature simulations.

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.

Generative Adversarial Neural Networks for Denoising Coherent Multidimensional Spectra

Ultrafast spectroscopy often involves measuring weak signals and long data acquisition times. Spectra are typically collected as a "pump-probe" spectrum by measuring differences in intensity across laser shots. Shot-to-shot intensity fluctuations are most often the primary source of noise in ultrafast spectroscopy. Here, we present a novel approach for denoising ultrafast two-dimensional infrared (2D IR) spectra using conditional generative adversarial neural networks (cGANNs). The cGANN approach is able to eliminate shot-to-shot noise and reconstruct the line shapes present in the noisy input spectrum. We present a general approach for training the cGANN using matched pairs of noisy and clean synthetic 2D IR spectra based on the Kubo-line shape model for a three-level system. Experimental shot-to-shot laser noise is added to synthetic spectra to recreate the noise profile present in measured experimental spectra. The cGANNs can recover line shapes from synthetic 2D IR spectra with signal-to-noise ratios as low as 2:1, while largely preserving the key features such as center frequencies, line widths, and diagonal elongation. In addition, we benchmark the performance of the cGANN using experimental 2D IR spectra of an ester carbonyl vibrational probe and demonstrate that, by applying the cGANN denoising approach, we can extract the frequency-frequency time correlation function (FFCF) from reconstructed spectra using a nodal-line slope analysis. Finally, we provide a set of practical guidelines for extending the denoising method to other coherent multidimensional spectroscopies.

The effect of hydration and dehydration on the conformation, assembling behavior and photoluminescence of PBLG

Hydration and dehydration play crucial roles in hydrophobic effects (HEs) and are yet to be understood. Poly(γ-benzyl-L-glutamate) (PBLG) homopolymers in THF/water with various water contents were investigated. We discovered that PBLG was hydrated at low water contents and adopted a helical conformation. The chain became dehydrated with increasing water content, which converted the PBLG₁₀₀ helix to a PPII-helix. The variation in the conformation resulted in an alteration of the self-assembled morphologies from fibers to particles. For PBLG₁₂ with a shorter chain, the chain underwent an α-to-β transition in the conformation due to dehydration as the water content increased, and correspondingly the morphologies varied from tapes to helical ribbons, and eventually to toroids at a higher water content. We also observed that this α-to-β transition is accompanied by an increase in intensity of the fluorescence, which is attributed to the through-space-conjugation of tightly packed phenyl groups within the β-sheet. The discovered effect of hydration and dehydration on the PBLG chain conformation, self-assembling behavior and optical function is essential for the innovation of polypeptide materials and understanding of water-mediated biological systems.

AIM: A Mapping Program for Infrared Spectroscopy of Proteins

Here, we present a new analysis program, AIM, that allows extracting the vibrational amide-I Hamiltonian using molecular dynamics trajectories for protein infrared spectroscopy modeling. The constructed Hamiltonians can be used as input for spectral calculations allowing the calculation of infrared absorption spectra, vibrational circular dichroism, and two-dimensional infrared spectra. These spectroscopies allow the study of the structure and dynamics of proteins. We will explain the essence of how AIM works and give examples of the information and spectra that can be obtained with the program using the Trypsin Inhibitor as an example. AIM is freely available from GitHub, and the package contains a demonstration allowing easy introduction to the use of the program.

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

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

Revealing the Relationship between Electric Fields and the Conformation of Oxytocin Using Quasi-Static Amide-I Two-Dimensional Infrared Spectra

It is reported that the cis/trans conformation change of the peptide hormone oxytocin plays an important role in its receptors and activation and the cis conformation does not lead to antagonistic activity. Motivated by recent experiments and theories, the quasi-static amide-I 2D IR spectra of oxytocin are investigated using DFT/B3LYP (D3)/6-31G (d, p) in combination with the isotope labeling method under different electric fields. The theoretical amide-I IR spectra and bond length of the disulfide bond are consistent with the experimental values, which indicates that the theoretical modes are reasonable. Our theoretical results demonstrate that the oxytocin conformation is transformed from the cis conformation to the trans conformation with the change of the direction of the electric field, which is confirmed by the distance of the backbone carbonyl oxygen of Cys6 and Pro7, the Ramachandran plot of Cys6 and Pro7, the dihedral angle of C_β-S-S-C_β, and the rmsd of the oxytocin backbone. Moreover, the trans conformation as the result of the turn in the vicinity of Pro7 has a tighter secondary spatial structure than the cis conformation, including stronger hydrogen bonds, longer γ-turn geometry involving five amino acids, and a more stable disulfide bond. Our work provides new insights into the relationship between the conformation, the activation of the peptide hormone oxytocin, and the electric fields.

Temperature-dependent local conformations and conformational distributions of cyanine dimer labeled single-stranded-double-stranded DNA junctions by 2D fluorescence spectroscopy

DNA replication and the related processes of genome expression require binding, assembly, and function of protein complexes at and near single-stranded (ss)-double-stranded (ds) DNA junctions. These central protein-DNA interactions are likely influenced by thermally induced conformational fluctuations of the DNA scaffold across an unknown distribution of functionally relevant states to provide regulatory proteins access to properly conformed DNA binding sites. Thus, characterizing the nature of conformational fluctuations and the associated structural disorder at ss-dsDNA junctions is critical for understanding the molecular mechanisms of these central biological processes. Here, we describe spectroscopic studies of model ss-dsDNA fork constructs that contain dimers of "internally labeled" cyanine (iCy3) chromophore probes that have been rigidly inserted within the sugar-phosphate backbones of the DNA strands. Our combined analyses of absorbance, circular dichroism, and two-dimensional fluorescence spectroscopy permit us to characterize the local conformational parameters and conformational distributions. We find that the DNA sugar-phosphate backbones undergo abrupt successive changes in their local conformations-initially from a right-handed and ordered DNA state to a disordered splayed-open structure and then to a disordered left-handed conformation-as the dimer probes are moved across the ss-dsDNA junction. Our results suggest that the sugar-phosphate backbones at and near ss-dsDNA junctions adopt specific position-dependent local conformations and exhibit varying extents of conformational disorder that deviate widely from the Watson-Crick structure. We suggest that some of these conformations can function as secondary-structure motifs for interaction with protein complexes that bind to and assemble at these sites.

Progress in Simulation Studies of Insulin Structure and Function

Insulin is a peptide hormone known for chiefly regulating glucose level in blood among several other metabolic processes. Insulin remains the most effective drug for treating diabetes mellitus. Insulin is synthesized in the pancreatic β-cells where it exists in a compact hexameric architecture although its biologically active form is monomeric. Insulin exhibits a sequence of conformational variations during the transition from the hexamer state to its biologically-active monomer state. The structural transitions and the mechanism of action of insulin have been investigated using several experimental and computational methods. This review primarily highlights the contributions of molecular dynamics (MD) simulations in elucidating the atomic-level details of conformational dynamics in insulin, where the structure of the hormone has been probed as a monomer, dimer, and hexamer. The effect of solvent, pH, temperature, and pressure have been probed at the microscopic scale. Given the focus of this review on the structure of the hormone, simulation studies involving interactions between the hormone and its receptor are only briefly highlighted, and studies on other related peptides (e.g., insulin-like growth factors) are not discussed. However, the review highlights conformational dynamics underlying the activities of reported insulin analogs and mimetics. The future prospects for computational methods in developing promising synthetic insulin analogs are also briefly highlighted.

Direct observation of peptide hydrogel self-assembly

The characterization of self-assembling molecules presents significant experimental challenges, especially when associated with phase separation or precipitation. Transparent window infrared (IR) spectroscopy leverages site-specific probes that absorb in the “transparent window” region of the biomolecular IR spectrum. Carbon–deuterium (C–D) bonds are especially compelling transparent window probes since they are non-perturbative, can be readily introduced site selectively into peptides and proteins, and their stretch frequencies are sensitive to changes in the local molecular environment. Importantly, IR spectroscopy can be applied to a wide range of molecular samples regardless of solubility or physical state, making it an ideal technique for addressing the solubility challenges presented by self-assembling molecules. Here, we present the first continuous observation of transparent window probes following stopped-flow initiation. To demonstrate utility in a self-assembling system, we selected the MAX1 peptide hydrogel, a biocompatible material that has significant promise for use in drug delivery and medical applications. C–D labeled valine was synthetically introduced into five distinct positions of the twenty-residue MAX1 β-hairpin peptide. Consistent with current structural models, steady-state IR absorption frequencies and linewidths of C–D bonds at all labeled positions indicate that these side chains occupy a hydrophobic region of the hydrogel and that the motion of side chains located in the middle of the hairpin is more restricted than those located on the hairpin ends. Following a rapid change in ionic strength to initiate self-assembly, the peptide absorption spectra were monitored as function of time, allowing determination of site-specific time constants. We find that within the experimental resolution, MAX1 self-assembly occurs as a cooperative process. These studies suggest that stopped-flow transparent window FTIR can be extended to other time-resolved applications, such as protein folding and enzyme kinetics.

To facilitate the characterization of phase-transitioning molecules, site-specific non-perturbative infrared probes are leveraged for continuous observation of the self-assembly of fibrils in a peptide hydrogel following stopped-flow initiation.

Distinct electrostatic frequency tuning rates for amide I and amide I' vibrations

Amide I spectroscopy probes the backbone C=O stretch vibrations of peptides and proteins. Amide I spectra are often collected in deuterated water (D₂O) since this provides a cleaner background in the amide I frequency range; such data are often referred to as amide I' spectra since deuteration induces changes in the mode structure, including a roughly ∼10 cm^-1 redshift. For biological samples, however, deuteration is often not possible. As amide I frequency maps are increasingly applied to quantitative protein structural analysis, this raises the interesting challenge of drawing direct connections between amide I and amide I' data. We here analyze amide I and amide I' peak frequencies for a series of dipeptides and related compounds. Changes in protonation state induce large electrostatic shifts in the peak frequencies, allowing us to amass a sizable library of data points for direct amide I/amide I' comparison. While we find an excellent linear correlation between amide I and amide I' peak frequencies, the deuteration-induced shift is smaller for more red-shifted vibrations, indicating different electrostatic tuning rates in the two solvents. H₂O/D₂O shifts were negligible for proline-containing dipeptides that lack exchangeable amide hydrogens, indicating that the intrinsic properties of the solvent do not strongly influence the H/D shift. These findings indicate that the distinct tuning rates observed for the two vibrations arise from modifications to the intrinsic properties of the amide bond and provide (at least for solvated dipeptides) a simple, linear "map" for translating between amide I and amide I' frequencies.

Structural Ensemble of the Insulin Monomer

Experimental evidence suggests that monomeric insulin exhibits significant conformational heterogeneity, and modifications of apparently disordered regions affect both biological activity and the longevity of pharmaceutical formulations, presumably through receptor binding and fibrillation/degradation, respectively. However, a microscopic understanding of conformational heterogeneity has been lacking. Here, we integrate all-atom molecular dynamics simulations with an analysis pipeline to investigate the structural ensemble of human insulin monomers. We find that 60% of the structures present at least one of the following elements of disorder: melting of the A-chain N-terminal helix, detachment of the B-chain N-terminus, and detachment of the B-chain C-terminus. We also observe partial melting and extension of the B-chain helix and significant conformational heterogeneity in the region containing the B-chain β-turn. We then estimate hydrogen-exchange protection factors for the sampled ensemble and find them in line with experimental results for KP-insulin, although the simulations underestimate the importance of unfolded states. Our results help explain the ready exchange of specific amide sites that appear to be protected in crystal structures. Finally, we discuss the implications for insulin function and stability.

Accurate prediction of mutation-induced frequency shifts in chlorophyll proteins with a simple electrostatic model

Photosynthetic pigment-protein complexes control local chlorophyll (Chl) transition frequencies through a variety of electrostatic and steric forces. Site-directed mutations can modify this local spectroscopic tuning, providing critical insight into native photosynthetic functions and offering the tantalizing prospect of creating rationally designed Chl proteins with customized optical properties. Unfortunately, at present, no proven methods exist for reliably predicting mutation-induced frequency shifts in advance, limiting the method's utility for quantitative applications. Here, we address this challenge by constructing a series of point mutants in the water-soluble chlorophyll protein of Lepidium virginicum and using them to test the reliability of a simple computational protocol for mutation-induced site energy shifts. The protocol uses molecular dynamics to prepare mutant protein structures and the charge density coupling model of Adolphs et al. [Photosynth. Res. 95, 197-209 (2008)] for site energy prediction; a graphical interface that implements the protocol automatically is published online at http://nanohub.org/tools/pigmenthunter. With the exception of a single outlier (presumably due to unexpected structural changes), we find that the calculated frequency shifts match the experiment remarkably well, with an average error of 1.6 nm over a 9 nm spread in wavelengths. We anticipate that the accuracy of the method can be improved in the future with more advanced sampling of mutant protein structures.

Lanthanide-dependent coordination interactions in lanmodulin: a 2D IR and molecular dynamics simulations study

The biological importance of lanthanides, and the early lanthanides (La³⁺-Nd³⁺) in particular, has only recently been recognized, and the structural principles underlying selective binding of lanthanide ions in biology are not yet well established. Lanmodulin (LanM) is a novel protein that displays unprecedented affinity and selectivity for lanthanides over most other metal ions, with an uncommon preference for the early lanthanides. Its utilization of EF-hand motifs to bind lanthanides, rather than the Ca²⁺ typically recognized by these motifs in other proteins, has led it to be used as a model system to understand selective lanthanide recognition. Two-dimensional infrared (2D IR) spectroscopy combined with molecular dynamics simulations were used to investigate LanM's selectivity mechanisms by characterizing local binding site geometries upon coordination of early and late lanthanides as well as calcium. These studies focused on the protein's uniquely conserved proline residues in the second position of each EF-hand binding loop. We found that these prolines constrain the EF-hands for strong coordination of early lanthanides. Substitution of this proline results in a more flexible binding site to accommodate a larger range of ions but also results in less compact coordination geometries and greater disorder within the binding site. Finally, we identify the conserved glycine in the sixth position of each EF-hand as a mediator of local binding site conformation and global secondary structure. Uncovering fundamental structure-function relationships in LanM informs the development of synthetic biology technologies targeting lanthanides in industrial applications.

Segregation on the nanoscale coupled to liquid water polyamorphism in supercooled aqueous ionic-liquid solution

The most intriguing hypothesis explaining many water anomalies is a metastable liquid-liquid phase transition (LLPT) at high pressure and low temperatures, experimentally hidden by homogeneous nucleation. Recent infrared spectroscopic experiments showed that upon addition of hydrazinium trifluoroacetate to water, the supercooled ionic solution undergoes a sharp, reversible LLPT at ambient pressure, possible offspring of that in pure water. Here, we calculate the temperature-dependent signature of the OH-stretching band, reporting on the low/high density phase of water, in neat water and in the same experimentally investigated ionic solution. The comparison between the infrared signature of the pure liquid and that of the ionic solution can be achieved only computationally, providing insight into the nature of the experimentally observed phase transition and allowing us to investigate the effects of ionic compounds on the high to low density supercooled liquid water transition. We show that the experimentally observed crossover behavior in the ionic solution can be reproduced only if the phase transition between the low- and high-density liquid states of water is coupled to a mixing-unmixing transition between the water component and the ions: at low temperatures, water and ions are separated and the water component is a low density liquid. At high temperatures, water and ions get mixed and the water component is a high-density liquid. The separation at low temperatures into ion-rich and ion-poor regions allows unveiling the polyamorphic nature of liquid water, leading to a crossover behavior resembling that observed in supercooled neat water under high pressure.

Pump Slice Amplitudes: A Simple and Robust Method for Connecting Two-Dimensional Infrared and Fourier Transform Infrared Spectra

Ultrafast two-dimensional infrared (2D IR) spectroscopy and Fourier transform infrared (FTIR) spectroscopy are often performed in tandem, with FTIR typically used to interpret and provide hypotheses for 2D IR experiments. Comparisons between 2D IR and FTIR spectra can also be used to examine the structure and orientation in systems of coupled vibrational chromophores. The most common method for comparing 2D IR and FTIR lineshapes, the diagonal slice method, contains significant artifacts when applied to oscillators with low anharmonicities. Here, we introduce a new technique, the pump slice amplitude (PSA) method, for relating 2D IR lineshapes to FTIR lineshapes and compare PSAs against diagonal slices using theoretical and experimental spectra. We find that PSAs are significantly more similar to FTIR lineshapes than diagonal slices in systems with low anharmonicity.

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.

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.

Computational IR Spectroscopy of Insulin Dimer Structure and Conformational Heterogeneity

We have investigated the structure and conformational dynamics of insulin dimer using a Markov state model (MSM) built from extensive unbiased atomistic molecular dynamics simulations and performed infrared spectral simulations of the insulin MSM to describe how structural variation within the dimer can be experimentally resolved. Our model reveals two significant conformations to the dimer: a dominant native state consistent with other experimental structures of the dimer and a twisted state with a structure that appears to reflect a ∼55° clockwise rotation of the native dimer interface. The twisted state primarily influences the contacts involving the C-terminus of insulin's B chain, shifting the registry of its intermolecular hydrogen bonds and reorganizing its side-chain packing. The MSM kinetics predict that these configurations exchange on a 14 μs time scale, largely passing through two Markov states with a solvated dimer interface. Computational amide I spectroscopy of site-specifically ¹³C¹⁸O labeled amides indicates that the native and twisted conformation can be distinguished through a series of single and dual labels involving the B24F, B25F, and B26Y residues. Additional structural heterogeneity and disorder is observed within the native and twisted states, and amide I spectroscopy can also be used to gain insight into this variation. This study will provide important interpretive tools for IR spectroscopic investigations of insulin structure and transient IR kinetics experiments studying the conformational dynamics of insulin dimer.

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.

Infrared spectroscopy probes ion binding geometries

Infrared (IR) spectroscopy is a well-established technique for probing the structure, behavior, and surroundings of molecules in their native environments. Its characteristics-most specifically high structural sensitivity, ready applicability to aqueous samples, and broad availability-make it a valuable enzymological technique, particularly for the interrogation of ion binding sites. While IR spectroscopy of the "garden variety" (steady state at room temperature with wild-type proteins) is versatile and powerful in its own right, the combination of IR spectroscopy with specialized experimental schemes for leveraging ultrafast time resolution, protein labeling, and other enhancements further extends this utility. This book chapter provides the fundamental physical background and literature context essential for harnessing IR spectroscopy in the general context of enzymology with specific focus on interrogation of ion binding. Studies of lanthanide ions binding to calmodulin are highlighted as illustrative examples of this process. Appropriate sample preparation, data collection, and spectral interpretation are discussed from a detail-oriented and practical perspective with the goal of facilitating the reader's rapid progression from reading words in a book to collecting and analyzing their own data in the lab.

Structure Changes of a Membrane Polypeptide under an Applied Voltage Observed with Surface-Enhanced 2D IR Spectroscopy

The structures of many membrane-bound proteins and polypeptides depend on the membrane potential. However, spectroscopically studying their structures under an applied field is challenging, because a potential is difficult to generate across more than a few bilayers. We study the voltage-dependent structures of the membrane-bound polypeptide, alamethicin, using a spectroelectrochemical cell coated with a rough, gold film to create surface plasmons. The plasmons sufficiently enhance the 2D IR signal to measure a single bilayer. The film is also thick enough to conduct current and thereby apply a potential. The 2D IR spectra resolve features from both 3₁₀- and α-helical structures and cross-peaks connecting the two. We observe changes in the peak intensity, not their frequencies, upon applying a voltage. A similar change occurs with pH, which is known to alter the angle of alamethicin relative to the surface normal. The spectra are modeled using a vibrational exciton Hamiltonian, and the voltage-dependent spectra are consistent with a change in angle of the 3₁₀- and α-helices in the membrane from 55 to 44°and from 31 to 60°, respectively. The 3₁₀- and α-helices are coupled by approximately 10 cm^-1. These experiments provide new structural information about alamethicin under a potential difference and demonstrate a technique that might be applied to voltage-gated membrane proteins and compared to molecular dynamics structures.

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

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

Catalytic Principles from Natural Enzymes and Translational Design Strategies for Synthetic Catalysts

As biocatalysts, enzymes are characterized by their high catalytic efficiency and strong specificity but are relatively fragile by requiring narrow and specific reactive conditions for activity. Synthetic catalysts offer an opportunity for more chemical versatility operating over a wider range of conditions but currently do not reach the remarkable performance of natural enzymes. Here we consider some new design strategies based on the contributions of nonlocal electric fields and thermodynamic fluctuations to both improve the catalytic step and turnover for rate acceleration in arbitrary synthetic catalysts through bioinspired studies of natural enzymes. With a focus on the enzyme as a whole catalytic construct, we illustrate the translational impact of natural enzyme principles to synthetic enzymes, supramolecular capsules, and electrocatalytic surfaces.

Windowless detection geometry for sum frequency scattering spectroscopy in the C-D and amide I regions

Understanding the structure and chemistry of nanoscopic surfaces is an important challenge for biointerface sciences. Sum frequency scattering (SFS) spectroscopy can specifically probe the surfaces of nanoparticles, vesicles, liposomes, and other materials relevant to biomaterial research, and, as a vibrational spectroscopy method, it can provide molecular level information about the surface chemistry. SFS is particularly promising to probe the structure of proteins, and other biological molecules, at nanoparticle surfaces. Here, amide I spectra can provide information about protein folding and orientation, while spectra in the C-D and C-H stretching regions allow experiments to determine the mode of interaction between particle surfaces and proteins. Methods used currently employ a closed liquid cell or cuvette, which works extremely well for C-H and phosphate regions but is often impeded in the amide I and C-D regions by a strong background signal that originates from the window material of the sample cells. Here, we discuss a windowless geometry for collecting background-free and high-fidelity SFS spectra in the amide I and C-D regions. We demonstrate the improvement in spectra quality by comparing SFS spectra of unextruded, multilamellar vesicles in a sample cuvette with those recorded using the windowless geometry. The sample geometry we propose will enable new experiments using SFS as a probe for protein-particle interactions.

A Machine Learning Protocol for Predicting Protein Infrared Spectra

Infrared (IR) absorption provides important chemical fingerprints of biomolecules. Protein secondary structure determination from IR spectra is tedious since its theoretical interpretation requires repeated expensive quantum-mechanical calculations in a fluctuating environment. Herein we present a novel machine learning protocol that uses a few key structural descriptors to rapidly predict amide I IR spectra of various proteins and agrees well with experiment. Its transferability enabled us to distinguish protein secondary structures, probe atomic structure variations with temperature, and monitor protein folding. This approach offers a cost-effective tool to model the relationship between protein spectra and their biological/chemical properties.

Ultrafast Dynamics at Lipid-Water Interfaces

Lipid membranes are more than just barriers between cell compartments; they provide molecular environments with a finely tuned balance between hydrophilic and hydrophobic interactions that enable proteins to dynamically fold and self-assemble to regulate biological function. Characterizing dynamics at the lipid-water interface is essential to understanding molecular complexities from the thermodynamics of liquid-liquid phase separation down to picosecond-scale reorganization of interfacial hydrogen-bond networks.Ultrafast vibrational spectroscopy, including two-dimensional infrared (2D IR) and vibrational sum-frequency generation (VSFG) spectroscopies, is a powerful tool to examine picosecond interfacial dynamics. Two-dimensional IR spectroscopy provides a bond-centered view of dynamics with subpicosecond time resolutions, as vibrational frequencies are highly sensitive to the local environment. Recently, 2D IR spectroscopy has been applied to carbonyl and phosphate vibrations intrinsically located at the lipid-water interface. Interface-specific VSFG spectroscopy probes the water vibrational modes directly, accessing H-bond strength and water organization at lipid headgroup positions. Signals in VSFG arise from the interfacial dipole contributions, directly probing headgroup ordering and water orientation to provide a structural view of the interface.In this Account we discuss novel applications of ultrafast spectroscopy to lipid membranes, a field that has experienced significant growth over the past decade. In particular, ultrafast experiments now offer a molecular perspective on increasingly complex membranes. The powerful combination of ultrafast, interface-selective spectroscopy and simulations opens up new routes to understanding multicomponent membranes and their function. This Account highlights key prevailing views that have emerged from recent experiments: (1) Water dynamics at the lipid-water interface are slow compared to those of bulk water as a result of disrupted H-bond networks near the headgroups. (2) Peptides, ions, osmolytes, and cosolvents perturb interfacial dynamics, indicating that dynamics at the interface are affected by bulk solvent dynamics and vice versa. (3) The interfacial environment is generally dictated by the headgroup structure and orientation, but hydrophobic interactions within the acyl chains also modulate interfacial dynamics. Ultrafast spectroscopy has been essential to characterizing the biophysical chemistry of the lipid-water interface; however, challenges remain in interpreting congested spectra as well as designing appropriate model systems to capture the complexity of a membrane environment.

Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction

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

Two-dimensional Raman spectroscopy of Lennard-Jones liquids via ring-polymer molecular dynamics

The simulation of multidimensional vibrational spectroscopy of condensed-phase systems including nuclear quantum effects is challenging since full quantum-mechanical calculations are still intractable for large systems comprising many degrees of freedom. Here, we apply the recently developed double Kubo transform (DKT) methodology in combination with ring-polymer molecular dynamics (RPMD) for evaluating multi-time correlation functions [K. A. Jung et al., J. Chem. Phys. 148, 244105 (2018)], providing a practical method for incorporating nuclear quantum effects in nonlinear spectroscopy of condensed-phase systems. We showcase the DKT approach in the simulation of the fifth-order two-dimensional (2D) Raman spectroscopy of Lennard-Jones liquids as a prototypical example, which involves nontrivial nonlinear spectroscopic observables of systems described by anharmonic potentials. Our results show that the DKT can faithfully reproduce the 2D Raman response of liquid xenon at high temperatures, where the system behaves classically. In contrast, liquid neon at low temperatures exhibits moderate but discernible nuclear quantum effects in the 2D Raman response compared to the responses obtained with classical molecular dynamics approaches. Thus, the DKT formalism in combination with RPMD simulations enables simulations of multidimensional optical spectroscopy of condensed-phase systems that partially account for nuclear quantum effects.

Insulin Dissociates by Diverse Mechanisms of Coupled Unfolding and Unbinding

The protein hormone insulin exists in various oligomeric forms, and a key step in binding its cellular receptor is dissociation of the dimer. This dissociation process and its corresponding association process have come to serve as paradigms of coupled (un)folding and (un)binding more generally. Despite its fundamental and practical importance, the mechanism of insulin dimer dissociation remains poorly understood. Here, we use molecular dynamics simulations, leveraging recent developments in umbrella sampling, to characterize the energetic and structural features of dissociation in unprecedented detail. We find that the dissociation is inherently multipathway with limiting behaviors corresponding to conformational selection and induced fit, the two prototypical mechanisms of coupled folding and binding. Along one limiting path, the dissociation leads to detachment of the C-terminal segment of the insulin B chain from the protein core, a feature believed to be essential for receptor binding. We simulate IR spectroscopy experiments to aid in interpreting current experiments and identify sites where isotopic labeling can be most effective for distinguishing the contributions of the limiting mechanisms.

Ultrafast Spectroscopy of Lipid-Water Interfaces: Transmembrane Crowding Drives H-Bond Dynamics

Biology takes place in crowded, heterogeneous environments, and it is therefore essential to account for crowding effects in our understanding of biophysical processes at the molecular level. Comparable to the cytosol, proteins occupy approximately 30% of the plasma membrane surface; thus, crowding should have an effect on the local structure and dynamics at the lipid-water interface. Using a combination of ultrafast two-dimensional infrared spectroscopy and molecular dynamics simulations, we quantify the effects of membrane peptide concentration on the picosecond interfacial H-bond dynamics. The measurements reveal a nonmonotonic dependence of water orientation and dynamics as a function of transmembrane peptide:lipid ratio. We observe three dynamical regimes: a "pure lipid-like" regime at low peptide concentrations, a bulk-like region at intermediate peptide concentrations where dynamics are faster by ∼20% compared to those of the pure lipid bilayer, and a crowded regime where high peptide concentrations slow dynamics by ∼50%.