Site-specific coupling of hydration water and protein flexibility studied in solution with ultrafast 2D-IR spectroscopy

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.

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.

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.

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.

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.

A Systematic Study on the Role of Hydrogen Bond Donors in Dictating the Dynamics of Choline-Based Deep Eutectic Solvents

Deep eutectic solvents, promising green alternatives to conventional solvents, consist of a hydrogen bond donor and a hydrogen bond acceptor. The hydrogen bonding components in deep eutectic solvents form an extended hydrogen bonding network, which can be tuned to specific applications by changing the hydrogen bond donors. In this work, we have changed the hydrogen bond donor from a diol to a dicarboxylic acid by systematically replacing a hydroxyl group with an acid group one at a time to investigate the solvation structure and dynamics of the deep eutectic systems. Using a combination of ultrafast vibrational spectroscopy and molecular dynamics simulations, we compared the spectral diffusion and orientational relaxation dynamics of three deep eutectic systems using the vibrational responses of a dissolved anion. Our results indicate that although the solvation structures are marginally different across the systems, distinct differences are present in the solvent fluctuation and solute reorientation dynamics. This work provides a detailed molecular understanding of carboxylic-acid-based deep eutectic systems and how they differ from alcohol-based deep eutectic systems.

Spatial Layouts of Low-Entropy Hydration Shells Guide Protein Binding

Protein-protein binding enables orderly biological self-organization and is therefore considered a miracle of nature. Protein‒protein binding is driven by electrostatic forces, hydrogen bonding, van der Waals force, and hydrophobic interactions. Among these physical forces, only hydrophobic interactions can be considered long-range intermolecular attractions between proteins due to the electrostatic shielding of surrounding water molecules. Low-entropy hydration shells around proteins drive hydrophobic attraction among them that essentially coordinate protein‒protein binding. Here, an innovative method is developed for identifying low-entropy regions of hydration shells of proteins by screening off pseudohydrophilic groups on protein surfaces and revealing that large low-entropy regions of the hydration shells typically cover the binding sites of individual proteins. According to an analysis of determined protein complex structures, shape matching between a large low-entropy hydration shell region of a protein and that of its partner at the binding sites is revealed as a universal law. Protein‒protein binding is thus found to be mainly guided by hydrophobic collapse between the shape-matched low-entropy hydration shells that is verified by bioinformatics analyses of hundreds of structures of protein complexes, which cover four test systems. A simple algorithm is proposed to accurately predict protein binding sites.

Heterogeneous and Allosteric Role of Surface Hydration for Protein-Ligand Binding

Atomistic-level understanding of surface hydration mediating protein-protein interactions and ligand binding has been a challenge due to the dynamic nature of water molecules near the surface. We develop a computational method to evaluate the solvation free energy based on the density map of the first hydration shell constructed from all-atom molecular dynamics simulation and use it to examine the binding of two intrinsically disordered ligands to their target protein domain. One ligand is from the human protein, and the other is from the 1918 Spanish flu virus. We find that the viral ligand incurs a 6.9 kcal/mol lower desolvation penalty upon binding to the target, which is consistent with its stronger binding affinity. The difference arises from the spatially fragmented and nonuniform water density profiles of the first hydration shell. In particular, residues that are distal from the ligand-binding site contribute to a varying extent to the desolvation penalty, among which the "entropy hotspot" residues contribute significantly. Thus, ligand binding alters hydration on remote sites in addition to affecting the binding interface. The nonlocal effect disappears when the conformational motion of the protein is suppressed. The present results elucidate the interplay between protein conformational dynamics and surface hydration. Our approach of measuring the solvation free energy based on the water density of the first hydration shell is tolerant of the conformational fluctuation of protein, and we expect it to be applicable to investigating a broad range of biomolecular interfaces.

Correlation between protein conformations and water structure and thermodynamics at high pressure: A molecular dynamics study of the Bovine Pancreatic Trypsin Inhibitor (BPTI) protein

Pressure-induced perturbation of a protein structure leading to its folding-unfolding mechanism is an important yet not fully understood phenomenon. The key point here is the role of water and its coupling with protein conformations as a function of pressure. In the current work, using extensive molecular dynamics simulation at 298 K, we systematically examine the coupling between protein conformations and water structures of pressures of 0.001, 5, 10, 15, 20 kbar, starting from (partially) unfolded structures of the protein Bovine Pancreatic Trypsin Inhibitor (BPTI). We also calculate localized thermodynamics at those pressures as a function of protein-water distance. Our findings show that both protein-specific and generic effects of pressure are operating. In particular, we found that (1) the amount of increase in water density near the protein depends on the protein structural heterogeneity; (2) the intra-protein hydrogen bond decreases with pressure, while the water-water hydrogen bond per water in the first solvation shell (FSS) increases; protein-water hydrogen bonds also found to increase with pressure, (3) with pressure hydrogen bonds of waters in the FSS getting twisted; and (4) water's tetrahedrality in the FSS decreases with pressure, but it is dependent on the local environment. Thermodynamically, at higher pressure, the structural perturbation of BPTI is due to the pressure-volume work, while the entropy decreases with the increase of pressure due to the higher translational and rotational rigidity of waters in the FSS. The local and subtle effects of pressure, found in this work, are likely to be typical of pressure-induced protein structure perturbation.

Frequency Changes in Terminal Alkynes Provide Strong, Sensitive, and Solvatochromic Raman Probes of Biochemical Environments

The C≡C stretching frequencies of terminal alkynes appear in the "clear" window of vibrational spectra, so they are attractive and increasingly popular as site-specific probes in complicated biological systems like proteins, cells, and tissues. In this work, we collected infrared (IR) absorption and Raman scattering spectra of model compounds, artificial amino acids, and model proteins that contain terminal alkyne groups, and we used our results to draw conclusions about the signal strength and sensitivity to the local environment of both aliphatic and aromatic terminal alkyne C≡C stretching bands. While the IR bands of alkynyl model compounds displayed surprisingly broad solvatochromism, their absorptions were weak enough that alkynes can be ruled out as effective IR probes. The same solvatochromism was observed in model compounds' Raman spectra, and comparisons to published empirical solvent scales (including a linear regression against four meta-aggregated solvent parameters) suggested that the alkyne C≡C stretching frequency mainly reports on local electronic interactions (i.e., short-range electron donor-acceptor interactions) with solvent molecules and neighboring functional groups. The strong solvatochromism observed here for alkyne stretching bands introduces an important consideration for Raman imaging studies based on these signals. Raman signals for alkynes (especially those that are π-conjugated) can be exceptionally strong and should permit alkynyl Raman signals to function as probes at very low concentrations, as compared to other widely used vibrational probe groups like azides and nitriles. We incorporated homopropargyl glycine into a transmembrane helical peptide via peptide synthesis, and we installed p-ethynylphenylalanine into the interior of the Escherichia coli fatty acid acyl carrier protein using a genetic code expansion technique. The Raman spectra from each of these test systems indicate that alkynyl C≡C bands can act as effective and unique probes of their local biomolecular environments. We provide guidance for the best possible future uses of alkynes as solvatochromic Raman probes, and while empirical explanations of the alkyne solvatochromism are offered, open questions about its physical basis are enunciated.

Resolving the Impact of Hydrogen Bonding on the Phylloquinone Cofactor through Two-Dimensional Infrared Spectroscopy

Two-dimensional infrared spectroscopy (2DIR) was applied to phylloquinone (PhQ), an important biological cofactor, to elucidate the impact of hydrogen bonding on the ultrafast dynamics and energetics of the carbonyl stretching modes. 2DIR measurements were performed on PhQ dissolved in hexanol, which served as the hydrogen bonding solvent, and hexane, which served as a non-hydrogen bonding control. Molecular dynamics simulations and quantum chemical calculations were performed to aid in spectral assignment and interpretation. From the position of the peaks in the 2DIR spectra, we extracted the transition frequencies for the fundamental, overtone, and combination bands of hydrogen bonded and non-hydrogen bonded carbonyl groups of PhQ in the 1635-1680 cm^-1 region. We find that hydrogen bonding to a single carbonyl group acts to decouple the two carbonyl units of PhQ. Through analysis of the time-resolved 2DIR data, we find that hydrogen bonding leads to faster vibrational relaxation as well as an increase in the inhomogeneous broadening of the carbonyl groups. Overall, this work demonstrates how hydrogen bonding to the carbonyl groups of PhQ presents in the 2DIR spectra, laying the groundwork to use PhQ as a 2DIR probe to characterize the ultrafast fluctuations in the local environment of natural photosynthetic complexes.

Cross-Correlated Motions in Azidolysozyme

The changes in the local and global dynamics of azide-labelled lysozyme compared with that of the wild type protein are quantitatively assessed for all alanine residues along the polypeptide chain. Although attaching -N3 to alanine residues has been considered to be a minimally invasive change in the protein it is found that depending on the location of the alanine residue, the local and global changes in the dynamics differ. For Ala92, the change in the cross-correlated motions are minimal, whereas attaching -N3 to Ala90 leads to pronounced differences in the local and global correlations as quantified by the cross-correlation coefficients of the Cα atoms. We also demonstrate that the spectral region of the asymmetric azide stretch distinguishes between alanine attachment sites, whereas changes in the low frequency, far-infrared region are less characteristic.

Decoding the Kinetic Pathways toward a Lipid/DNA Complex of Alkyl Alcohol Cationic Lipids Formed in a Microfluidic Channel

Complexes of cationic liposomes with DNA have emerged as promising nonviral vectors for delivering genetic information into cells for gene therapy. Kinetics of the liposome/DNA complex (lipoplex) formation on a millisecond time scale are studied by monitoring time evolution of fluorescence of 8-anilino-1-naphthalene sulfonic acid (ANS) and ethidium bromide (EtBr) in a continuous flow microfluidic channel coupled to a fluorescence microscope. The formation of lipoplexes between calf thymus DNA and liposomes based on two novel cationic lipids (Lip1810 and Lip1814) are found to follow a two-step process with kinetic constants for the Lip1814/DNA complex (k₁ = 1120-1383 s^-1, k₂ = 0.227-1.45 s^-1) being significantly different from those (k₁ = 68.53-98.5 s^-1, k₂ = 32.3-60.19 s^-1) corresponding to formation of the Lip1810/DNA complex. The kinetic pathway leading to the formation of Lip1814/DNA complex is diffusion-controlled whereas the formation of Lip1810/DNA complex occurs by a conformational rearrangement-controlled pathway. The observed difference in the kinetics of lipoplex formation likely originates from different structures of the lipid/DNA complexes.

Seeking Solvation: Exploring the Role of Protein Hydration in Silk Gelation

The mechanism by which arthropods (e.g., spiders and many insects) can produce silk fibres from an aqueous protein (fibroin) solution has remained elusive, despite much scientific investigation. In this work, we used several techniques to explore the role of a hydration shell bound to the fibroin in native silk feedstock (NSF) from Bombyx mori silkworms. Small angle X-ray and dynamic light scattering (SAXS and DLS) revealed a coil size (radius of gyration or hydrodynamic radius) around 12 nm, providing considerable scope for hydration. Aggregation in dilute aqueous solution was observed above 65 °C, matching the gelation temperature of more concentrated solutions and suggesting that the strength of interaction with the solvent (i.e., water) was the dominant factor. Infrared (IR) spectroscopy indicated decreasing hydration as the temperature was raised, with similar changes in hydration following gelation by freezing or heating. It was found that the solubility of fibroin in water or aqueous salt solutions could be described well by a relatively simple thermodynamic model for the stability of the protein hydration shell, which suggests that the affected water is enthalpically favoured but entropically penalised, due to its reduced (vibrational or translational) dynamics. Moreover, while the majority of this investigation used fibroin from B. mori, comparisons with published work on silk proteins from other silkworms and spiders, globular proteins and peptide model systems suggest that our findings may be of much wider significance.

Using Molecular Simulation to Guide Protein Engineering for Biocatalysis in Organic Solvents

Biocatalysis in organic solvents (OSs) is very appealing for the industry in producing bulk and/or fine chemicals, such as pharmaceuticals, biodiesel, and fragrances. The poor performance of enzymes in OSs (e.g., reduced activity, insufficient stability, and deactivation) negates OSs' excellent solvent properties. Molecular dynamics (MD) simulations provide a complementary method to study the relationship between enzymes dynamics and the stability in OSs. Here we describe computational procedure for MD simulation of enzymes in OSs with an example of Bacillus subtilis lipase A (BSLA) in dimethyl sulfoxide (DMSO) cosolvent with software GROMACS. We discuss main essential practical issues considered (such as choice of force field, parameterization, simulation setup, and trajectory analysis). The core part of this protocol (enzyme-OS system setup, analysis of structural-based and solvation-based observables) is transferable to other enzymes and any OS systems. Combining with experimental studies, the obtained molecular knowledge is most likely to guide researchers to access rational protein engineering approaches to tailor OS resistant enzymes and expand the scope of biocatalysis in OS media. Finally, we discuss potential solutions to overcome the remaining challenges of computational biocatalysis in OSs and briefly draw future directions for further improvement in this field.

Universal dynamical onset in water at distinct material interfaces

Interfacial water remains liquid and mobile much below 0 °C, imparting flexibility to the encapsulated materials to ensure their diverse functions at subzero temperatures. However, a united picture that can describe the dynamical differences of interfacial water on different materials and its role in imparting system-specific flexibility to distinct materials is lacking. By combining neutron spectroscopy and isotope labeling, we explored the dynamics of water and the underlying substrates independently below 0 °C across a broad range of materials. Surprisingly, while the function-related anharmonic dynamical onset in the materials exhibits diverse activation temperatures, the surface water presents a universal onset at a common temperature. Further analysis of the neutron experiment and simulation results revealed that the universal onset of water results from an intrinsic surface-independent relaxation: switching of hydrogen bonds between neighboring water molecules with a common energy barrier of ∼35 kJ mol⁻¹.

We demonstrated that the dynamical onset of interfacial water is an intrinsic property of water itself, resulting from a surface independent relaxation process in water with an approximately universal energy barrier of ∼35 kJ mol⁻¹.

Site-Specific Interrogation of Protein Structure and Stability

To execute their function or activity, proteins need to possess variability in local electrostatic environment, solvent accessibility, structure, and stability. However, assessing any protein property in a site-specific manner is not easy since native spectroscopic signals often lack the needed specificity. One strategy that overcomes this limitation is to use unnatural amino acids that exhibit distinct spectroscopic features. In this chapter, we describe several such unnatural amino acids (UAAs) and their respective applications in site-specific interrogation of protein structure and stability using standard biophysical methods, including circular dichroism (CD), infrared (IR), and fluorescence spectroscopies.

Protein-Bath Coupling of an Internal Reaction Coordinate at Intermediate Time Scales

Thermally activated barrier-crossing processes are central to protein reaction kinetics. A determining factor for such kinetics is the extent to which the protein's motions are coupled to the surrounding bath. It is understood that slow large-scale conformational motions are strongly coupled to the environment, while fast librational motions are uncoupled. However, less is known about protein-bath coupling of reaction coordinates located on the interior of a protein and with dynamics on intermediate time scales. In this work, we use single molecule 2D fluorescence lifetime correlation spectroscopy to study the microsecond chemical reaction occurring in the chromophore pocket of eGFP. The equilibrium reaction involves a dihedral rotation of a glutamic acid residue and a rearrangement of the local hydrogen-bonding network surrounding the endogenous chromophore, with no accompanying large-scale conformational changes. We observe that the internal chemical reaction is coupled to the solvent viscosity, though the scaling deviates from Kramers' behavior. We attribute this deviation to the internal friction of the protein, which weakens the protein-solvent coupling at high viscosity and intermediate time scales.

Ultrafast vibrational dynamics of a solute correlates with dynamics of the solvent

Two-dimensional infrared (2D-IR) spectroscopy is used to measure the spectral dynamics of the metal carbonyl complex cyclopentadienyl manganese tricarbonyl (CMT) in a series of linear alkyl nitriles. 2D-IR spectroscopy provides direct readout of solvation dynamics through spectral diffusion, probing the decay of frequency correlation induced by fluctuations of the solvent environment. 2D-IR simultaneously monitors intramolecular vibrational energy redistribution (IVR) among excited vibrations, which can also be influenced by the solvent through the spectral density rather than the dynamical friction underlying solvation. Here, we report that the CMT vibrational probe reveals solvent dependences in both the spectral diffusion and the IVR time scales, where each slows with increased alkyl chain length. In order to assess the degree to which solute-solvent interactions can be correlated with bulk solvent properties, we compared our results with low-frequency dynamics obtained from optical Kerr effect (OKE) spectroscopy-performed by others-on the same nitrile solvent series. We find excellent correlation between our spectral diffusion results and the orientational dynamics time scales from OKE. We also find a correlation between our IVR time scales and the amplitudes of the low-frequency spectral densities evaluated at the 90-cm^-1 energy difference, corresponding to the gap between the two strong vibrational modes of the carbonyl probe. 2D-IR and OKE provide complementary perspectives on condensed phase dynamics, and these findings provide experimental evidence that at least at the level of dynamical correlations, some aspects of a solute vibrational dynamics can be inferred from properties of the solvent.

The Structural and Dynamical Properties of the Hydration of SNase Based on a Molecular Dynamics Simulation

The dynamics of protein–water fluctuations are of biological significance. Molecular dynamics simulations were performed in order to explore the hydration dynamics of staphylococcal nuclease (SNase) at different temperatures and mutation levels. A dynamical transition in hydration water (at ~210 K) can trigger larger-amplitude fluctuations of protein. The protein–water hydrogen bonds lost about 40% in the total change from 150 K to 210 K, while the Mean Square Displacement increased by little. The protein was activated when the hydration water in local had a comparable trend in making hydrogen bonds with protein– and other waters. The mutations changed the local chemical properties and the hydration exhibited a biphasic distribution, with two time scales. Hydrogen bonding relaxation governed the local protein fluctuations on the picosecond time scale, with the fastest time (24.9 ps) at the hydrophobic site and slowest time (40.4 ps) in the charged environment. The protein dynamic was related to the water’s translational diffusion via the relaxation of the protein–water’s H-bonding. The structural and dynamical properties of protein–water at the molecular level are fundamental to the physiological and functional mechanisms of SNase.

Transparent window 2D IR spectroscopy of proteins

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

Protein Dynamics by Two-Dimensional Infrared Spectroscopy

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

Direct comparison of amplitude and geometric measures of spectral inhomogeneity using phase-cycled 2D-IR spectroscopy

Two-dimensional infrared (2D-IR) spectroscopy provides access to equilibrium dynamics with the extraction of the frequency-fluctuation correlation function (FFCF) from the measured spectra. Several different methods of obtaining the FFCF from experimental spectra, such as the center line slope (CLS), ellipticity, phase slope, and nodal line slope, all depend on the geometrical nature of the 2D line shape and necessarily require spectral extent in order to achieve a measure of the FFCF. Amplitude measures, on the other hand, such as the inhomogeneity index, rely only on signal amplitudes and can, in principle, be computed using just a single point in a 2D spectrum. With a pulse shaper-based 2D-IR spectrometer, in conjunction with phase cycling, we separate the rephasing and nonrephasing signals used to determine the inhomogeneity index. The same measured data provide the absorptive spectrum, needed for the CLS. Both methods are applied to two model molecular systems: tungsten hexacarbonyl (WCO₆) and methylcyclopentadienyl manganese tricarbonyl [Cp'Mn(CO)₃, MCMT]. The three degenerate IR modes of W(CO)₆ lack coherent modulation or noticeable intramolecular vibrational redistribution (IVR) and are used to establish a baseline comparison. The two bands of the MCMT tripod complex include intraband coherences and IVR as well as likely internal torsional motion on a few-picosecond time scale. We find essentially identical spectral diffusion, but faster, non-equilibrium dynamics lead to differences in the FFCFs extracted with the two methods. The inhomogeneity index offers an advantage in cases where spectra are complex and energy transfer can mimic line shape changes due to frequency fluctuations.

Site-selective dynamics of azidolysozyme

The spectroscopic response of and structural dynamics around all azido-modified alanine residues (AlaN₃) in lysozyme are characterized. It is found that AlaN₃ is a positionally sensitive probe for the local dynamics, covering a frequency range of ∼15 cm^-1 for the center frequency of the line shape. This is consistent with findings from selective replacements of amino acids in PDZ2, which reported a frequency span of ∼10 cm^-1 for replacements of Val, Ala, or Glu by azidohomoalanine. For the frequency fluctuation correlation functions, the long-time decay constants τ₂ range from ∼1 to ∼10 ps, which compares with experimentally measured correlation times of 3 ps. Attaching azide to alanine residues can yield dynamics that decays to zero on the few ps time scale (i.e., static component Δ₀ ∼ 0 ps^-1) or to a remaining, static contribution of ∼0.5 ps^-1 (corresponding to 2.5 cm^-1), depending on the local environment on the 10 ps time scale. The magnitude of the static component correlates qualitatively with the degree of hydration of the spectroscopic probe. Although attaching azide to alanine residues is found to be structurally minimally invasive with respect to the overall protein structure, analysis of the local hydrophobicity indicates that the hydration around the modification site differs for modified and unmodified alanine residues, respectively.

Protein flexibility reduces solvent-mediated friction barriers of ligand binding to a hydrophobic surface patch

Solvent fluctuations have been explored in detail for idealized and rigid hydrophobic model systems, but so far it has remained unclear how internal protein motions and their coupling to the surrounding solvent affect the dynamics of ligand binding to biomolecular surfaces. Here, molecular dynamics simulations were used to elucidate the solvent-mediated binding of a model ligand to the hydrophobic surface patch of ubiquitin. The ligand's friction profiles reveal pronounced long-time correlations and enhanced friction in the vicinity of the protein, similar to idealized hydrophobic surfaces. Interestingly, these effects are shaped by internal protein motions. Protein flexibility modulates water density fluctuations near the hydrophobic surface patch and smooths out the friction profile of ligand binding.

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.

Mapping solvation heterogeneity in live cells by hyperspectral stimulated Raman scattering microscopy

Water provides a dynamic matrix in which all biochemical processes occur in living organisms. The structure and dynamics of intracellular water constitute the cornerstone for understanding all aspects of cellular function. Fundamentally, direct visualization of subcellular solvation heterogeneity is essential but remains challenging with commonly used nuclear magnetic resonance methods due to poor spatial resolution. To explore this question, we demonstrate a vibrational-shift imaging approach by combining the spectral-focusing hyperspectral stimulated Raman scattering technique with an environmentally sensitive nitrile probe. The sensing ability of a near-infrared nitrile-containing molecule is validated in the solution phase, microscopic droplets, and cellular environments. Finally, we quantitatively measure the subcellular solvation variance between the cytoplasm (29.5%, S.E. 1.8%) and the nucleus (57.3%, S.E. 1.0%), which is in good agreement with previous studies. This work sheds light on heterogeneous solvation in live systems using coherent Raman microscopy and opens up new avenues to explore environmental variance in complex systems with high spatiotemporal resolution.

Decoupling between the translation and rotation of water in the proximity of a protein molecule

The interaction between water and biomacromolecules is of fundamental interest in biophysics, biochemistry and physical chemistry. By combining neutron scattering and molecular dynamics simulations on a perdeuterated protein at a series of hydration levels, we demonstrated that the translational motion of water is slowed down more significantly than its rotation, when water molecules approach the protein molecule. Further analysis of the simulation trajectories reveals that the observed decoupling results from the fact that the translational motion of water is more correlated over space and more retarded by the charged/polar residues and spatial confinement on the protein surface, than the rotation. Moreover, around the stable protein residues (with smaller atomic fluctuations), water exhibits more decoupled dynamics, indicating a connection between the observed translation-rotation decoupling in hydration water and the local stability of the protein molecule.

Hydrophilic dangling chain interfacial segregation in polyurethane networks at aqueous interfaces and its underlying mechanisms: molecular dynamics simulations

Polymer networks with hydrophilic dangling chains are ideal candidates for many submerged applications, e.g., protein non-adhesive coatings with non-fouling behavior. The dangling chains segregate from the polymer network towards the water and form a brush-like structure at the interface. Several factors such as the polymer network structure, dangling chain length, and water/dangling chain interaction may all affect the interfacial performance of the polymer. Therefore, we employed a Martini based coarse-grained (CG) molecular dynamics (MD) simulation to elucidate the influences of the abovementioned parameters on dangling chain interfacial segregation. We built up several polyurethane (PU) networks based on poly(tetra methylene glycol) (PTMG), as a macrodiol, and methoxy poly(ethylene glycol) (mPEG), as a dangling chain, with varying molecular weights. We found out that the macrodiol/dangling chain length ratio considerably smaller than one impedes the migration of dangling chains towards the water interface, while the dangling chain hydrophilicity and length determine the polymer interfacial layer density/thickness. Then, we artificially changed the dangling chain affinity to water from an intermediate to a very attractive water/dangling chain interaction. We justified that a brush-like structure forms in two consecutive steps: first, a longitudinal, and then a lateral migration of dangling chains in water. The latter step results in a uniform interfacial layer over the polymer interface that mainly occurs in the case of the attractive water/dangling chain interaction.

Phonon-mediated lipid raft formation in biological membranes

Short-wavelength collective molecular motions, also known as phonons, have recently attracted much interest in revealing dynamic properties of biological membranes through the use of neutron and X-ray scattering, infrared and Raman spectroscopies, and molecular dynamics simulations. Experimentally detecting unique vibrational patterns such as, shear phonon excitations, viscoelastic crossovers, transverse acoustic phonon gaps, and continuous and truncated optical phonon modes in cellular membranes, to name a few, has proven non-trivial. Here, we review recent advances in liquid thermodynamics that have resulted in the development of the phonon theory of liquids. The theory has important predictions regarding the shear vibrational spectra of fluids, namely the emergence of viscoelastic crossovers and transverse acoustic phonon gaps. Furthermore, we show that these vibrational patterns are common in soft (non-crystalline) materials, including, but not limited to liquids, colloids, liquid crystals (mesogens), block copolymers, and biological membranes. The existence of viscoelastic crossovers and acoustic phonon gaps define the self-diffusion properties of cellular membranes and provide a molecular picture of the transient nature of lipid rafts (Bolmatov et al., 2020). Importantly, the timescales (picoseconds) for the formation and dissolution of transient lipid rafts match the lifetime of the formation and breakdown of interfacial water hydrogen bonds. Apart from acoustic propagating phonon modes, biological membranes can also support more energetic non-propagating optical phonon excitations, also known as standing waves or breathing modes. Importantly, optical phonons can be truncated due to the existence of finite size nanodomains made up of strongly correlated lipid-cholesterol molecular pairs. These strongly coupled molecular pairs can serve as nucleation centers for the formation of stable rafts at larger length scales, due to correlations of spontaneous fluctuations (Onsager's regression hypothesis). Finally and importantly, molecular level viscoelastic crossovers, acoustic phonon gaps, and continuous and truncated optical phonon modes may offer insights as to how lipid-lipid and lipid-protein interactions enable biological function.

Protein Preferential Solvation in Water:Glycerol Mixtures

For proteins in solvent mixtures, the relative abundances of each solvent in their solvation shell have a critical impact on their properties. Preferential solvation of a series of proteins in water-glycerol mixtures is studied here over a broad range of solvent compositions via classical molecular dynamics simulations. Our simulation results reveal that the differences between shell and bulk compositions exhibit dramatic changes with solvent composition, temperature, and protein nature. In contrast with the simple and widely used picture where glycerol is completely excluded from the protein interface, we show that for aqueous solutions with less than 50% glycerol in volume, protein solvation shells have approximately the same composition as the bulk solvent and proteins are in direct contact with glycerol. We further demonstrate that at high glycerol concentration, glycerol depletion from the solvation shell is due to an entropic factor arising from the reduced accessibility of bulky glycerol molecules in protein cavities. The resulting molecular picture is important to understand protein activity and cryopreservation in mixed aqueous solvents.

Dipole-dipole interactions between tryptophan side chains and hydration water molecules dominate the observed dynamic stokes shift of lysozyme

The fluorescence intensity of tryptophan residues in hen egg-white lysozyme was measured up to 500 ps after the excitation by irradiation pulses at 290 nm. From the time-dependent variation of fluorescence intensity in a wavelength range of 320-370 nm, the energy relaxation in the dynamic Stokes shift was reconstructed as the temporal variation in wavenumber of the estimated fluorescence maximum. The relaxation was approximated by two exponential curves with decay constants of 1.2 and 26.7 ps. To interpret the relaxation, a molecular dynamics simulation of 75 ns was conducted for lysozyme immersed in a water box. From the simulation, the energy relaxation in the electrostatic interactions of each tryptophan residue was evaluated by using a scheme derived from the linear response theory. Dipole-dipole interactions between each of the Trp62 and Trp123 residues and hydration water molecules displayed an energy relaxation similar to that experimentally observed regarding time constants and magnitudes. The side chains of these residues were partly or fully exposed to the solvent. In addition, by inspecting the variation in dipole moments of the hydration water molecules around lysozyme, it was suggested that the observed relaxation could be attributed to the orientational relaxation of hydration water molecules participating in the hydrogen-bond network formed around each of the two tryptophan residues.

Protein Solvent Shell Structure Provides Rapid Analysis of Hydration Dynamics

The solvation layer surrounding a protein is clearly an intrinsic part of protein structure-dynamics-function, and our understanding of how the hydration dynamics influences protein function is emerging. We have recently reported simulations indicating a correlation between regional hydration dynamics and the structure of the solvation layer around different regions of the enzyme Candida antarctica lipase B, wherein the radial distribution function (RDF) was used to calculate the pairwise entropy, providing a link between dynamics (diffusion) and thermodynamics (excess entropy) known as Rosenfeld scaling. Regions with higher RDF values/peaks in the hydration layer (the first peak, within 6 Å of the protein surface) have faster diffusion in the hydration layer. The finding thus hinted at a handle for rapid evaluation of hydration dynamics at different regions on the protein surface in molecular dynamics simulations. Such an approach may move the analysis of hydration dynamics from a specialized venture to routine analysis, enabling an informatics approach to evaluate the role of hydration dynamics in biomolecular function. This paper first confirms that the correlation between regional diffusive dynamics and hydration layer structure (via water center of mass around protein side-chain atom RDF) is observed as a general relationship across a set of proteins. Second, it seeks to devise an approach for rapid analysis of hydration dynamics, determining the minimum amount of information and computational effort required to get a reliable value of hydration dynamics from structural data in MD simulations based on the protein-water RDF. A linear regression model using the integral of the hydration layer in the water-protein RDF was found to provide statistically equivalent apparent diffusion coefficients at the 95% confidence level for a set of 92 regions within five different proteins. In summary, RDF analysis of 10 ns of data after simulation convergence is sufficient to accurately map regions of fast and slow hydration dynamics around a protein surface. Additionally, it is anticipated that a quick look at protein-water RDFs, comparing peak heights, will be useful to provide a qualitative ranking of regions of faster and slower hydration dynamics at the protein surface for rapid analysis when investigating the role of solvent dynamics in protein function.

Site-Specific 1D and 2D IR Spectroscopy to Characterize the Conformations and Dynamics of Protein Molecular Recognition

Proteins exist as ensembles of interconverting states on a complex energy landscape. A complete, molecular-level understanding of their function requires knowledge of the populated states and thus the experimental tools to characterize them. Infrared (IR) spectroscopy has an inherently fast time scale that can capture all states and their dynamics with, in principle, bond-specific spatial resolution, and 2D IR methods that provide richer information are becoming more routine. Although application of IR spectroscopy for investigation of proteins is challenged by spectral congestion, the issue can be overcome by site-specific introduction of amino acid side chains that have IR probe groups with frequency-resolved absorptions, which furthermore enables selective characterization of different locations in proteins. Here, we briefly introduce the biophysical methods and summarize the current progress toward the study of proteins. We then describe our efforts to apply site-specific 1D and 2D IR spectroscopy toward elucidation of protein conformations and dynamics to investigate their involvement in protein molecular recognition, in particular mediated by dynamic complexes: plastocyanin and its binding partner cytochrome f, cytochrome P450s and substrates or redox partners, and Src homology 3 domains and proline-rich peptide motifs. We highlight the advantages of frequency-resolved probes to characterize specific, local sites in proteins and uncover variation among different locations, as well as the advantage of the fast time scale of IR spectroscopy to detect rapidly interconverting states. In addition, we illustrate the greater insight provided by 2D methods and discuss potential routes for further advancement of the field of biomolecular 2D IR spectroscopy.

Ultrafast energy relaxation dynamics of amide I vibrations coupled with protein-bound water molecules

The influence of hydration water on the vibrational energy relaxation in a protein holds the key to understand ultrafast protein dynamics, but its detection is a major challenge. Here, we report measurements on the ultrafast vibrational dynamics of amide I vibrations of proteins at the lipid membrane/H₂O interface using femtosecond time-resolved sum frequency generation vibrational spectroscopy. We find that the relaxation time of the amide I mode shows a very strong dependence on the H₂O exposure, but not on the D₂O exposure. This observation indicates that the exposure of amide I bond to H₂O opens up a resonant relaxation channel and facilitates direct resonant vibrational energy transfer from the amide I mode to the H₂O bending mode. The protein backbone motions can thus be energetically coupled with protein-bound water molecules. Our findings highlight the influence of H₂O on the ultrafast structure dynamics of proteins.

Heterogeneous and Highly Dynamic Interface in Plastocyanin-Cytochrome f Complex Revealed by Site-Specific 2D-IR Spectroscopy

Transient protein complexes are crucial for sustaining dynamic cellular processes. The complexes of electron-transfer proteins are a notable example, such as those formed by plastocyanin (Pc) and cytochrome f (cyt f) in the photosynthetic apparatus. The dynamic and heterogeneous nature of these complexes, however, makes their study challenging. To better elucidate the complex of Nostoc Pc and cyt f, 2D-IR spectroscopy coupled to site-specific labeling with cyanophenylalanine infrared (IR) probes was employed to characterize how the local environments at sites along the surface of Pc were impacted by cyt f binding. The results indicate that Pc most substantially engages with cyt f via the hydrophobic patch around the copper redox site. Complexation with cyt f led to an increase in inhomogeneous broadening of the probe absorptions, reflective of increased heterogeneity of interactions with their environment. Notably, most of the underlying states interconverted very rapidly (1 to 2 ps), suggesting a complex with a highly mobile interface. The data support a model of the complex consisting of a large population of an encounter complex. Additionally, the study demonstrates the application of 2D-IR spectroscopy with site-specifically introduced probes to reveal new quantitative insight about dynamic biochemical systems.

Site-specific 2D IR spectroscopy: a general approach for the characterization of protein dynamics with high spatial and temporal resolution

The conformational heterogeneity and dynamics of protein side chains contribute to function, but investigating exactly how is hindered by experimental challenges arising from the fast timescales involved and the spatial heterogeneity of protein structures. The potential of two-dimensional infrared (2D IR) spectroscopy for measuring conformational heterogeneity and dynamics with unprecedented spatial and temporal resolution has motivated extensive effort to develop amino acids with functional groups that have frequency-resolved absorptions to serve as probes of their protein microenvironments. We demonstrate the full advantage of the approach by selective incorporation of the probe p-cyanophenylalanine at six distinct sites in a Src homology 3 domain and the application of 2D IR spectroscopy to site-specifically characterize heterogeneity and dynamics and their contribution to cognate ligand binding. The approach revealed a wide range of microenvironments and distinct responses to ligand binding, including at the three adjacent, conserved aromatic residues that form the recognition surface of the protein. Molecular dynamics simulations performed for all the labeled proteins provide insight into the underlying heterogeneity and dynamics. Similar application of 2D IR spectroscopy and site-selective probe incorporation will allow for the characterization of heterogeneity and dynamics of other proteins, how heterogeneity and dynamics are affected by solvation and local structure, and how they might contribute to biological function.

Overhauser Dynamic Nuclear Polarization for the Study of Hydration Dynamics, Explained

We outline the physical properties of hydration water that are captured by Overhauser Dynamic Nuclear Polarization (ODNP) relaxometry and explore the insights that ODNP yields about the water and the surface that this water is coupled to. As ODNP relies on the pairwise cross-relaxation between the electron spin of a spin probe and a proton nuclear spin of water, it captures the dynamics of single-particle diffusion of an ensemble of water molecules moving near the spin probe. ODNP principally utilizes the same physics as other nuclear magnetic resonance (NMR) relaxometry (i.e., relaxation measurement) techniques. However, in ODNP, electron paramagnetic resonance (EPR) excites the electron spins probes and their high net polarization acts as a signal amplifier. Furthermore, it renders ODNP parameters highly sensitive to water moving at rates commensurate with the EPR frequency of the spin probe (typically 10GHz). Also, ODNP selectively enhances the NMR signal contributions of water moving within close proximity to the spin label. As a result, ODNP can capture ps-ns movements of hydration waters with high sensitivity and locality, even in samples with protein concentrations as dilute as 10 µM. To date, the utility of the ODNP technique has been demonstrated for two major applications: the characterization of the spatial variation in the properties of the hydration layer of proteins or other surfaces displaying topological diversity, and the identification of structural properties emerging from highly disordered proteins and protein domains. The former has been shown to correlate well with the properties of hydration water predicted by MD simulations and has been shown capable of evaluating the hydrophilicity or hydrophobicity of a surface. The latter has been demonstrated for studies of an interhelical loop of proteorhodopsin, the partial structure of α-synuclein embedded at the lipid membrane surface, incipient structures adopted by tau proteins en route to fibrils, and the structure and hydration profile of a transmembrane peptide. This chapter focuses on offering a mechanistic understanding of the ODNP measurement and the molecular dynamics encoded in the ODNP parameters. In particular, it clarifies how the electron-nuclear dipolar coupling encodes information about the molecular dynamics in the nuclear spin self-relaxation and, more importantly, the electron-nuclear spin cross-relaxation rates. The clarification of the molecular dynamics underlying ODNP should assist in establishing a connection to theory and computer simulation that will offer far richer interpretations of ODNP results in future studies.