Understanding the Vibrational Structure and Ultrafast Dynamics of the Metal Carbonyl Precatalyst [Mn(ppy)(CO)4]

The solution phase structure, vibrational spectroscopy, and ultrafast relaxation dynamics of the precatalyst species [Mn(ppy)(CO)4] (1) in solution have been investigated using ultrafast two-dimensional infrared (2D-IR) spectroscopy. By comparing 2D-IR data with the results of anharmonic density functional theory (DFT) calculations, we establish an excellent agreement between measured and predicted inter-mode couplings of the carbonyl stretching vibrational modes of 1 that relates to the atomic displacements of axial and equatorial ligands in the modes and the nature of the molecular orbitals involved in M-CO bonding. Measurements of IR pump-probe spectra and 2D-IR spectra as a function of waiting time reveal the presence of ultrafast (few ps) intramolecular vibrational energy redistribution between carbonyl stretching modes prior to vibrational relaxation. The vibrational relaxation times of the CO-stretching modes of 1 are found to be relatively solvent-insensitive, suggestive of limited solvent-solute interactions in the ground electronic state. Overall, these data provide a detailed picture of the complex potential energy surface, bonding and vibrational dynamics of 1, establishing a fundamental basis for the next steps in understanding and modulating precatalyst behavior.

Understanding the [NiFe] Hydrogenase Active Site Environment through Ultrafast Infrared and 2D-IR Spectroscopy of the Subsite Analogue K[CpFe(CO)(CN)2] in Polar and Protic Solvents

The [CpFe(CO)(CN)2]- unit is an excellent structural model for the Fe(CO)(CN)2 moiety of the active site found in [NiFe] hydrogenases. Ultrafast infrared (IR) pump-probe and 2D-IR spectroscopy have been used to study K[CpFe(CO)(CN)2] (M1) in a range of protic and polar solvents and as a dry film. Measurements of anharmonicity, intermode vibrational coupling strength, vibrational relaxation time, and solvation dynamics of the CO and CN stretching modes of M1 in H2O, D2O, methanol, dimethyl sulfoxide, and acetonitrile reveal that H-bonding to the CN ligands plays an important role in defining the spectroscopic characteristics and relaxation dynamics of the Fe(CO)(CN)2 unit. Comparisons of the spectroscopic and dynamic data obtained for M1 in solution and in a dry film with those obtained for the enzyme led to the conclusion that the protein backbone forms an important part of the bimetallic active site environment via secondary coordination sphere interactions.

Why are Vibrational Lines Narrow in Proteins?

The vibrational Stark effect in proteins yields line shifts indicative of strong internal electric fields up to a few volts per angstrom. These values are supported by numerical simulations of proteins. The simulations also show a significant breadth of field fluctuations translating to inhomogeneous broadening of vibrational lines. According to fluctuation-dissipation arguments, strong internal fields should lead to broad lines. Experimentally reported vibrational lines in proteins are, however, very narrow. This disconnect is explained here in terms of the insufficient (nonergodic) sampling of the protein's configurations on the lifetime of the vibrational probe. The slow component of the electric field fluctuations in proteins relaxes on the time scale of tens of nanoseconds and is dynamically frozen on the vibrational lifetime.

Understanding the structure and dynamics of hydrogenases by ultrafast and two-dimensional infrared spectroscopy

Hydrogenases are valuable model enzymes for sustainable energy conversion approaches using H2, but rational utilization of these base-metal biocatalysts requires a detailed understanding of the structure and dynamics of their complex active sites. The intrinsic CO and CN- ligands of these metalloenzymes represent ideal chromophores for infrared (IR) spectroscopy, but structural and dynamic insight from conventional IR absorption experiments is limited. Here, we apply ultrafast and two-dimensional (2D) IR spectroscopic techniques, for the first time, to study hydrogenases in detail. Using an O2-tolerant [NiFe] hydrogenase as a model system, we demonstrate that IR pump-probe spectroscopy can explore catalytically relevant ligand bonding by accessing high-lying vibrational states. This ultrafast technique also shows that the protein matrix is influential in vibrational relaxation, which may be relevant for energy dissipation from the active site during fast reaction steps. Further insights into the relevance of the active site environment are provided by 2D-IR spectroscopy, which reveals equilibrium dynamics and structural constraints imposed on the H2-accepting intermediate of [NiFe] hydrogenases. Both techniques offer new strategies for uniquely identifying redox-structural states in complex catalytic mixtures via vibrational quantum beats and 2D-IR off-diagonal peaks. Together, these findings considerably expand the scope of IR spectroscopy in hydrogenase research, and new perspectives for the characterization of these enzymes and other (bio-)organometallic targets are presented.

Solvent Quality Controls Macromolecular Structural Dynamics of a Dendrimeric Hydrogenase Model

We report a spectroscopic investigation of the ultrafast dynamics of the second-generation poly(aryl ether) dendritic hydrogenase model using two-dimensional infrared (2D-IR) spectroscopy to probe the metal carbonyl vibrations of the dendrimer and a reference small molecule, [Fe(μ-S)(CO)₃]₂. We find that the structural dynamics of the dendrimer are reflected in a slow phase of the spectral diffusion, which is absent from [Fe(μ-S)(CO)₃]_2, and we relate the slow phase to the quality of the solvent for poly(aryl ether) dendrimers. We observe a solvent-dependent modulation of the initial phase of vibrational relaxation of the carbonyl groups, which we attribute to an inhibition of solvent assistance in the intramolecular vibrational redistribution process for the dendrimer. There is also a clear solvent dependence of the vibrational frequencies of both the dendrimer and [Fe(μ-S)(CO)₃]₂. Our data represent the first 2D-IR study of a dendritic complex and provide insight into the solvent dependence of molecular conformation in solution and the ultrafast dynamics of moderately sized, conformationally mobile compounds.

Ultrafast Photodynamics of Cyano-Functionalized [FeFe] Hydrogenase Model Compounds

[FeFe] hydrogenases are efficient enzymes that produce hydrogen gas under mild conditions. Synthetic model compounds containing all CO or mixed CO/PMe3 ligands were previously studied by us and others with ultrafast ultraviolet or visible pump-infrared probe spectroscopy in an effort to better understand the function and interactions of the active site with light. Studies of anionic species containing cyano groups, which more closely match the biological active site, have been elusive. In this work, two model compounds dissolved in room-temperature acetonitrile solution were examined: [Fe2(μ-S2C3H6)(CO)4(CN)2]2- (1) and [Fe2(μ-S2C2H4)(CO)4(CN)2]2- (2). These species exhibit long-lived transient signals consistent with loss of one CO ligand with potential isomerization of newly formed ground electronic state photoproducts, as previously observed with all-CO and CO/PMe3-containing models. We find no evidence for fast (ca. 150 ps) relaxation seen in the all-CO and CO/PMe3 compounds because of the absence of the metal-to-metal charge transfer band in the cyano-functionalized models. These results indicate that incorporation of cyano ligands may significantly alter the electronic properties and photoproducts produced immediately after photoexcitation, which may influence the catalytic activity of model compounds when attached to photosensitizers.

Electron-Hole-Pair-Induced Vibrational Energy Relaxation of Rhenium Catalysts on Gold Surfaces

A combination of time-resolved vibrational spectroscopy and density functional theory techniques have been applied to study the vibrational energy relaxation dynamics of the Re(4,4'-dicyano-2,2'-bipyridine)(CO)₃Cl (Re(CO)₃Cl) catalyst for CO₂ to CO conversion bound to gold surfaces. The kinetics of vibrational relaxation exhibits a biexponential decay including an ultrafast initial relaxation and complete recovery of the ground vibrational state. Ab initio molecular dynamics simulations and time-dependent perturbation theory reveal the former to be due to vibrational population exchange between CO stretching modes and the latter to be a combination of intramolecular vibrational relaxation (IVR) and electron-hole pair (EHP)-induced energy transfer into the gold substrate. EHP-induced energy transfer from the Re(CO)₃Cl adsorbate into the gold surface occurs on the same time scale as IVR of Re(CO)₃Cl in aprotic solvents. Therefore, it is expected to be particularly relevant to understanding the reduced catalytic activity of the homogeneous catalyst when anchored to a metal surface.

Dynamic Flexibility of Hydrogenase Active Site Models Studied with 2D-IR Spectroscopy

Hydrogenase enzymes enable organisms to use H₂ as an energy source, having evolved extremely efficient biological catalysts for the reversible oxidation of molecular hydrogen. Small-molecule mimics of these enzymes provide both simplified models of the catalysis reactions and potential artificial catalysts that might be used to facilitate a hydrogen economy. We have studied two diiron hydrogenase mimics, μ-pdt-[Fe(CO)₃]₂ and μ-edt-[Fe(CO)₃]₂ (pdt = propanedithiolate, edt = ethanedithiolate), in a series of alkane solvents and have observed significant ultrafast spectral dynamics using two-dimensional infrared (2D-IR) spectroscopy. Since solvent fluctuations in nonpolar alkanes do not lead to substantial electrostatic modulations in a solute's vibrational mode frequencies, we attribute the spectral diffusion dynamics to intramolecular flexibility. The intramolecular origin is supported by the absence of any measurable solvent viscosity dependence, indicating that the frequency fluctuations are not coupled to the solvent motional dynamics. Quantum chemical calculations reveal a pronounced coupling between the low-frequency torsional rotation of the carbonyl ligands and the terminal CO stretching vibrations. The flexibility of the CO ligands has been proposed to play a central role in the catalytic reaction mechanism, and our results highlight that the CO ligands are highly flexible on a picosecond time scale.

Encapsulating Subsite Analogues of the [FeFe]-Hydrogenases in Micelles Enables Direct Water Interactions

Encapsulation of subsite analogues of the [FeFe]-hydrogenase enzymes in supramolecular structures has been shown to dramatically increase their catalytic ability, but the molecular basis for this enhancement remains unclear. We report the results of experiments employing infrared absorption, ultrafast infrared pump-probe, and 2D-IR spectroscopy to investigate the molecular environment of Fe2(pdt)(CO)6 (pdt: propanedithiolate) [1] encapsulated in the dispersed alkane phase of a heptane-dodecyltrimethylammonium bromide-water microemulsion. It is demonstrated that 1 is partitioned between two molecular environments, one that closely resembles bulk heptane solution and a second that features direct hydrogen-bonding interactions with water molecules that penetrate the surfactant shell. Our results demonstrate that the extent of water access to the normally water-insoluble subsite analogue 1 can be tuned with micelle size, while IR spectroscopy provides a straightforward tool that can be used to measure and fine-tune the chemical environment of catalyst species in self-assembled structures.

Snapshot of the equilibrium dynamics of a drug bound to HIV-1 reverse transcriptase

The anti-AIDS drug rilpivirine undergoes conformational changes to bind HIV-1 reverse transcriptase (RT), which is an essential enzyme for the replication of HIV. These changes allow it to retain potency against mutations that otherwise would render the enzyme resistant. Here we report that water molecules play an essential role in this binding process. Femtosecond experiments and theory expose the molecular level dynamics of rilpivirine bound to HIV-1 RT. Two nitrile substituents, one on each arm of the drug, are used as vibrational probes of the structural dynamics within the binding pocket. Two-dimensional vibrational echo spectroscopy reveals that one nitrile group is unexpectedly hydrogen-bonded to a mobile water molecule, not identified in previous X-ray structures. Ultrafast nitrile-water dynamics are confirmed by simulations. A higher (1.51 Å) resolution X-ray structure also reveals a water-drug interaction network. Maintenance of a crucial anchoring hydrogen bond may help retain the potency of rilpivirine against pocket mutations despite the structural variations they cause.

Solvent-hindered intramolecular vibrational redistribution

Ultrafast two-dimensional infrared spectroscopy and molecular dynamics simulations of Mn(2)(CO)(10) in a series of linear alcohols reveal that the rate of intramolecular vibrational redistribution among the terminal carbonyl stretches is dictated by the average number of hydrogen bonds formed between the solute and solvent. The presence of hydrogen bonds was found to hinder vibrational redistribution between eigenstates, while leaving the overall T(1) relaxation rate unchanged.

Solvent-dependent spectral diffusion in a hydrogen bonded "vibrational aggregate"

Two-dimensional infrared spectroscopy (2DIR) is used to measure the viscosity-dependent spectral diffusion of a model vibrational probe, Mn(2)(CO)(10) (dimanganese decacarbonyl, DMDC), in a series of alcohols with time scales ranging from 2.67 ps in methanol to 5.33 ps in 1-hexanol. Alcohol-alkane solvent mixtures were found to produce indistinguishable linear IR spectra, while still demonstrating viscosity-dependent spectral diffusion. Using a vibrational exciton model to characterize the inhomogeneous energy landscape, several analogies emerge with multichromophoric electronic systems, such as J-aggregates and light-harvesting protein complexes. An excitonic, local vibrational mode Hamiltonian parametrized to reproduce the vibrational structure of DMDC serves as a starting point from which site energies (i.e., local carbonyl frequencies) are given Gaussian distributed disorder. The model gives excellent agreement with both the linear IR spectrum and the inhomogeneous widths extracted from 2DIR, indicating the system can be considered to be a "vibrational aggregate." This model naturally leads to exchange narrowing due to disorder-induced exciton localization, producing line widths consistent with our 1D and 2D measurements. Further, the diagonal disorder alone effectively reduces the molecular symmetry, leading to the appearance of Raman bands in the IR spectrum in accord with the measurements. Here, we show that the static inhomogeneity of the excitonic model with disorder successfully captures the essential details of the 1D spectrum while predicting the degree of IR activity of forbidden modes as well as the inhomogeneous widths and relative magnitudes of the transition moments.

Temperature dependence of vibrational frequency fluctuation of N(3) (-) in D(2)O

We have studied the temperature dependence of the vibrational frequency fluctuation of the antisymmetric stretching mode of N(3) (-) in D(2)O by three-pulse infrared (IR) photon echo experiments. IR pump-probe measurements were also carried out to investigate the population relaxation and the orientational relaxation of the same band. It was found that the time-correlation function (TCF) of the frequency fluctuation of this mode is well described by a biexponential function with a quasistatic term. The faster decay component has a time constant of about 0.1 ps, and the slower component varies from 1.4 to 1.1 ps in the temperature range from 283 to 353 K. This result indicates that liquid dynamics related to the frequency fluctuation are not highly sensitive to temperature. We discuss the relationship between the temperature dependence of the vibrational frequency fluctuation and that of the molecular motion of the system to investigate the molecular origin of the frequency fluctuation of the solute. We compare the temperature dependence of the frequency fluctuation with that of other dynamics such as dielectric relaxation of water. In contrast to the Debye dielectric relaxation time of D(2)O, the two time constants of the TCF of the frequency fluctuation do not exhibit strong temperature dependence. We propose a simple theoretical model for the frequency fluctuation in solutions based on perturbation theory and the dipole-dipole interaction between the vibrational mode of the solute and the solvent molecules. This model suggests that the neighboring solvent molecules in the vicinity of the solute play an important role in the frequency fluctuation. We suggest that the picosecond component of the frequency fluctuation results from structural fluctuation of the hydrogen-bonding network in water.