Effect of solvent polarity on the vibrational dephasing dynamics of the nitrosyl stretch in an Fe(II) complex revealed by 2D IR spectroscopy

Reply to "Comment on: 'Isolating Vibrational Polariton 2D-IR Transmission Spectra'"

In a Comment on our recent Letter, the authors take issue with our method of refining 2D-IR transmission spectra to remove a background contribution that arises from nonpolaritonic molecules in the cavity. In our response to their Comment, we describe how our approach was motivated by the previous work of the authors, and we present a spatially dependent molecule-cavity Tavis-Cummings model that can account for the significant response from localized molecules with nonzero oscillator strengths. The telltale signature of the localized molecule response is the spectral diffusion dynamics of the bare W(CO)₆ molecules in the polar butyl acetate solvent. Inhomogeneous broadening is absent from polaritonic states due to the extreme degree of exchange narrowing in coupling very large numbers of molecules to a cavity mode.

Isolating Polaritonic 2D-IR Transmission Spectra

Strong coupling between vibrational transitions in molecules within a resonant optical microcavity leads to the formation of collective, delocalized vibrational polaritons. There are many potential applications of "polaritonic chemistry", ranging from modified chemical reactivity to quantum information processing. One challenge in obtaining the polaritonic response is removing a background contribution due to the uncoupled molecules that generate an ordinary 2D-IR spectrum whose amplitude is filtered by the polariton transmission spectrum. We show that most features in 2D-IR spectra of vibrational polaritons can be explained by a linear superposition of this background signal and the true polariton response. Through a straightforward correction procedure, in which the filtered bare-molecule 2D-IR spectrum is subtracted from the measured cavity response, we recover the polaritonic spectrum.

Probing Ligand Effects on the Ultrafast Dynamics of Copper Complexes via Midinfrared Pump-Probe and 2DIR Spectroscopies

The effects of ligand structural variation on the ultrafast dynamics of a series of copper coordination complexes were investigated using polarization-dependent mid-IR pump-probe spectroscopy and two-dimensional infrared (2DIR) spectroscopy. The series consists of three copper complexes [(^R3P₃tren)Cu^IIN₃]BAr₄^F (1^PR3, ^R3P₃tren = tris[2-(phosphiniminato)ethyl]amine, BAr₄^F = tetrakis(pentafluorophenyl)borate) where the number of methyl and phenyl groups in the PR₃ ligand are systematically varied across the series (PR₃ = PMe₃, PMe₂Ph, PMePh₂). The asymmetric stretching mode of azide in the 1^PR3 series is used as a vibrational probe of the small-molecule binding site. The results of the pump-probe measurements indicate that the vibrational energy of azide dissipates through intramolecular pathways and that the bulkier phenyl groups lead to an increase in the spatial restriction of the diffusive reorientation of bound azide. From 2DIR experiments, we characterize the spectral diffusion of the azide group and find that an increase in the number of phenyl groups maps to a broader inhomogeneous frequency distribution (Δ₂). This indicates that an increase in the steric bulk of the secondary coordination sphere acts to create more distinct configurations in the local environment that are accessible to the azide group. This work demonstrates how ligand structural variation affects the ultrafast dynamics of a small molecular group bound to the metal center, which could provide insight into the structure-function relationship of the copper coordination complexes and transition-metal coordination complexes in general.

Excited-state vibration-polariton transitions and dynamics in nitroprusside

Strong cavity coupling to molecular vibrations creates vibration-polaritons capable of modifying chemical reaction kinetics, product branching ratios, and charge transfer equilibria. However, the mechanisms impacting these molecular processes remain elusive. Furthermore, even basic elements determining the spectral properties of polaritons, such as selection rules, transition moments, and lifetimes are poorly understood. Here, we use two-dimensional infrared and filtered pump–probe spectroscopy to report clear spectroscopic signatures and relaxation dynamics of excited vibration-polaritons formed from the cavity-coupled NO band of nitroprusside. We apply an extended multi-level quantum Rabi model that predicts transition frequencies and strengths that agree well with our experiment. Notably, the polariton features decay ~3–4 times slower than the polariton dephasing time, indicating that they support incoherent population, a consequence of their partial matter character.

Here the authors report spectroscopy and dynamics of cavity coupled NO band of sodium nitroprusside using 2D infrared and transient spectroscopy employing pump-probe technique. They find signatures of third-order nonlinearity, incoherent and strong coupling effects of vibrational polaritons.

Linear and Nonlinear Infrared Spectroscopies Reveal Detailed Solute-Solvent Dynamic Interactions of a Nitrosyl Ruthenium Complex in Solution

In this work, the solvation of a nitrosyl ruthenium complex, [(CH₃)₄N][RuCl₃(qn)(NO)] (with qn = deprotonated 8-hydroxyquinoline), which is a potential NO-releasing molecule in the bio-environment, was studied in two bio-friendly solvents, namely deuterated dimethyl sulfoxide (dDMSO) and water (D₂O). A blue-shifted NO stretching frequency was observed in water with respect to that in dDMSO, which was believed to be due to ligand-solvent hydrogen-bonding interactions, one N═O···D and particularly three Ru-Cl···D, that show competing effects on the NO bond length. The dynamic differences of the NO stretch in these two solvents were further revealed by transient pump-probe IR and two-dimensional IR results: faster vibrational relaxation and faster spectral diffusion (SD) were observed in D₂O, confirming stronger solvent-solute interaction and also faster solvent structural dynamics in D₂O than in DMSO. Further, a significant non-decaying residual in the SD dynamics was observed in D₂O but not in DMSO, suggesting the formation of a stable solvation shell in water due to strong multi-site ligand-solvent hydrogen-bonding interactions, which is in agreement with the observed blue-shifted NO stretching frequency. This work demonstrates that small solvent molecules such as water can form a relatively rigid solvation shell for certain transition metal complexes due to cooperative ligand-solvent interactions and show slower dynamics.

Amide I SFG Spectral Line Width Probes the Lipid-Peptide and Peptide-Peptide Interactions at Cell Membrane In Situ and in Real Time

The balance of lipid-peptide and peptide-peptide interactions at cell membrane is essential to a large variety of cellular processes. In this study, we have experimentally demonstrated for the first time that sum frequency generation vibrational spectroscopy can be used to probe the peptide-peptide and lipid-peptide interactions in cell membrane in situ and in real time by determination of the line width of amide I band of protein backbone. Using a "benchmark" model of α-helical WALP23, it is found that the dominated lipid-peptide interaction causes a narrow line width of the amide I band, whereas the peptide-peptide interaction can markedly broaden the line width. When WALP23 molecules insert into the lipid bilayer, a quite narrow line width of the amide I band is observed because of the lipid-peptide interaction. In contrast, when the peptide lies down on the bilayer surface, the line width of amide I band becomes very broad owing to the peptide-peptide interaction. In terms of the real-time change in the line width, the transition from peptide-peptide interaction to lipid-peptide interaction is monitored during the insertion of WALP23 into 1,2-dipalmitoyl- sn-glycero-3-phospho-(1'- rac-glycerol) (DPPG) lipid bilayer. The dephasing time of a pure α-helical WALP23 in 1-palmitoyl-2-oleoyl- sn-glycero-3-phospho-(1'- rac-glycerol) and DPPG bilayer is determined to be 2.2 and 0.64 ps, respectively. The peptide-peptide interaction can largely accelerate the dephasing time.

Solvent exchange in preformed photocatalyst-donor precursor complexes determines efficiency

In homogeneous photocatalytic reduction of CO2, it is widely assumed that the primary electron transfer from the sacrificial donor to the catalyst is diffusion controlled, thus little attention has been paid to optimizing this step. We present spectroscopic evidence that the precursor complex is preformed, driven by preferential solvation, and two-dimensional infrared spectroscopy reveals triethanolamine (donor)/tetrahydrofuran (solvent) exchange in the photocatalyst's solvation shell, reaching greatest magnitude at the known optimal concentration (∼20% v/v TEOA in THF) for catalytically reducing CO2 to CO. Transient infrared absorption shows the appearance of the singly reduced catalyst on an ultrafast (<70 ps) time scale, consistent with non-diffusion controlled electron transfer within the preformed precursor complex. Identification of preferential catalyst-cosolvent interactions suggests a revised paradigm for the primary electron transfer, while illuminating the pivotal importance of solvent exchange in determining the overall efficiency of the photocycle.

Ultrafast kinetics of linkage isomerism in Na2[Fe(CN)5NO] aqueous solution revealed by time-resolved photoelectron spectroscopy

The kinetics of ultrafast photoinduced structural changes in linkage isomers is investigated using Na2[Fe(CN)5NO] as a model complex. The buildup of the metastable side-on configuration of the NO ligand, as well as the electronic energy levels of ground, excited, and metastable states, has been revealed by means of time-resolved extreme UV (XUV) photoelectron spectroscopy in aqueous solution, aided by theoretical calculations. Evidence of a short-lived intermediate state in the isomerization process and its nature are discussed, finding that the complete isomerization process occurs in less than 240 fs after photoexcitation.

Differentiating Two Nitrosylruthenium Isomeric Complexes by Steady-State and Ultrafast Infrared Spectroscopies

The [Ru(II)-NO⁺] group affects the structure and chemical reactivity of nitrosylruthenium(II) complexes. A characteristic infrared absorption band due to the nitrosyl (NO) stretching motion is shown in the frequency region 1800-1900 cm^-1. In this work, linear infrared (IR) and nonlinear IR methods, including pump-probe and two-dimensional (2D) IR, were utilized to study the structures and dynamics of two isomeric nitrosylruthenium complexes [Ru(OAc)(2mqn)₂NO] (H2mqn = 2-methyl-8-quinolinol) in cis and trans isomeric configurations in a weak polar solvent (CDCl₃). Using the NO stretching mode as a vibrational probe, information about local structural dynamics of the Ru complex as well as solvent fluctuation dynamics was obtained. In particular, a "structured" solvent environment is believed to form in the vicinity of the NO group in the case of the cis isomer with the aid of a neighboring OAc ligand, which is the reason for more efficient vibrational relaxation but more inhomogeneously distributed solvent and thus associated slower spectral diffusion. Our results also suggest a more anharmonic potential surface for the NO stretching mode in the less stable trans isomer.

High Sensitivity Nanosecond Mid-Infrared Transient Absorption Spectrometer Enabling Low Excitation Density Measurements of Electronic Materials

A dispersive nanosecond transient absorption instrument was developed to enable rapid time-resolved and steady-state measurements in the mid-infrared (mid-IR) region for thin films without the need for gated integrators or lock-in amplifiers. Two detectors are used depending on the experimental needs (100 MHz and 16 MHz) with time resolution from nano-millisecond and spectral coverage from 1000-5000 cm^-1 (2000-10 000 nm). The instrument utilizes flexible digitization resolution (8 bit to 14 bit) to enable high sensitivity (10^-5) measurements on thin films under low excitation (<50 µJ/cm²). We highlight the instrument's improvement over prior state-of-the-art time-resolved capabilities by measuring transient species (e.g., polarons) under extremely low energy densities (<5 µJ/cm²) in less than 10 minutes to achieve high fidelity signals. Additionally, to highlight the spectral capabilities we study two optoelectronic materials for which we resolve vibrational features as small as 10 µOD.

Structural dynamics of nitrosylruthenium isomeric complexes studied with steady-state and transient pump-probe infrared spectroscopies

The characteristic nitrosyl stretching (NO) in the region of 1800-1900cm(-1) was used to study the geometric and ligand effect on two nitrosylruthenium complexes, namely [Ru(OAc)(2QN)2NO] (QN=2-chloro-8-quinolinol (H2cqn) or QN=2-methyl-8-quinolinol (H2mqn)). The NO stretching frequency (νNO) was found in the following order: νcis-1 (2cqn)>νcis-2 (2cqn)>νcis-1 (2mqn)>νtrans (2mqn). The results exhibited a spectral sensitivity of the NO mode to both charge distribution and ligand arrangement, which was supported by ab initio computations and natural bond orbital (NBO) analyses. Further, the vibrational population of the vibrationally excited NO stretching mode was found to relax on the order of 7-10ps, showing less than 30% variation from one isomer to another, which were explained on the basis of NO local structures and solute-solvent interactions in these isomeric nitrosylruthenium complexes.

Two-dimensional vibrational-electronic spectroscopy

Two-dimensional vibrational-electronic (2D VE) spectroscopy is a femtosecond Fourier transform (FT) third-order nonlinear technique that creates a link between existing 2D FT spectroscopies in the vibrational and electronic regions of the spectrum. 2D VE spectroscopy enables a direct measurement of infrared (IR) and electronic dipole moment cross terms by utilizing mid-IR pump and optical probe fields that are resonant with vibrational and electronic transitions, respectively, in a sample of interest. We detail this newly developed 2D VE spectroscopy experiment and outline the information contained in a 2D VE spectrum. We then use this technique and its single-pump counterpart (1D VE) to probe the vibrational-electronic couplings between high frequency cyanide stretching vibrations (νCN) and either a ligand-to-metal charge transfer transition ([Fe(III)(CN)6](3-) dissolved in formamide) or a metal-to-metal charge transfer (MMCT) transition ([(CN)5Fe(II)CNRu(III)(NH3)5](-) dissolved in formamide). The 2D VE spectra of both molecules reveal peaks resulting from coupled high- and low-frequency vibrational modes to the charge transfer transition. The time-evolving amplitudes and positions of the peaks in the 2D VE spectra report on coherent and incoherent vibrational energy transfer dynamics among the coupled vibrational modes and the charge transfer transition. The selectivity of 2D VE spectroscopy to vibronic processes is evidenced from the selective coupling of specific νCN modes to the MMCT transition in the mixed valence complex. The lineshapes in 2D VE spectra report on the correlation of the frequency fluctuations between the coupled vibrational and electronic frequencies in the mixed valence complex which has a time scale of 1 ps. The details and results of this study confirm the versatility of 2D VE spectroscopy and its applicability to probe how vibrations modulate charge and energy transfer in a wide range of complex molecular, material, and biological systems.

Investigating vibrational relaxation in cyanide-bridged transition metal mixed-valence complexes using two-dimensional infrared and infrared pump-probe spectroscopies

Using polarization-selective two-dimensional infrared (2D IR) and infrared pump-probe spectroscopies, we study vibrational relaxation of the four cyanide stretching (νCN) vibrations found in [(NH3)5Ru(III)NCFe(II)(CN)5](-) (FeRu) dissolved in D2O or formamide and [(NC)5Fe(II)CNPt(IV)(NH3)4NCFe(II)(CN)5](4-) (FePtFe) dissolved in D2O. These cyanide-bridged transition metal complexes serve as models for understanding the role high frequency vibrational modes play in metal-to-metal charge transfers over a bridging ligand. However, there is currently little information about vibrational relaxation and dephasing dynamics of the anharmonically coupled νCN modes in the electronic ground state of these complexes. IR pump-probe experiments reveal that the vibrational lifetimes of the νCN modes are ∼2 times faster when FeRu is dissolved in D2O versus formamide. They also reveal that the vibrational lifetimes of the νCN modes of FePtFe in D2O are almost four times as long as for FeRu in D2O. Combined with mode-specific relaxation dynamics measured from the 2D IR experiments, the IR pump-probe experiments also reveal that intramolecular vibrational relaxation is occurring in all three systems on ∼1 ps timescale. Center line slope dynamics, which have been shown to be a measure of the frequency-frequency correlation function, reveal that the radial, axial, and trans νCN modes exhibit a ∼3 ps timescale for frequency fluctuations. This timescale is attributed to the forming and breaking of hydrogen bonds between each mode and the solvent. The results presented here along with our previous work on FeRu and FePtFe reveal a picture of coupled anharmonic νCN modes where the spectral diffusion and vibrational relaxation dynamics depend on the spatial localization of the mode on the molecular complex and its specific interaction with the solvent.

Two-Dimensional Infrared Study of (13)C-Natural Abundant Vibrational Transition Reveals Intramolecular Vibrational Redistribution Rather than Fluxional Exchange in Mn(CO)5Br

In this work, molecular-symmetry enhanced (13)CO natural abundant isotopic infrared transition was identified in Mn(CO)5Br dissolved in CCl4 by FTIR spectroscopy. Diagonal and associated off-diagonal two-dimensional IR (2D IR) peaks of the (13)CO-species were found to be spectrally separated from the all-(12)CO species, allowing a direct probe of the (13)C natural abundant ensemble. Temperature-dependent FTIR experiment showed no evidence of ligand exchange in the metal carbonyl complex. Intramolecular vibrational redistribution dynamics among the CO stretching vibrational states were extracted using population-time dependent 2D IR diagonal and off-diagonal peaks for both radial mono-(13)CO and all-(12)CO isotopomers. This work demonstrates the potential use of natural abundant isotopic molecular species as a probe for revealing equilibrium and nonequilibrium structural dynamics in condensed-phase molecular systems.

Applications of two-dimensional infrared spectroscopy

Two-dimensional infrared (2D IR) spectroscopy has recently emerged as a powerful tool with applications in many areas of scientific research. The inherent high time resolution coupled with bond-specific spatial resolution of IR spectroscopy enable direct characterization of rapidly interconverting species and fast processes, even in complex systems found in chemistry and biology. In this minireview, we briefly outline the fundamental principles and experimental procedures of 2D IR spectroscopy. Using illustrative example studies, we explain the important features of 2D IR spectra and their capability to elucidate molecular structure and dynamics. Primarily, this minireview aims to convey the scope and potential of 2D IR spectroscopy by highlighting select examples of recent applications including the use of innate or introduced vibrational probes for the study of nucleic acids, peptides/proteins, and materials.

Ultrafast vibrational and structural dynamics of dimeric cyclopentadienyliron dicarbonyl examined by infrared spectroscopy

Equilibrium and ultrafast structural dynamics of a classic transition metal carbonyl compound were revealed by linear and nonlinear infrared methods.

Deactivation mechanism of a novel AIE-active naphthalimide derivative in more polar solutions

We used experimental and computational methods to unveil the deactivation mechanism of a new AIE-active naphthalimide derivative in polar solutions.

The effect on structural and solvent water molecules of substrate binding to ferric horseradish peroxidase

Ultrafast, multi-dimensional infrared spectroscopy, in the form of 2D-IR and pump–probe measurements, has been employed to investigate the effect of substrate binding on the structural dynamics of the horseradish peroxidase (HRP) enzyme. Using nitric oxide bound to the ferric haem of HRP as a sensitive probe of local dynamics, we report measurements of the frequency fluctuations (spectral diffusion) and vibrational lifetime of the NO stretching mode with benzohydroxamic acid (BHA) located in the substrate-binding position at the periphery of the haem pocket, in both D₂O and H₂O solvents. The results reveal that, with BHA bound to the enzyme, the local structural dynamics are insensitive to H/D exchange. These results are in stark contrast to those found in studies of the substrate-free enzyme, which demonstrated that the local chemical and dynamic environment of the haem ligand is influenced by water molecules. In light of the large changes in solvent accessibility caused by substrate binding, we discuss the potential for varying roles for the solvent in the haem pocket of HRP at different stages along the reaction coordinate of the enzymatic mechanism.

Ultrafast infrared spectroscopy reveals water-mediated coherent dynamics in an enzyme active site

Ultrafast infrared spectroscopy provides insights into the dynamic nature of water in the active sites of catalase and peroxidase enzymes.

Understanding the impact of fast dynamics upon the chemical processes occurring within the active sites of proteins and enzymes is a key challenge that continues to attract significant interest, though direct experimental insight in the solution phase remains sparse. Similar gaps in our knowledge exist in understanding the role played by water, either as a solvent or as a structural/dynamic component of the active site. In order to investigate further the potential biological roles of water, we have employed ultrafast multidimensional infrared spectroscopy experiments that directly probe the structural and vibrational dynamics of NO bound to the ferric haem of the catalase enzyme from Corynebacterium glutamicum in both H₂O and D₂O. Despite catalases having what is believed to be a solvent-inaccessible active site, an isotopic dependence of the spectral diffusion and vibrational lifetime parameters of the NO stretching vibration are observed, indicating that water molecules interact directly with the haem ligand. Furthermore, IR pump–probe data feature oscillations originating from the preparation of a coherent superposition of low-frequency vibrational modes in the active site of catalase that are coupled to the haem ligand stretching vibration. Comparisons with an exemplar of the closely-related peroxidase enzyme family shows that they too exhibit solvent-dependent active-site dynamics, supporting the presence of interactions between the haem ligand and water molecules in the active sites of both catalases and peroxidases that may be linked to proton transfer events leading to the formation of the ferryl intermediate Compound I. In addition, a strong, water-mediated, hydrogen bonding structure is suggested to occur in catalase that is not replicated in peroxidase; an observation that may shed light on the origins of the different functions of the two enzymes.

Photoisomerization and structural dynamics of two nitrosylruthenium complexes: a joint study by NMR and nonlinear IR spectroscopies

In this work, the photoisomerization and structural dynamics of two isomeric nitrosylruthenium(ii) complexes [Ru(OAc)(2cqn)2NO] (H2cqn = 2-chloro-8-quinolinol) in CDCl3 and DMSO are examined using NMR and IR spectroscopic methods. The two N atoms in the 2cqn ligand are in trans position in the synthesized cis-1 isomer, while they are in cis position in the cis-2 isomer. Kinetics monitored by NMR spectroscopy shows that the rate constant of photoisomerization from cis-2 to cis-1 isomer depends on the wavelength of irradiation and solvent polarity; it proceeds faster on irradiating near the absorption peak in the UV-Vis region, and also in more polar solvents (DMSO). Density functional theory computation indicates that the peculiarity of [Ru(ii)-NO(+)] group affects the structure and reactivity of the nitrosylruthenium complexes. Using the nitrosyl stretching (νNO) to be vibrational probe, the structural dynamics and structural distributions of the cis-1 and cis-2 isomers are examined by steady-state linear infrared and ultrafast two-dimensional infrared (2D IR) spectroscopies. The structural and photochemical aspects of the observed spectroscopic parameters are discussed in terms of solute-solvent interactions for the two nitrosylruthenium complexes.

Vibrational nano-spectroscopic imaging correlating structure with intermolecular coupling and dynamics

Molecular self-assembly, the function of biomembranes and the performance of organic solar cells rely on nanoscale molecular interactions. Understanding and control of such materials have been impeded by difficulties in imaging their properties with the desired nanometre spatial resolution, attomolar sensitivity and intermolecular spectroscopic specificity. Here we implement vibrational scattering-scanning near-field optical microscopy with high spectral precision to investigate the structure–function relationship in nano-phase separated block copolymers. A vibrational resonance is used as a sensitive reporter of the local chemical environment and we image, with few nanometre spatial resolution and 0.2 cm⁻¹ spectral precision, solvatochromic Stark shifts and line broadening correlated with molecular-scale morphologies. We discriminate local variations in electric fields between nano-domains with quantitative agreement with dielectric continuum models. This ability to directly resolve nanoscale morphology and associated intermolecular interactions can form a basis for the systematic control of functionality in multicomponent soft matter systems.

Quantifying intermolecular coupling and local morphology is important to understand soft matter systems. Pollard et al. show how multispectral vibrational near-field optical microscopy can be used to image molecular-scale morphology and intermolecular interactions with nanometre spatial resolution.

Linear absorption and two-dimensional infrared spectra of N-methylacetamide in chloroform revisited: polarizability and multipole effects

The effect of solvent polarizability and multipole effects on the amide I vibrational spectra of a peptide unit is investigated. Four molecular dynamics force fields of increasing complexity for the solvent are used to model both the linear absorption and two-dimensional infrared spectra. It is observed that, at least in chloroform solution, the predicted solvent shift is considerably improved when accounting for the polarizabiltiy and multipole effects. The latter are typically connected with halogen bonding. Significant deviations are still observed for more sensitive line shape parameters such as the spectral width and line skewness. However, the findings demonstrate that previously observed deviations have an origin in the force field treatment rather than in the electrostatic mapping procedure frequently employed to simulate linear absorption and two-dimensional infrared spectroscopy.

Calculations of the electric fields in liquid solutions

The electric field created by a condensed-phase environment is a powerful and convenient descriptor for intermolecular interactions. Not only does it provide a unifying language to compare many different types of interactions, but it also possesses clear connections to experimental observables, such as vibrational Stark effects. We calculate here the electric fields experienced by a vibrational chromophore (the carbonyl group of acetophenone) in an array of solvents of diverse polarities using molecular dynamics simulations with the AMOEBA polarizable force field. The mean and variance of the calculated electric fields correlate well with solvent-induced frequency shifts and band broadening, suggesting Stark effects as the underlying mechanism of these key solution-phase spectral effects. Compared to fixed-charge and continuum models, AMOEBA was the only model examined that could describe nonpolar, polar, and hydrogen bonding environments in a consistent fashion. Nevertheless, we found that fixed-charge force fields and continuum models were able to replicate some results of the polarizable simulations accurately, allowing us to clearly identify which properties and situations require explicit polarization and/or atomistic representations to be modeled properly, and to identify for which properties and situations simpler models are sufficient. We also discuss the ramifications of these results for modeling electrostatics in complex environments, such as proteins.

pH-dependent picosecond structural dynamics in the distal pocket of nitrophorin 4 investigated by 2D IR spectroscopy

Nitrophorin 4 (NP4) belongs to a family of pH-sensitive, nitric oxide (NO) transporter proteins that undergo a large structural change from a closed to an open conformation at high pH to allow for NO delivery. Measuring the pH-dependent structural dynamics in NP4-NO around the ligand binding site is crucial for developing a mechanistic understanding of NO binding and release. In this study, we use coherent two-dimensional infrared (2D IR) spectroscopy to measure picosecond structural dynamics sampled by the nitrosyl stretch in NP4-NO as a function of pH at room temperature. Our results show that both the closed and open conformers of the protein are present at low (pD 5.1) and high (pD 7.9) pH conditions. The closed and open conformers are characterized by two frequencies of the nitrosyl stretching vibration labeled A0 and A1, respectively. Analysis of the 2D IR line shapes reveals that at pD 5.1, the closed conformer experiences structural fluctuations arising from solvation dynamics on a ∼3 ps time scale. At pD 7.9, both the open and closed conformers exhibit fluctuations on a ∼1 ps time scale. At both pD conditions, the closed conformers maintain a static distribution of structures within the experimental time window of 100 ps. This is in contrast to the open conformer, which is able to interconvert among its substates on a ∼100 ps time scale. Our results directly measure the time scales of solvation dynamics in the distal pocket, the flexibility of the open conformation at high pH, and the rigidity of the closed conformers at both pH conditions. We discuss how the pH-dependent equilibrium structural fluctuations of the nitrosyl ligand measured in this study are related to the uptake and delivery of nitric oxide in NP4.