Infrared Probes for Studying the Structure and Dynamics of Biomolecules

Scattering Elimination in 2D IR Immune from Detector Artifacts

Highly scattering samples, such as polymer droplets or solid-state powders, are difficult to study via coherent two-dimensional infrared (2D IR) spectroscopy. Previously, researchers have employed (quasi-) phase cycling, local-oscillator chopping, and polarization control to reduce scattering, but the latter method poses a limit on polarization-dependent measurements. Here, we present a method for Scattering Elimination Immune from Detector Artifacts (SEIFDA) in pump-probe 2D IR experiments. Our method extends the negative probe delay method of removing scattering from pump-probe spectroscopy to 2D experiments. SEIFDA works well for all polarizations when combined with the optimized noise reduction scheme to remove additive and multiplicative noise. We demonstrate that our method can be employed with any polarization scheme and reliably lowers the scattering at parallel polarization to comparable levels to the conventional 8-frame phase cycling with probe chopping (8FPCPC) at perpendicular polarization. Our system can acquire artifact free spectra in parallel polarization when the signal intensity is as little as 5% of the intensity of the interference between the pump pulses scattered into the detector. It reduces the time required to characterize the scattering term by at least 50% over 8FPCPC. Through detailed analysis of detector nonlinearity, we show that the performance of 8FPCPC can be improved by incorporating nonlinear correction factors, but it is still worse than that of SEIFDA. Application of SEIFDA to study the encapsulation of Nile red in polymer droplets demonstrates that this method will be very useful for probing highly scattering systems.

SCN as a local probe of protein structural dynamics

The dynamics of lysozyme is probed by attaching -SCN to all alanine residues. The one-dimensional infrared spectra exhibit frequency shifts in the position of the maximum absorption of 4 cm-1, which is consistent with experiments in different solvents and indicates moderately strong interactions of the vibrational probe with its environment. Isotopic substitution 12C → 13C leads to a redshift by -47 cm-1, which agrees quantitatively with experiments for CN-substituted copper complexes in solution. The low-frequency, far-infrared part of the protein spectra contains label-specific information in the difference spectra when compared with the wild type protein. Depending on the position of the labels, local structural changes are observed. For example, introducing the -SCN label at Ala129 leads to breaking of the α-helical structure with concomitant change in the far-infrared spectrum. Finally, changes in the local hydration of SCN-labeled alanine residues as a function of time can be related to the reorientation of the label. It is concluded that -SCN is potentially useful for probing protein dynamics, both in the high-frequency part (CN-stretch) and in the far-infrared part of the spectrum.

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

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

Fermi Resonance of the N-D Stretching Mode Probing the Local Hydrogen-Bonding Environment in Proteins

Local H-bonding interactions are crucial for proteins to undergo various structural transitions and form different secondary structures. However, identifying slight distinctions in the local H-bonding of proteins is rather challenging. Here, we demonstrate that the Fermi resonance of the N-D stretching mode can provide an effective probe for the localized H-bonding environment of proteins both at the surface/interface and in the bulk. Using sum frequency generation vibrational spectroscopy and infrared spectroscopy, we established a correlation between the Fermi resonance of the N-D mode and protein secondary structures. The H-bond of N-D···C═O splits the N-D modes into two peaks (∼2410 and ∼2470 cm-1). The relative strength ratio (R) between the ∼2410 cm-1 peak and the ∼2470 cm-1 peak is very sensitive to H-bond strength and protein secondary structure. R is less than 1 for α-helical peptides, while R is greater than 1 for β-sheet peptides. For R < 2.5, both α-helical/loop structures and β-sheet structures exhibit almost identical Fermi coupling strengths (W = 28 cm-1). For R > 2.5, W decreases from 28 to 14 cm-1 and depends on the aggregation degree of the β-sheet oligomers or fibrils. The initial local H-bonding status impacts the misfolding dynamics of proteins at the lipid bilayer interface.

Unnatural Amino Acids for Biological Spectroscopy and Microscopy

Due to advances in methods for site-specific incorporation of unnatural amino acids (UAAs) into proteins, a large number of UAAs with tailored chemical and/or physical properties have been developed and used in a wide array of biological applications. In particular, UAAs with specific spectroscopic characteristics can be used as external reporters to produce additional signals, hence increasing the information content obtainable in protein spectroscopic and/or imaging measurements. In this Review, we summarize the progress in the past two decades in the development of such UAAs and their applications in biological spectroscopy and microscopy, with a focus on UAAs that can be used as site-specific vibrational, fluorescence, electron paramagnetic resonance (EPR), or nuclear magnetic resonance (NMR) probes. Wherever applicable, we also discuss future directions.

Early Aggregation Mechanism of SOD128-38 Based on Force Field Parameter of 5-Cyano-Tryptophan

The aggregation of superoxide dismutase 1 (SOD1) results in amyloid deposition and is involved in familial amyotrophic lateral sclerosis, a fatal motor neuron disease. There have been extensive studies of its aggregation mechanism. Noncanonical amino acid 5-cyano-tryptophan (5-CN-Trp), which has been incorporated into the amyloid segments of SOD1 as infrared probes to increase the structural sensitivity of IR spectroscopy, is found to accelerate the overall aggregation rate and potentially modulate the aggregation process. Despite these observations, the underlying mechanism remains elusive. Here, we optimized the force field parameters of 5-CN-Trp and then used molecular dynamics simulation along with the Markov state model on the SOD128-38 dimer to explore the kinetics of key intermediates in the presence and absence of 5-CN-Trp. Our findings indicate a significantly increased probability of protein aggregate formation in 5CN-Trp-modified ensembles compared to wildtype. Dimeric β-sheets of different natures were observed exclusively in the 5CN-Trp-modified peptides, contrasting with wildtype simulations. Free-energy calculations and detailed analyses of the dimer structure revealed augmented interstrand interactions attributed to 5-CN-Trp, which contributed more to peptide affinity than any other residues. These results explored the key events critical for the early nucleation of amyloid-prone proteins and also shed light on the practice of using noncanonical derivatives to study the aggregation mechanism.

Perturbation of Fermi Resonance on Hydrogen-Bonded > C═O: 2D IR Studies of Small Ester Probes

We utilized linear and 2D infrared spectroscopy to analyze the carbonyl stretching modes of small esters in different solvents. Particularly noteworthy were the distinct carbonyl spectral line shapes in aqueous solutions, prompting our investigation of the underlying factors responsible for these differences. Through our experimental and theoretical calculations, we identified the presence of the hydrogen-bond-induced Fermi resonance as the primary contributor to the varied line shapes of small esters in aqueous solutions. Furthermore, our findings revealed that the skeletal deformation mode plays a crucial role in the Fermi resonance for all small esters. Specifically, the first overtone band of the skeletal deformation mode intensifies when hydrogen bonds form with the carbonyl group of esters, whereas such coupling is rare in aprotic organic solvents. These spectral insights carry significant implications for the utilization of esters as infrared probes in both biological and chemical systems.

Nitrile Substituents at the Conjugated Dipyridophenazine Moiety as Infrared Redox Markers in Electrochemically Reduced Heteroleptic Ru(II) Polypyridyl Complexes

Ruthenium(II) complexes [Ru(tap)2(NN)]2+ (tap = 1,4,5,8-tetraazaphenanthrene, NN = 11-cyano-dipyrido[3,2-a:2',3'-c]phenazine (11-CN-dppz) and 11,12-dicyano-dipyrido[3,2-a:2',3'-c]phenazine (11,12-CN-dppz)) feature the C≡N groups as infrared (IR)-active redox markers. They were studied by cyclic voltammetry, UV-vis, and IR spectroelectrochemistry (SEC), and density functional theory calculations to assign the four 1e- reduction waves R1-R4 observed in dichloromethane. Generally, the NN ligands are reduced first (R1). For [Ru(tap)2(11,12-CN-dppz)]2+, R1 is sufficiently separated from R2 and delocalized over both tap ligands. Accordingly, IR SEC conducted at R1 shows a large red shift of the νs,as(C≡N) modes by -18/-28 cm-1, accompanied by a 4-fold enhancement of the νs(C≡N) intensity, comparably with reference data for free 11,12-CN-dppz. The first tap-based reduction of spin-doublet [Ru(tap)2(11,12-CN-dppz)]+ to spin-triplet [Ru(tap)2(11,12-CN-dppz)] at R2 decreased ν(C≡N) by merely -2 cm-1, while the intensity enhancement reached an overall factor of 8. Comparably, a red shift of ν(C≡N) by -27 cm-1 resulted from the 1e- reduction of [Ru(tap)2(11-CN-dppz)]2+ at R1 (poorly resolved from R2), and the intensity enhancement was roughly 3-fold. Concomitant 1e- reductions of the tap ligands (R2 and R3) caused only minor ν(C≡N) shifts of -3 cm-1 and increased the absorbance by overall factors of 6.5 and 8, respectively.

Site specifically probing the unfolding process of human telomere i-motif DNA using vibrationally enhanced alkynyl stretch

The microscopic unfolding process of a cytosine-rich DNA forming i-motif by hemi-protonated base pairs is related to gene regulation. However, the detailed thermal unfolding mechanism and the protonation/deprotonation status of site-specific cytosine in DNA in a physiological environment are still obscure. To address this issue, a vibration-enhanced CC probe tagged on 5'E terminal cytosine of human telomere i-motif DNA was examined using linear and nonlinear infrared (IR) spectroscopies and quantum-chemistry calculations. The CC probe extended into the major groove of the i-motif was found using nonlinear IR results only to introduce a minor steric effect on both steady-state structure and local structure dynamics; however, its IR absorption profile effectively reports the cleavage of the hemi-protonated base pair of C1-C13 upon the unfolding with C1 remaining protonated. The temperature mid-point (Tm) of the local transition reported using the CC tag was slightly lower than the Tm of global transition, and the enthalpy of the former exceeds 60% of the global transition. It is shown that the base-pair unraveling is noncooperative, with outer base pairs breaking first and being likely the rate limiting step. Our results offered an in-depth understanding of the macroscopic unfolding characteristics of the i-motif DNA and provided a nonlinear IR approach to monitoring the local structural transition and dynamics of DNA and its complexes.

Triple-Bond Vibrations: Emerging Applications in Energy and Biological Sciences

Triple bonds, such as that formed between two carbon atoms (i.e., C≡C) or that formed between one carbon atom and one nitrogen atom (i.e., C≡N), afford unique chemical bonding and hence vibrational characteristics. As such, they are not only frequently used to construct molecules with tailored chemical and/or physical properties but also employed as vibrational probes to provide site-specific chemical and/or physical information at the molecular level. Herein, we offer our perspective on the emerging applications of various triple-bond vibrations in energy and biological sciences with a focus on C≡C and C≡N triple bonds.

Anisotropic dynamics of an interfacial enzyme active site observed using tethered substrate analogs and ultrafast 2D IR spectroscopy

Enzymes accelerate the rates of biomolecular reactions by many orders of magnitude compared to bulk solution, and it is widely understood that this catalytic effect arises from a combination of polar pre-organization and electrostatic transition state stabilization. A number of recent reports have also implicated ultrafast (femtosecond-picosecond) timescale motions in enzymatic activity. However, complications arising from spatially-distributed disorder, the occurrence of multiple substrate binding modes, and the influence of hydration dynamics on solvent-exposed active sites still confound many experimental studies. Here we use ultrafast two-dimensional infrared (2D IR) spectroscopy and covalently-tethered substrate analogs to examine dynamical properties of the promiscuous Pyrococcus horikoshii ene-reductase (PhENR) active site in two binding configurations mimicking proposed "inactive" and "reactive" Michaelis complexes. Spectral diffusion measurements of aryl-nitrile substrate analogs reveal an end-to-end tradeoff between fast (sub-ps) and slow (>5 ps) motions. Fermi resonant aryl-azide analogs that sense interactions of coupled oscillators are described. Lineshape and quantum beat analyses of these probes reveal characteristics that correlate with aryl-nitrile frequency fluctuation correlation functions parameters, demonstrating that this anisotropy is an intrinsic property of the water-exposed active site, where countervailing gradients of fast dynamics and disorder in the reactant ground state are maintained near the hydration interface. Our results suggest several plausible factors leading to state-selective rate enhancement and promiscuity in PhENR. This study also highlights a strategy to detect perturbations to vibrational modes outside the transparent window of the mid-IR spectrum, which may be extended to other macromolecular systems.

Effects of the Intramolecular Group and Solvent on Vibrational Coupling Modes and Strengths of Fermi Resonances in Aryl Azides: A DFT Study of 4-Azidotoluene and 4-Azido-N-phenylmaleimide

The high transition dipole strength of the azide asymmetric stretch makes aryl azides good candidates as vibrational probes (VPs). However, aryl azides have complex absorption profiles due to Fermi resonances (FRs). Understanding the origin and the vibrational modes involved in FRs of aryl azides is critically important toward developing them as VPs for studies of protein structures and structural changes in response to their surroundings. As such, we studied vibrational couplings in 4-azidotoluene and 4-azido-N-phenylmaleimide in two solvents, N,N-dimethylacetamide and tetrahydrofuran, to explore the origin and the effects of intramolecular group and solvent on the FRs of aryl azides using density functional theory (DFT) calculations with the B3LYP functional and seven basis sets, 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p), and 6-311++G(df,pd). Two combination bands consisting of the azide symmetric stretch and another mode form strong FRs with the azide asymmetric stretch for both molecules. The FR profile was altered by replacing the methyl group with maleimide. Solvents change the relative peak position and intensity more significantly for 4-azido-N-phenylmaleimide, which makes it a more sensitive VP. Furthermore, the DFT results indicate that a comparison among the results from different basis sets can be used as a means to predict more reliable vibrational spectra.

C≡N Stretching Frequency as a Convenient Reporter of Charge Separation in Molecular Systems

Previously, several studies have shown that, for a set of structurally related nitrile compounds, there could be a linear relationship between the total charge on the nitrile group (qCN) and its stretching frequency (νCN). However, it is unclear whether the corresponding frequency and charge properties of structurally different nitrile compounds can be described by a single linear νCN-qCN relationship. Herein, we compute the qCN magnitudes of a large number of nitrile-containing molecules whose νCN values cover a spectral range of ca. 200 cm-1 and are measured under different experimental conditions. Our results reveal that νCN indeed exhibits a linear dependence on qCN, with a slope of 637 ± 30 cm-1/charge. Because the nitrile moiety is a commonly used building block in electronic donor-acceptor (D-A) molecular systems, we believe that this linear relationship will find utility in a wide range of applications where such D-A constructs are used, such as in organic photovoltaic assemblies. In addition, we apply this linear relationship to characterize the degree of charge transfer upon photoexcitation of two indole derivatives, 5-cyanoindole and 6-cyanoindole, and are able to show that in both cases, the fluorescence emission arises from a charge-transfer or La state.

Combining fragmentation method and high-performance computing: Geometry optimization and vibrational spectra of proteins

Exploring the structures and spectral features of proteins with advanced quantum chemical methods is an uphill task. In this work, a fragment-based molecular tailoring approach (MTA) is appraised for the CAM-B3LYP/aug-cc-pVDZ-level geometry optimization and vibrational infrared (IR) spectra calculation of ten real proteins containing up to 407 atoms and 6617 basis functions. The use of MTA and the inherently parallel nature of the fragment calculations enables a rapid and accurate calculation of the IR spectrum. The applicability of MTA to optimize the protein geometry and evaluate its IR spectrum employing a polarizable continuum model with water as a solvent is also showcased. The typical errors in the total energy and IR frequencies computed by MTA vis-à-vis their full calculation (FC) counterparts for the studied protein are 5-10 millihartrees and 5 cm-1, respectively. Moreover, due to the independent execution of the fragments, large-scale parallelization can also be achieved. With increasing size and level of theory, MTA shows an appreciable advantage in computer time as well as memory and disk space requirement over the corresponding FCs. The present study suggests that the geometry optimization and IR computations on the biomolecules containing ∼1000 atoms and/or ∼15 000 basis functions using MTA and HPC facility can be clearly envisioned in the near future.

Photoluminescent ellipsometric circular dichroism

A novel spectroscopic technique, photoluminescent ellipsometric circular dichroism (PECD), which distinguishes all radiative electronic transitions related to molecular chiral centers. Additionally, it is proposed as complementary to the ellipsometric Raman spectroscopy (ERS) technique, thus establishing a relationship between vibrational modes and electronic transitions, associated with molecular chiral centers. In this way, PECD turns into a powerful technique for chiral material characterization. The PECD technique was performed on a chiral oligomer (1R,2R)-diiminocyclohexane, and its derivative polymer. A complete photophysical characterization in solution was performed to corroborate the new PECD technique.

Dissecting amide-I vibrations in histidine dipeptide

The amide-I vibrational characteristics and conformational preferences of the model compound - histidine dipeptide (Ac-His-NHCH₃, HISD) in gas phase and solution have been revealed with the help of ab initio calculations and wavefunction analyses. The Gibbs free energy surfaces (FESs) of solvated HISD were smoothed by solvent effect to exhibit different structural populations concerning various external environments. It was shown that the most stable conformations of HISD in CHCl₃ and gas phase are C7eq, while those in DMSO and water are β and PPII, respectively. Compared with ALAD, the number of accessible conformational states on these FESs was predicted to be reduced due to the steric effect of imidazole group. The two amide-I normal modes of HISD were found to have intrinsically secondary structural dependencies, and be sensitive to surrounding environments. The average amide-I_a frequencies of HISD isomers in these environments were predicted to be almost the same as those of ALAD, while the amide-I_b mean frequencies were estimated to be lower than ALAD due to the intramolecular interactions between the imidazole group and amino-terminal amide unit. The good linear correlations between amide-I frequencies and the atomic electrostatic potentials (ESPs) of amide groups were also found to interpret the solvent-induced amide-I frequency shifts of HISD at the electronic structure level. These results allow us to gain a deep understanding of amide-I vibrations of HISD, and would be helpful for the site-specific conformational monitoring and spectral interpretation of solvated polypeptides.

Distinguishing between the Electrostatic Effects and Explicit Ion Interactions in a Stark Probe

Vibrational Stark probes are incisive tools for measuring local electric fields in a wide range of chemical environments. The interpretation of the frequency shift often gets complicated due to the specific interactions of the probe, such as hydrogen bonding and Lewis bonding. Therefore, it is important to distinguish between the pure electrostatic response and the response due to such specific interactions. Here we report a molecular system that is sensitive to both the Stark effect from a single ion and the explicit Lewis bonding of ions with the probe. The molecule consists of a crown ether with an appended benzonitrile. The crown captures cations of various charges, and the electric field from the ions is sensed by the benzonitrile probe. Additionally, the lone pair of the benzonitrile can engage in Lewis interactions with some of the ions by donating partial charge density to the ions. Our system exhibits both of these effects and therefore is a suitable test bed for distinguishing between the pure electrostatic and the Lewis interactions. Our computational results show that the electrostatic influence of the ion is operative at large distances, while the Lewis interaction becomes important only within distances that permit orbital overlap. Our results may be useful for using the nitrile probe for measuring electrostatic and coordination effects in complex ionic environments such as the electrode-electrolyte interfaces.

Inconsistent hydrogen bond-mediated vibrational coupling of amide I

Using infrared spectroscopy and density functional theory (DFT) calculations, we scrutinized an amide (dimethylformamide) as a "model" compound to interpret the interactions of amide 1 with different phenol derivatives (para-chlorophenol (PCP) and para-cresol (CP)) as "model guest molecules". We established the involvement of amide I in vibrational coupling with symmetric and asymmetric C[double bond, length as m-dash]C modes of different phenolic derivatives and how their coupling was dependent upon different guest aromatic phenolic compounds. Interestingly, substitution of phenol perturbed the pattern of vibrational coupling with amide I. The symmetric and asymmetric C[double bond, length as m-dash]C modes of PC were coupled significantly with amide 1. For PCP, the symmetric C[double bond, length as m-dash]C mode coupled significantly, but the asymmetric mode coupled negligibly, with amide I. Here, we reveal the nature of vibrational coupling based on the structure of a guest molecule hydrogen-bonded with amide I. Our conclusions could be valuable for depiction of the unusual dynamics of coupled amide-I modes as well as the dependency of vibrational coupling on altered factors.

Ion-specific binding of cations to the carboxylate and of anions to the amide of alanylalanine

Studies of ion-specific effects on oligopeptides have aided our understanding of Hofmeister effects on proteins, yet the use of different model peptides and different experimental sensitivities have led to conflicting conclusions. To resolve these controversies, we study a small model peptide, L-Alanyl-L-alanine (2Ala), carrying all fundamental chemical protein motifs: C-terminus, amide bond, and N-terminus. We elucidate the effect of GdmCl, LiCl, KCl, KI, and KSCN by combining dielectric relaxation, nuclear magnetic resonance (¹H-NMR), and (two-dimensional) infrared spectroscopy. Our dielectric results show that all ions reduce the rotational mobility of 2Ala, yet the magnitude of the reduction is larger for denaturing cations than for anions. The NMR chemical shifts of the amide group are particularly sensitive to denaturing anions, indicative of anion-amide interactions. Infrared experiments reveal that LiCl alters the spectral homogeneity and dynamics of the carboxylate, but not the amide group. Interaction of LiCl with the negatively charged pole of 2Ala, the COO^- group, can explain the marked cationic effect on dipolar rotation, while interaction of anions between the poles, at the amide, only weakly perturbs dipolar dynamics. As such, our results provide a unifying view on ions' preferential interaction sites at 2Ala and help rationalize Hofmeister effects on proteins.

Probing the local structure and dynamics of nucleotides using vibrationally enhanced alkynyl stretching

Monitoring the site-specific local structure and dynamics of polynucleotides and DNA is important for understanding their biological functions. However, structurally characterizing these biomolecules with high time resolution has been known to be experimentally challenging. In this work, several 5-silylethynyl-2'-deoxynucleosides and 5-substituted phenylethynyl-2'-deoxynucleosides on the basis of deoxycytidine (dC) and deoxythymidine (dT) were synthesized, in which the alkynyl group shows intensified CC stretching vibration with infrared transition dipole moment magnitude close to that of typical CO stretching, and exhibits structural sensitivities in both vibrational frequency and spectral width. In particular, 5-trimethylsilylethynyl-2'-dC (^TMSEdC, molecule 1a) was examined in detail using femtosecond nonlinear IR spectroscopy. The solvent dependent CC stretching frequency of 1a can be reasonably interpreted mainly as the hydrogen-bonding effect between the solvent and cytosine base ring structure. Transient 2D IR and pump-probe IR measurements of 1a carried out comparatively in two aprotic solvents (DMSO and THF) and one protic solvent (MeOH) further reveal solvent dependent ultrafast vibrational properties, including diagonal anharmonicity, spectral diffusion, vibrational relaxation and anisotropy dynamics. These observed sensitivities are rooted in an extended π-conjugation of the base ring structure in which the CC group is actively involved. Our results show that the intensified CC stretching vibration can potentially provide a site-specific IR probe for monitoring the equilibrium and ultrafast structural dynamics of polynucleotides.

Molecular Simulation-Guided Spectroscopy of Imidazolium-Based Ionic Liquids and Effects of Methylation on Ion-Cage and -Pair Dynamics

Classical molecular dynamics simulations were performed to assess an atomistic interpretation of the ion-probe structural interactions in two typical ionic liquids (ILs), 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide [BMIm][NTf₂] and 1-butyl-2,3-dimethylimidazolium bis(trifluoromethylsulfonyl)imide [BDimIm][NTf₂] through computational ultrafast spectroscopy. The nitrile stretching vibrations of the thiocyanate anion, [SCN]^-, serve as the local mode of the ultrafast system dynamics within the imidazolium-based ionic liquid environment. The wavelet transform of classical trajectories determines the time-varying fluctuating frequencies and the stretch spectral signatures of SCN^- in the normalized distribution. However, computational modeling of the two-dimensional (2D) spectra from the wavelet-derived vibrational frequencies yields time evolution of the local molecular structure along with the varied time-dependent dynamics of the spectral diffusion process. We calculated the frequency-frequency correlation functions (FFCFs), time correlations associated with the ion-pair and -cage dynamics, and mean square displacements as a function of time, depicting diffusive dynamics. The calculated results based on the pair correlation functions and the distribution of atomic density suggest that the hydrogen and methylated carbon at the two-position of the imidazolium ring of [BMIm] and [BDimIm] cations, respectively, strongly interact with the probe through the N of the thiocyanate anion rather than the S atom. The center-of-mass center-of-mass (COM-COM) cation-probe radial distribution functions (RDFs) in conjunction with the site-specific structural analysis further reveal well-structured interactions of the thiocyanate ion and [BMIm]⁺ cation rather than the [BDimIm] cation. In contrast, the anion-probe COM-COM RDFs depict weak interactive associations within the vibrational probe [SCN]^- and [NTf₂]^- ions. Methylation at the two-position of the imidazolium ring predicts slower structural reorganization and breaking and reformation dynamics of the ion pairs and cages within the ionic liquid framework.

Development of Efficient Computational Analysis of Difference Infrared and Raman Spectroscopies

Computational analysis of difference spectra between two analogous systems is a challenging issue, as reliable estimation of a tiny difference spectrum requires an extraordinary precision of the two original spectra. We have developed an alternative new method to calculate the difference spectra in background-free conditions, which greatly improved the efficiency of computation. In this paper we report further improvement by using efficient parallel implementation and the time correlation formula based on time derivative quantities. As a consequence, the present work achieved further remarkable acceleration in the calculations of difference infrared and Raman spectra in the order of magnitude, and thereby allowed us with analyzing these difference spectra at a practical cost of computation.

Biomacromolecule-Assisted Screening for Reaction Discovery and Catalyst Optimization

Reaction discovery and catalyst screening lie at the heart of synthetic organic chemistry. While there are efforts at de novo catalyst design using computation/artificial intelligence, at its core, synthetic chemistry is an experimental science. This review overviews biomacromolecule-assisted screening methods and the follow-on elaboration of chemistry so discovered. All three types of biomacromolecules discussed─enzymes, antibodies, and nucleic acids─have been used as "sensors" to provide a readout on product chirality exploiting their native chirality. Enzymatic sensing methods yield both UV-spectrophotometric and visible, colorimetric readouts. Antibody sensors provide direct fluorescent readout upon analyte binding in some cases or provide for cat-ELISA (Enzyme-Linked ImmunoSorbent Assay)-type readouts. DNA biomacromolecule-assisted screening allows for templation to facilitate reaction discovery, driving bimolecular reactions into a pseudo-unimolecular format. In addition, the ability to use DNA-encoded libraries permits the barcoding of reactants. All three types of biomacromolecule-based screens afford high sensitivity and selectivity. Among the chemical transformations discovered by enzymatic screening methods are the first Ni(0)-mediated asymmetric allylic amination and a new thiocyanopalladation/carbocyclization transformation in which both C-SCN and C-C bonds are fashioned sequentially. Cat-ELISA screening has identified new classes of sydnone-alkyne cycloadditions, and DNA-encoded screening has been exploited to uncover interesting oxidative Pd-mediated amido-alkyne/alkene coupling reactions.

Ligand Dynamics Time Scales Identify the Surface-Ligand Interactions in Thiocyanate-Capped Cadmium Sulfide Nanocrystals

The nanocrystal surface, which acts as an interface between the semiconductor lattice and the capping ligands, plays a significant role in the attractive photophysical properties of semiconductor nanocrystals for use in a wide range of applications. Replacing the long-chain organic ligands with short inorganic variants improves the conductivity and carrier mobility of nanocrystal-based devices. However, our current understanding of the interactions between the inorganic ligands and the nanocrystals is obscure due to the lack of experiments to directly probe the inorganic ligands. Herein, using two-dimensional infrared spectroscopy, we show that the variations in the inorganic ligand dynamics within the heterogeneous nanocrystal ensemble can identify the diversities in the inorganic ligand-nanocrystal interactions. The ligand dynamics time scale in SCN^- capped CdS nanocrystals identifies three distinct ligand populations and provides molecular insight into the nanocrystal surface. Our results demonstrate that the SCN^- ligands engage in a dynamic equilibrium and stabilize the nanocrystals by neutralizing the surface charges through both direct binding and electrostatic interaction.

Origin of thiocyanate spectral shifts in water and organic solvents

Vibrational spectroscopy is a useful technique for probing chemical environments. The development of models that can reproduce the spectra of nitriles and azides is valuable because these probes are uniquely suited for investigating complex systems. Empirical vibrational spectroscopic maps are commonly employed to obtain the instantaneous vibrational frequencies during molecular dynamics simulations but often fail to adequately describe the behavior of these probes, especially in its transferability to a diverse range of environments. In this paper, we demonstrate several reasons for the difficulty in constructing a general-purpose vibrational map for methyl thiocyanate (MeSCN), a model for cyanylated biological probes. In particular, we found that electrostatics alone are not a sufficient metric to categorize the environments of different solvents, and the dominant features in intermolecular interactions in the energy landscape vary from solvent to solvent. Consequently, common vibrational mapping schemes do not cover all essential interaction terms adequately, especially in the treatment of van der Waals interactions. Quantum vibrational perturbation (QVP) theory, along with a combined quantum mechanical and molecular mechanical potential for solute-solvent interactions, is an alternative and efficient modeling technique, which is compared in this paper, to yield spectroscopic results in good agreement with experimental FTIR. QVP has been used to analyze the computational data, revealing the shortcomings of the vibrational maps for MeSCN in different solvents. The results indicate that insights from QVP analysis can be used to enhance the transferability of vibrational maps in future studies.

Population Dynamics of Stretching Excitations of p-Azido-phenylalanine Incorporated in Calmodulin-Peptide Complexes

We genetically incorporated the unnatural amino acid p-azido-phenylalanine (AzF) into the ubiquitous Ca²⁺ sensor protein calmodulin (CaM) in complex with different peptides to explore the response of the azido stretching line shape to varying binding motifs with femtosecond infrared spectroscopy. The dynamic response of the azido stretching mode varies in different CaM-peptide complexes. We model these dynamics as coherent excitations of Fermi resonances and extract a lifetime of the azido stretching vibration of about 1 ps. The resulting model parameters are commensurate with the linear infrared absorption lineshapes which suggests that the conformation-sensitive vibrational lineshape could be composed of Fermi resonances that differ between the protein-peptide complexes.

Local Charge Distributions, Electric Dipole Moments, and Local Electric Fields Influence Reactivity Patterns and Guide Regioselectivities in α-Ketoglutarate-Dependent Non-heme Iron Dioxygenases

Non-heme iron dioxygenases catalyze vital processes for human health related to the biosynthesis of essential products and the biodegradation of toxic metabolites. Often the natural product biosyntheses by these non-heme iron dioxygenases is highly regio- and chemoselective, which are commonly assigned to tight substrate-binding and positioning. However, recent high-level computational modeling has shown that substrate-binding and positioning is only part of the story and long-range electrostatic interactions can play a major additional role.In this Account, we review and summarize computational viewpoints on the high regio- and chemoselectivity of α-ketoglutarate-dependent non-heme iron dioxygenases and how external perturbations affect the catalysis. In particular, studies from our groups have shown that often a regioselectivity in enzymes can be accomplished by stabilization of the rate-determining transition state for the reaction through external charges, electric dipole moments, or local electric field effects. Furthermore, bond dissociation energies in molecules are shown to be influenced by an electric field effect, and through targeting a specific bond in an electric field, this can lead to an unusually specific reaction. For instance, in the carbon-induced starvation protein, we studied two substrate-bound conformations and showed that regardless of what C-H bond of the substrate is closest to the iron(IV)-oxo oxidant, the lowest hydrogen atom abstraction barrier is always for the pro-S C²-H abstraction due to an induced dipole moment of the protein that weakens this bond. In another example of the hygromycin biosynthesis enzyme, an oxidative ring-closure reaction in the substrate forms an ortho-δ-ester ring. Calculations on this enzyme show that the selectivity is guided by a protonated lysine residue in the active site that, through its positive charge, triggers a low energy hydrogen atom abstraction barrier. A final set of examples in this Account discuss the viomycin biosynthesis enzyme and the 2-(trimethylammonio)ethylphosphonate dioxygenase (TmpA) enzyme. Both of these enzymes are shown to possess a significant local dipole moment and local electric field effect due to charged residues surrounding the substrate and oxidant binding pockets. The protein dipole moment and local electric field strength changes the C-H bond strengths of the substrate as compared to the gas-phase triggers the regioselectivity of substrate activation. In particular, we show that in the gas phase and in a protein environment C-H bond strengths are different due to local electric dipole moments and electric field strengths. These examples show that enzymes have an intricately designed structure that enables a chemical reaction under ambient conditions through the positioning of positively and negatively charged residues that influence and enhance reaction mechanisms. These computational insights create huge possibilities in bioengineering to apply local electric field and dipole moments in proteins to achieve an unusual selectivity and specificity and trigger a fit-for-purpose biocatalyst for unique biotransformations.

Quantum vibration perturbation approach with polyatomic probe in simulating infrared spectra

The quantitative prediction of vibrational spectra of chromophore molecules in solution is challenging and numerous methods have been developed. In this work, we present a quantum vibration perturbation (QVP) approach, which is a procedure that combines molecular quantum vibration and molecular dynamics with perturbation theory. In this framework, an initial Newtonian molecular dynamics simulation is performed, followed by a substitution process to embed molecular quantum vibrational wave functions into the trajectory. The instantaneous vibrational frequency shift at each time step is calculated using the Rayleigh-Schrödinger perturbation theory, where the perturbation operator is the difference in the vibrational potential between the reference chromophore and the perturbed chromophore in the environment. Semi-classical statistical mechanics is employed to obtain the spectral lineshape function. We validated our method using HCOOH·nH₂O (n = 1-2) clusters and HCOOH aqueous solution as examples. The QVP method can be employed for rapid prediction of the vibrational spectrum of a specific mode in solution.

Adsorbed Water Structure on Acrylate-Based Biocompatible Polymer Surface

The role of water in the excellent biocompatibility of the acrylate-based polymers widely used for antibiofouling coating material has been realized previously. Here, we report femtosecond mid-infrared pump-probe spectroscopy of the OD stretch band of HOD molecule adsorbed on highly biocompatible poly(2-methoxyethyl) acrylate [PMEA] and poorly biocompatible poly(2-phenoxyethyl) acrylate [PPEA], both of which reveal that there are two water species with significantly different vibrational lifetime. PMEA interacts more strongly with water than PPEA through the H-bonding interaction between carbonyl (C═O) and water. The vibrational lifetime of the OD stretch in PPEA is notably longer by factors of 3 and 7 than those in PMEA and bulk water, respectively. The IR-pump visible-probe photothermal imaging further unravels substantial spatial overlap between polymer CO group and water for hydrated PMEA and a significant difference in surface morphology than those in PPEA, which exhibits the underlying relationships among polymer-water interaction, surface morphology, and biocompatibility.

Ultrafast vibrational dynamics of the tyrosine ring mode and its application to enkephalin insertion into phospholipid membranes as probed by two-dimensional infrared spectroscopy

Enkephalins are small opioid peptides whose binding conformations are catalyzed by phospholipid membranes. Binding to opioid receptors is determined by the orientation of tyrosine and phenylalanine side chains. In this work, we investigate the effects of different charged phospholipid headgroups on the insertion of the tyrosine side chain into a lipid bilayer using a combination of 2D IR spectroscopy, anharmonic DFT calculations, and third order response function modeling. The insertion is probed by using the ∼1515 cm^-1 tyrosine ring breathing mode, which we found exhibits rich vibrational dynamics on the picosecond timescale. These dynamics include rapid intramolecular vibrational energy redistribution (IVR), where some of the energy ends up in a dark state that shows up as an anharmonically shifted combination band. The waiting-time dependent 2D IR spectra also show an unusual line shape distortion that affects the extraction of the frequency-frequency correlation function (FFCF), which is the dynamic observable of interest that reflects the tyrosine side chain's insertion into the lipid bilayer. We proposed three models to account for this distortion: a hot-state exchange model, a local environment dependent IVR model, and a coherence transfer model. A qualitative analysis of these models suggests that the local environment dependent IVR rate best explains the line shape distortion, while the coherence transfer model best reproduced the effects on the FFCF. Even with these complex dynamics, we found that the tyrosine ring mode's FFCF is qualitatively correlated with the degree of insertion expected from the different phospholipid headgroups.

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.

Local Electric Fields in Aqueous Electrolytes

Vibrational Stark shifts were explored in aqueous solutions of organic molecules with carbonyl- and nitrile-containing constituents. In many cases, the vibrational resonances from these moieties shifted toward lower frequency as salt was introduced into solution. This is in contrast to the blue-shift that would be expected based upon Onsager's reaction field theory. Salts containing well-hydrated cations like Mg²⁺ or Li⁺ led to the most pronounced Stark shift for the carbonyl group, while poorly hydrated cations like Cs⁺ had the greatest impact on nitriles. Moreover, salts containing I^- gave rise to larger Stark shifts than those containing Cl^-. Molecular dynamics simulations indicated that cations and anions both accumulate around the probe in an ion- and probe-dependent manner. An electric field was generated by the ion pair, which pointed from the cation to the anion through the vibrational chromophore. This resulted from solvent-shared binding of the ions to the probes, consistent with their positions in the Hofmeister series. The "anti-Onsager" Stark shifts occur in both vibrational spectroscopy and fluorescence measurements.

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.

Site-by-site tracking of signal transduction in an azidophenylalanine-labeled bacteriophytochrome with step-scan FTIR spectroscopy

Signal propagation in photosensory proteins is a complex and multidimensional event. Unraveling such mechanisms site-specifically in real time is an eligible but a challenging goal. Here, we elucidate the site-specific events in a red-light sensing phytochrome using the unnatural amino acid azidophenylalanine, vibrationally distinguishable from all other protein signals. In canonical phytochromes, signal transduction starts with isomerization of an excited bilin chromophore, initiating a multitude of processes in the photosensory unit of the protein, which eventually control the biochemical activity of the output domain, nanometers away from the chromophore. By implementing the label in prime protein locations and running two-color step-scan FTIR spectroscopy on the Deinococcus radiodurans bacteriophytochrome, we track the signal propagation at three specific sites in the photosensory unit. We show that a structurally switchable hairpin extension, a so-called tongue region, responds to the photoconversion already in microseconds and finalizes its structural changes concomitant with the chromophore, in milliseconds. In contrast, kinetics from the other two label positions indicate that the site-specific changes deviate from the chromophore actions, even though the labels locate in the chromophore vicinity. Several other sites for labeling resulted in impaired photoswitching, low structural stability, or no changes in the difference spectrum, which provides additional information on the inner dynamics of the photosensory unit. Our work enlightens the multidimensionality of the structural changes of proteins under action. The study also shows that the signaling mechanism of phytochromes is accessible in a time-resolved and site-specific approach by azido probes and demonstrates challenges in using these labels.

Advanced Materials for Energy-Water Systems: The Central Role of Water/Solid Interfaces in Adsorption, Reactivity, and Transport

The structure, chemistry, and charge of interfaces between materials and aqueous fluids play a central role in determining properties and performance of numerous water systems. Sensors, membranes, sorbents, and heterogeneous catalysts almost uniformly rely on specific interactions between their surfaces and components dissolved or suspended in the water-and often the water molecules themselves-to detect and mitigate contaminants. Deleterious processes in these systems such as fouling, scaling (inorganic deposits), and corrosion are also governed by interfacial phenomena. Despite the importance of these interfaces, much remains to be learned about their multiscale interactions. Developing a deeper understanding of the molecular- and mesoscale phenomena at water/solid interfaces will be essential to driving innovation to address grand challenges in supplying sufficient fit-for-purpose water in the future. In this Review, we examine the current state of knowledge surrounding adsorption, reactivity, and transport in several key classes of water/solid interfaces, drawing on a synergistic combination of theory, simulation, and experiments, and provide an outlook for prioritizing strategic research directions.

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.

Dynamics underlying hydroxylation selectivity of cytochrome P450cam

Structural heterogeneity and the dynamics of the complexes of enzymes with substrates can determine the selectivity of catalysis; however, fully characterizing how remains challenging as heterogeneity and dynamics can vary at the spatial level of an amino acid residue and involve rapid timescales. We demonstrate the nascent approach of site-specific two-dimensional infrared (IR) spectroscopy to investigate the archetypical cytochrome P450, P450cam, to better delineate the mechanism of the lower regioselectivity of hydroxylation of the substrate norcamphor in comparison to the native substrate camphor. Specific locations are targeted throughout the enzyme by selectively introducing cyano groups that have frequencies in a spectrally isolated region of the protein IR spectrum as local vibrational probes. Linear and two-dimensional IR spectroscopy were applied to measure the heterogeneity and dynamics at each probe and investigate how they differentiate camphor and norcamphor recognition. The IR data indicate that the norcamphor complex does not fully induce a large-scale conformational change to a closed state of the enzyme adopted in the camphor complex. Additionally, a probe directed at the bound substrate experiences rapidly interconverting states in the norcamphor complex that explain the hydroxylation product distribution. Altogether, the study reveals large- and small-scale structural heterogeneity and dynamics that could contribute to selectivity of a cytochrome P450 and illustrates the approach of site-selective IR spectroscopy to elucidate protein dynamics.

Optical and Infrared Spectroelectrochemical Studies of CN-Substituted Bipyridyl Complexes of Ruthenium(II)

Ruthenium(II) polypyridyl complexes [Ru(CN-Me-bpy)_x(bpy)_3-x]²⁺ (CN-Me-bpy = 4,4'-dicyano-5,5'-dimethyl-2,2'-bipyridine, bpy = 2,2'-bipyridine, and x = 1-3, abbreviated as 1²⁺, 2²⁺, and 3²⁺) undergo four (1²⁺) or five (2²⁺ and 3²⁺) successive one-electron reduction steps between -1.3 and -2.75 V versus ferrocenium/ferrocene (Fc⁺/Fc) in tetrahydrofuran. The CN-Me-bpy ligands are reduced first, with successive one-electron reductions in 2²⁺ and 3²⁺ being separated by 150-210 mV; reduction of the unsubstituted bpy ligand in 1²⁺ and 2²⁺ occurs only when all CN-Me-bpy ligands have been converted to their radical anions. Absorption spectra of the first three reduction products of each complex were measured across the UV, visible, near-IR (NIR), and mid-IR regions and interpreted with the help of density functional theory calculations. Reduction of the CN-Me-bpy ligand shifts the ν(C≡N) IR band by ca. -45 cm^-1, enhances its intensity ∼35 times, and splits the symmetrical and antisymmetrical modes. Semireduced complexes containing two and three CN-derivatized ligands 2⁺, 3⁺, and 3⁰ show distinct ν(C≡N) features due to the presence of both CN-Me-bpy and CN-Me-bpy^•-, confirming that each reduction is localized on a single ligand. NIR spectra of 1⁰, 1^-, and 2^- exhibit a prominent band attributable to the CN-Me-bpy^•- moiety between 6000 and 7500 cm^-1, whereas bpy^•--based absorption occurs between 4500 and 6000 cm^-1; complexes 2⁺, 3⁺, and 3⁰ also exhibit a band at ca. 3300 cm^-1 due to a CN-Me-bpy^•- → CN-Me-bpy interligand charge-transfer transition. In the UV-vis region, the decrease of π → π* intraligand bands of the neutral ligands and the emergence of the corresponding bands of the radical anions are most diagnostic. The first reduction product of 1²⁺ is spectroscopically similar to the lowest triplet metal-to-ligand charge-transfer excited state, which shows pronounced NIR absorption, and its ν(C≡N) IR band is shifted by -38 cm^-1 and 5-7-fold-enhanced relative to the ground state.

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.

Spectroscopic Probe Molecule Selection Using Quantum Theory, First-Principles Calculations, and Machine Learning

Probe molecule vibrational spectra have a long history of being used to characterize materials including metals, oxides, metal-organic frameworks, and even human proteins. Furthermore, recent advances in machine learning have enabled computationally generated spectra to aid in detailed characterization of complex surfaces with probe molecules. Despite widespread use of probe molecules, the science of probe molecule selection is underdeveloped. Here, we develop physical concepts, including orbital interaction energy and the energy overlap integral, to explain and predict the ability of probe molecules to discriminate structural descriptors. We resolve the crystal orbital overlap population (COOP) to specific molecular orbitals and quantify their bonding character, which directly influences vibrational frequencies. Using only a single adsorbate calculation from density function theory (DFT), we compute the interaction energy of individual adsorbate molecular orbitals with adsorption site atomic orbitals across many different sites. Combining the molecular orbital resolved COOP and changes in orbital interaction energy enables probe molecule selection for improved discrimination of various sites. We demonstrate these concepts by comparing the predicted effectiveness of carbon monoxide (CO), nitric oxide (NO), and ethylene (C₂H₄) to probe Pt adsorption sites. Finally, using a previously developed machine learning framework, we show that models trained on hundreds of thousands of C₂H₄ spectra, computed from DFT, which regress surface binding-type and generalized coordination number, outperform those trained using CO and NO spectra. A python package, pDOS_overlap, for implementing the electron density-based analysis on any combination of adsorbates and materials, is also made available.

Controls of nature: Secondary, tertiary, and quaternary structure of the enamel protein amelogenin in solution and on hydroxyapatite

Amelogenin, a protein critical to enamel formation, is presented as a model for understanding how the structure of biomineralization proteins orchestrate biomineral formation. Amelogenin is the predominant biomineralization protein in the early stages of enamel formation and contributes to the controlled formation of hydroxyapatite (HAP) enamel crystals. The resulting enamel mineral is one of the hardest tissues in the human body and one of the hardest biominerals in nature. Structural studies have been hindered by the lack of techniques to evaluate surface adsorbed proteins and by amelogenin's disposition to self-assemble. Recent advancements in solution and solid state nuclear magnetic resonance (NMR) spectroscopy, atomic force microscopy (AFM), and recombinant isotope labeling strategies are now enabling detailed structural studies. These recent studies, coupled with insights from techniques such as CD and IR spectroscopy and computational methodologies, are contributing to important advancements in our structural understanding of amelogenesis. In this review we focus on recent advances in solution and solid state NMR spectroscopy and in situ AFM that reveal new insights into the secondary, tertiary, and quaternary structure of amelogenin by itself and in contact with HAP. These studies have increased our understanding of the interface between amelogenin and HAP and how amelogenin controls enamel formation.

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.

Proving and Probing the Presence of the Elusive C-H⋅⋅⋅O Hydrogen Bond in Liquid Solutions at Room Temperature

Hydrogen bonds (H bonds) play a major role in defining the structure and properties of many substances, as well as phenomena and processes. Traditional H bonds are ubiquitous in nature, yet the demonstration of weak H bonds that occur between a highly polarized C-H group and an electron-rich oxygen atom, has proven elusive. Detailed here are linear and nonlinear IR spectroscopy experiments that reveal the presence of H bonds between the chloroform C-H group and an amide carbonyl oxygen atom in solution at room temperature. Evidence is provided for an amide solvation shell featuring two clearly distinguishable chloroform arrangements that undergo chemical exchange with a time scale of about 2 ps. Furthermore, the enthalpy of breaking the hydrogen bond is found to be 6-20 kJ mol^-1 . Ab-initio computations support the findings of two distinct solvation shells formed by three chloroform molecules, where one thermally undergoes hydrogen-bond making and breaking.

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.