Observation of coherent reaction dynamics in heme proteins

Kinetic Basis for the Design of Azobenzene-Based Photoswitchable A2a Adenosine Receptor Ligands

Photoisomerization of ligands is a key process in the field of photopharmacology. Thus, the kinetics and efficiency of this initial photoreaction are of great importance but can be influenced by the molecular environment of the binding pocket and the resulting confinement of the reaction pathway. In this study, we investigated the photoisomerization of an azobenzene derivative of the anti-Parkinson's drug istradefylline. To identify the impact of the binding pocket, the ligand was examined in solution and bound to its target protein, the A2a adenosine receptor (A2aR), belonging to the family of G protein-coupled receptors (GPCRs). Although the overall efficiency of isomerization is reduced when the ligand is bound, the initial photoreaction experiences little influence from the binding pocket. However, protein-coupled motion promotes a longer-lived excited-state population and thus leads to a reduction in efficiency. The results provide the kinetic basis for a photoswitchable GPCR ligand and demonstrate the influence of the binding pocket on fundamental photochemistry.

Vibrational-Mode-Selective Modulation of Electronic Excitation

Vibrational-mode-selective modulation of electronic excitation is conducted with a synchronized femtosecond (fs) visible (vis) pulse and a picosecond (ps) infrared (IR) pulse. The mechanism of modulation of vibrational and vibronic relaxation behavior of excited state is investigated with ultrafast vis/IR, IR/IR, and vis-IR/IR transient spectroscopy, optical gating experiments and theoretical calculations. An organic molecule, 4'-(N,N-dimethylamino)-3-methoxyflavone (DMA3MHF) is chosen as the model system. Upon 1608 cm-1 excitation, the skeleton stretching vibration of DMA3MHF is energized, which can significantly change the shape of the absorption, facilitate the radiative decay and promote emission from vibrational excited states. As results, a remarkable enhancement and a slight blueshift in fluorescence are observed. The mode-selective modulation of electronic excitation is not limited in luminescence or photophysics. It is expected to be widely applicable in tuning many photochemical processes.

Watching atomically resolved structural dynamics-Blinded by the light

In a recent issue of Nature, Barends et al.1 studied the photodissociation of carboxymyoglobin with ultrafast laser pump-probe serial femtosecond crystallography experiments. They observed significant differences in heme protein structural dynamics for biologically relevant 1-photon excitation relative to high excitation leading to the absorption of several photons per heme.

Ultrafast proton transfer of the aqueous phenol radical cation

Proton transfer (PT) reactions are fundamental to numerous chemical and biological processes. While sub-picosecond PT involving electronically excited states has been extensively studied, little is known about ultrafast PT triggered by photoionization. Here, we employ femtosecond optical pump-probe spectroscopy and quantum dynamics calculations to investigate the ultrafast proton transfer dynamics of the aqueous phenol radical cation (PhOH˙⁺). Analysis of the vibrational wave packet dynamics reveals unusually short dephasing times of 0.18 ± 0.02 ps and 0.16 ± 0.02 ps for the PhOH˙⁺ O-H wag and bend frequencies, respectively, suggestive of ultrafast PT occurring on the ∼0.1 ps timescale. The reduced potential energy surface obtained from ab initio calculations shows that PT is barrierless when it is coupled to the intermolecular hindered translation between PhOH˙⁺ and the proton-acceptor water molecule. Quantum dynamics calculations yield a lifetime of 193 fs for PhOH˙⁺, in good agreement with the experimental results and consistent with the PT reaction being mediated by the intermolecular O⋯O stretch. These results suggest that photoionization can be harnessed to produce photoacids that undergo ultrafast PT. In addition, they also show that PT can serve as an ultrafast deactivation channel for limiting the oxidative damage potential of radical cations.

Spectroscopic observation and ultrafast coherent vibrational dynamics of the aqueous phenylalanine radical

The phenylalanine radical (Phe˙) has been proposed to mediate biological electron transport (ET) and exhibit long-lived electronic coherences following attosecond photoionization. However, the coupling of ultrafast structural reorganization to the oxidation/ionization of biomolecules such as phenylalanine remains unexplored. Moreover, studies of ET involving Phe˙ are hindered by its hitherto unobserved electronic spectrum. Here, we report the spectroscopic observation and coherent vibrational dynamics of aqueous Phe˙, prepared by sub-6 fs photodetachment of phenylalaninate anions. Sub-picosecond transient absorption spectroscopy reveals the ultraviolet absorption signature of Phe˙. Ultrafast structural reorganization drives coherent vibrational motion involving nine fundamental frequencies and one overtone. DFT calculations rationalize the absence of the decarboxylation reaction, a photodegradation pathway previously identified for Phe˙. Our findings guide the interpretation of future attosecond experiments aimed at elucidating coherent electron motion in photoionized aqueous biomolecules and pave way for the spectroscopic identification of Phe˙ in studies of biological ET.

Mapping Structural Dynamics of Proteins with Femtosecond Stimulated Raman Spectroscopy

The structure-function relationships of biomolecules have captured the interest and imagination of the scientific community and general public since the field of structural biology emerged to enable the molecular understanding of life processes. Proteins that play numerous functional roles in cellular processes have remained in the forefront of research, inspiring new characterization techniques. In this review, we present key theoretical concepts and recent experimental strategies using femtosecond stimulated Raman spectroscopy (FSRS) to map the structural dynamics of proteins, highlighting the flexible chromophores on ultrafast timescales. In particular, wavelength-tunable FSRS exploits dynamic resonance conditions to track transient-species-dependent vibrational motions, enabling rational design to alter functions. Various ways of capturing excited-state chromophore structural snapshots in the time and/or frequency domains are discussed. Continuous development of experimental methodologies, synergistic correlation with theoretical modeling, and the expansion to other nonequilibrium, photoswitchable, and controllable protein systems will greatly advance the chemical, physical, and biological sciences.

Unveiling coupled electronic and vibrational motions of chromophores in condensed phases

The quest for capturing molecular movies of functional systems has motivated scientists and engineers for decades. A fundamental understanding of electronic and nuclear motions, two principal components of the molecular Schrödinger equation, has the potential to enable the de novo rational design for targeted functionalities of molecular machines. We discuss the development and application of a relatively new structural dynamics technique, femtosecond stimulated Raman spectroscopy with broadly tunable laser pulses from the UV to near-IR region, in tracking the coupled electronic and vibrational motions of organic chromophores in solution and protein environments. Such light-sensitive moieties hold broad interest and significance in gaining fundamental knowledge about the intramolecular and intermolecular Hamiltonian and developing effective strategies to control macroscopic properties. Inspired by recent experimental and theoretical advances, we focus on the in situ characterization and spectroscopy-guided tuning of photoacidity, excited state proton transfer pathways, emission color, and internal conversion via a conical intersection.

A molecular movie of ultrafast singlet fission

The complex dynamics of ultrafast photoinduced reactions are governed by their evolution along vibronically coupled potential energy surfaces. It is now often possible to identify such processes, but a detailed depiction of the crucial nuclear degrees of freedom involved typically remains elusive. Here, combining excited-state time-domain Raman spectroscopy and tree-tensor network state simulations, we construct the full 108-atom molecular movie of ultrafast singlet fission in a pentacene dimer, explicitly treating 252 vibrational modes on 5 electronic states. We assign the tuning and coupling modes, quantifying their relative intensities and contributions, and demonstrate how these modes coherently synchronise to drive the reaction. Our combined experimental and theoretical approach reveals the atomic-scale singlet fission mechanism and can be generalized to other ultrafast photoinduced reactions in complex systems. This will enable mechanistic insight on a detailed structural level, with the ultimate aim to rationally design molecules to maximise the efficiency of photoinduced reactions.

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

A Nuclear Resonance Vibrational Spectroscopic Study of Oxy Myoglobins Reconstituted with Chemically Modified Heme Cofactors: Insights into the Fe-O2 Bonding and Internal Dynamics of the Protein

The molecular mechanism of O₂ binding to hemoglobin (Hb) and myoglobin (Mb) is a long-standing issue in the field of bioinorganic and biophysical chemistry. The nature of Fe-O₂ bond in oxy Hb and Mb had been extensively investigated by resonance Raman spectroscopy, which assigned the Fe-O₂ stretching bands at ∼570 cm^-1. However, resonance Raman assignment of the vibrational mode had been elusive due to the spectroscopic selection rule and to the limited information available about the ground-state molecular structure. Thus, nuclear resonance vibrational spectroscopy was applied to oxy Mbs reconstituted with ⁵⁷Fe-labeled native heme cofactor and two chemically modified ones. This advanced spectroscopy in conjunction with DFT analyses gave new insights into the nature of the Fe-O₂ bond of oxy heme by revealing the effect of heme peripheral substitutions on the vibrational dynamics of heme Fe atom, where the main Fe-O₂ stretching band of the native protein was characterized at ∼420 cm^-1.

Spiers Memorial Lecture. Introductory lecture: the impact of structure on photoinduced processes in nucleic acids and proteins

Light is an important environmental variable and most organisms have evolved means to sense, exploit or avoid it and to repair detrimental effects on their genome. In general, light absorption is the task of specific chromophores, however other biomolecules such as oligonucleotides also do so which can result in undesired outcomes such as mutations and cancer. Given the biological importance of light-induced processes and applications for imaging, optogenetics, photodynamic therapy or photovoltaics, there is a great interest in understanding the detailed molecular mechanisms of photoinduced processes in proteins and nucleic acids. The processes are typically characterized by time-resolved spectroscopic approaches or computation, inferring structural information on transient species from stable ground state structures. Recently, however, structure determination of excited states or other short-lived species has become possible with the advent of X-ray free-electron lasers. This review gives an overview of the impact of structure on the understanding of photoinduced processes in macromolecules, focusing on systems presented at this Faraday Discussion meeting.

What Can Be Learned from Nuclear Resonance Vibrational Spectroscopy: Vibrational Dynamics and Hemes

Nuclear resonance vibrational spectroscopy (NRVS; also known as nuclear inelastic scattering, NIS) is a synchrotron-based method that reveals the full spectrum of vibrational dynamics for Mössbauer nuclei. Another major advantage, in addition to its completeness (no arbitrary optical selection rules), is the unique selectivity of NRVS. The basics of this recently developed technique are first introduced with descriptions of the experimental requirements and data analysis including the details of mode assignments. We discuss the use of NRVS to probe ⁵⁷Fe at the center of heme and heme protein derivatives yielding the vibrational density of states for the iron. The application to derivatives with diatomic ligands (O₂, NO, CO, CN^–) shows the strong capabilities of identifying mode character. The availability of the complete vibrational spectrum of iron allows the identification of modes not available by other techniques. This permits the correlation of frequency with other physical properties. A significant example is the correlation we find between the Fe–Im stretch in six-coordinate Fe(XO) hemes and the trans Fe–N(Im) bond distance, not possible previously. NRVS also provides uniquely quantitative insight into the dynamics of the iron. For example, it provides a model-independent means of characterizing the strength of iron coordination. Prediction of the temperature-dependent mean-squared displacement from NRVS measurements yields a vibrational “baseline” for Fe dynamics that can be compared with results from techniques that probe longer time scales to yield quantitative insights into additional dynamical processes.

Tuning calcium biosensors with a single-site mutation: structural dynamics insights from femtosecond Raman spectroscopy

Structural dynamics governing the emission properties of a single-site mutant of fluorescent-protein-based calcium biosensors are elucidated by femtosecond stimulated Raman spectroscopy.

Kinetic Control of O2 Reactivity in H-NOX Domains

Transient absorption, resonance Raman, and vibrational coherence spectroscopies are used to investigate the mechanisms of NO and O2 binding to WT Tt H-NOX and its P115A mutant. Vibrational coherence spectra of the oxy complexes provide clear evidence for the enhancement of an iron-histidine mode near 217 cm(-1) following photoexcitation, which indicates that O2 can be dissociated in these proteins. However, the quantum yield of O2 photolysis is low, particularly in the wild type (≲3%). Geminate recombination of O2 and NO in both of these proteins is very fast (∼1.4 × 10(11) s(-1)) and highly efficient. We show that the distal heme pocket of the H-NOX system forms an efficient trap that limits the O2 off-rate and determines the overall affinity. The distal pocket hydrogen bond, which appears to be stronger in the P115A mutant, may help retard the O2 ligand from escaping into the solvent following either photoinduced or thermal dissociation. This, along with a strengthening of the Fe-O2 bond that is correlated with the significant heme ruffing and saddling distortions, explains the unusually high O2 affinity of WT Tt H-NOX and the even higher affinity found in the P115A mutant.

Panoramic portrait of primary molecular events preceding excited state proton transfer in water

Primary events that power ultrafast excited state proton transfer in water are revealed to involve coupled intermolecular and intramolecular motions.

Photochemistry powers numerous processes from luminescence and human vision, to light harvesting. However, the elucidation of multidimensional photochemical reaction coordinates on molecular timescales remains challenging. We developed wavelength-tunable femtosecond stimulated Raman spectroscopy to simultaneously achieve pre-resonance enhancement for transient reactant and product species of the widely used photoacid pyranine undergoing excited-state proton transfer (ESPT) reaction in solution. In the low-frequency region, the 280 cm^–1 ring deformation mode following 400 nm photoexcitation exhibits pronounced intensity oscillations on the sub-picosecond timescale due to anharmonic vibrational coupling to the 180 cm^–1 hydrogen-bond stretching mode only in ESPT-capable solvents, indicating a primary event of functional relevance. This leads to the contact ion pair formation on the 3 ps timescale before diffusion-controlled separation. The intermolecular 180 cm^–1 mode also reveals vibrational cooling time constants, ∼500 fs and 45 ps in both H₂O and D₂O, which differ from ESPT time constants of ∼3/8 and 90/250 ps in H₂O/D₂O, respectively. Spectral results using H₂¹⁸O further substantiate the functional role of the intermolecular 180 cm^–1 mode in modulating the distance between proton donor and acceptor and forming the transient ion pair. The direct observation of molecular structural evolution across a wide spectral region during photochemical reactions enriches our fundamental understanding of potential energy surface and holds the key to advancing energy and biological sciences with exceptional atomic and temporal precision.

3D Motions of Iron in Six-Coordinate {FeNO}(7) Hemes by Nuclear Resonance Vibration Spectroscopy

The vibrational spectrum of a six-coordinate nitrosyl iron porphyrinate, monoclinic [Fe(TpFPP)(1-MeIm)(NO)] (TpFPP=tetra-para-fluorophenylporphyrin; 1-MeIm=1-methylimidazole), has been studied by oriented single-crystal nuclear resonance vibrational spectroscopy (NRVS). The crystal was oriented to give spectra perpendicular to the porphyrin plane and two in-plane spectra perpendicular or parallel to the projection of the FeNO plane. These enable assignment of the FeNO bending and stretching modes. The measurements reveal that the two in-plane spectra have substantial differences that result from the strongly bonded axial NO ligand. The direction of the in-plane iron motion is found to be largely parallel and perpendicular to the projection of the bent FeNO on the porphyrin plane. The out-of-plane Fe-N-O stretching and bending modes are strongly mixed with each other, as well as with porphyrin ligand modes. The stretch is mixed with v50 as was also observed for dioxygen complexes. The frequency of the assigned stretching mode of eight Fe-X-O (X=N, C, and O) complexes is correlated with the Fe-XO bond lengths. The nature of highest frequency band at ≈560 cm(-1) has also been examined in two additional new derivatives. Previously assigned as the Fe-NO stretch (by resonance Raman), it is better described as the bend, as the motion of the central nitrogen atom of the FeNO group is very large. There is significant mixing of this mode. The results emphasize the importance of mode mixing; the extent of mixing must be related to the peripheral phenyl substituents.

Atomistic Simulation of Adiabatic Reactive Processes Based on Multi-State Potential Energy Surfaces

The adiabatic reactive molecular dynamics (ARMD) method provides a framework to study chemical reactions using molecular dynamics simulations with minimal computational overhead. Here, ARMD is generalized to an arbitrary reactive process between two states in which reactants and products can be treated by an atomistic force field. The implementation is described, and the method is applied to two systems: the kinetics of NO rebinding to myoglobin (Mb) as a validation system and the conformational transition in neuroglobin (Ngb) which explores the full functionality of ARMD. For MbNO, the nonexponential kinetics observed both in experiment and earlier ARMD studies is reproduced. Furthermore, the sensitivity of the results with respect to the asymptotic separation between the two potential energy surfaces (NO bound and unbound) is studied.

Direct observation of ultrafast collective motions in CO myoglobin upon ligand dissociation

The hemoprotein myoglobin is a model system for the study of protein dynamics. We used time-resolved serial femtosecond crystallography at an x-ray free-electron laser to resolve the ultrafast structural changes in the carbonmonoxy myoglobin complex upon photolysis of the Fe-CO bond. Structural changes appear throughout the protein within 500 femtoseconds, with the C, F, and H helices moving away from the heme cofactor and the E and A helices moving toward it. These collective movements are predicted by hybrid quantum mechanics/molecular mechanics simulations. Together with the observed oscillations of residues contacting the heme, our calculations support the prediction that an immediate collective response of the protein occurs upon ligand dissociation, as a result of heme vibrational modes coupling to global modes of the protein.

Coherent Nuclear Wave Packets in Q States by Ultrafast Internal Conversions in Free Base Tetraphenylporphyrin

Persistence of vibrational coherence in electronic transition has been noted especially in biochemical systems. Here, we report the dynamics between electronic excited states in free base tetraphenylporphyrin (H2TPP) by time-resolved fluorescence with high time resolution. Following the photoexcitation of the B state, ultrafast internal conversion occurs to the Qx state directly as well as via the Qy state. Unique and distinct coherent nuclear wave packet motions in the Qx and Qy states are observed through the modulation of the fluorescence intensity in time. The instant, serial internal conversions from the B to the Qy and Qx states generate the coherent wave packets. Theory and experiment show that the observed vibrational modes involve the out-of-plane vibrations of the porphyrin ring that are strongly coupled to the internal conversion of H2TPP.

Observing heme doming in myoglobin with femtosecond X-ray absorption spectroscopy

We report time-resolved X-ray absorption measurements after photolysis of carbonmonoxy myoglobin performed at the LCLS X-ray free electron laser with nearly 100 fs (FWHM) time resolution. Data at the Fe K-edge reveal that the photoinduced structural changes at the heme occur in two steps, with a faster (∼70 fs) relaxation preceding a slower (∼400 fs) one. We tentatively attribute the first relaxation to a structural rearrangement induced by photolysis involving essentially only the heme chromophore and the second relaxation to a residual Fe motion out of the heme plane that is coupled to the displacement of myoglobin F-helix.

Comprehensive Fe-ligand vibration identification in {FeNO}6 hemes

Oriented single-crystal nuclear resonance vibrational spectroscopy (NRVS) has been used to obtain all iron vibrations in two {FeNO}(6) porphyrinate complexes, five-coordinate [Fe(OEP)(NO)]ClO4 and six-coordinate [Fe(OEP)(2-MeHIm)(NO)]ClO4. A new crystal structure was required for measurements of [Fe(OEP)(2-MeHIm)(NO)]ClO4, and the new structure is reported herein. Single crystals of both complexes were oriented to be either parallel or perpendicular to the porphyrin plane and/or axial imidazole ligand plane. Thus, the FeNO bending and stretching modes can now be unambiguously assigned; the pattern of shifts in frequency as a function of coordination number can also be determined. The pattern is quite distinct from those found for CO or {FeNO}(7) heme species. This is the result of unchanging Fe-N(NO) bonding interactions in the {FeNO}(6) species, in distinct contrast to the other diatomic ligand species. DFT calculations were also used to obtain detailed predictions of vibrational modes. Predictions were consistent with the intensity and character found in the experimental spectra. The NRVS data allow the assignment and observation of the challenging to obtain Fe-Im stretch in six-coordinate heme derivatives. NRVS data for this and related six-coordinate hemes with the diatomic ligands CO, NO, and O2 reveal a strong correlation between the Fe-Im stretch and Fe-N(Im) bond distance that is detailed for the first time.

Population-controlled impulsive vibrational spectroscopy: background- and baseline-free Raman spectroscopy of excited electronic states

We have developed the technique of population-controlled impulsive vibrational spectroscopy (PC-IVS) aimed at providing high-quality, background-free Raman spectra of excited electronic states and their dynamics. Our approach consists of a modified transient absorption experiment using an ultrashort (<10 fs) pump pulse with additional electronic excitation and control pulses. The latter allows for the experimental isolation of excited-state vibrational coherence and, hence, vibrational spectra. We illustrate the capabilities of PC-IVS by reporting the Raman spectra of well-established molecular systems such as the carotenoid astaxanthin and trans-stilbene and present the first excited-state Raman spectra of the retinal protonated Schiff base chromophore in solution. Our approach, illustrated here with impulsive vibrational spectroscopy, is equally applicable to transient and even multidimensional infrared and electronic spectroscopies to experimentally isolate spectroscopic signatures of interest.

Excited state structural events of a dual-emission fluorescent protein biosensor for Ca²⁺ imaging studied by femtosecond stimulated Raman spectroscopy

Fluorescent proteins (FPs) are luminescent biomolecules that emit characteristic hues upon irradiation. A group of calmodulin (CaM)-green FP (GFP) chimeras have been previously engineered to enable the optical detection of calcium ions (Ca(2+)). We investigate one of these genetically encoded Ca(2+) biosensors for optical imaging (GECOs), GEM-GECO1, which fluoresces green without Ca(2+) but blue with Ca(2+), using femtosecond stimulated Raman spectroscopy (FSRS). The time-resolved FSRS data (<800 cm(-1)) reveal that initial structural evolution following 400 nm photoexcitation involves small-scale coherent proton motions on both ends of the chromophore two-ring system with a <250 fs time constant. Upon Ca(2+) binding, the chromophore adopts a more twisted conformation in the protein pocket with increased hydrophobicity, which inhibits excited-state proton transfer (ESPT) by effectively trapping the protonated chromophore in S1. Both the chromophore photoacidity and local environment form the ultrafast structural dynamics basis for the dual-emission properties of GEM-GECO1. Its photochemical transformations along multidimensional reaction coordinates are evinced by distinct stages of FSRS spectral evolution, particularly related to the ∼460 and 504 cm(-1) modes. The direct observation of lower frequency modes provides crucial information about the nuclear motions preceding ESPT, which enriches our understanding of photochemistry and enables the rational design of new biosensors.

Excited-state structural dynamics of a dual-emission calmodulin-green fluorescent protein sensor for calcium ion imaging

Fluorescent proteins (FPs) have played a pivotal role in bioimaging and advancing biomedicine. The versatile fluorescence from engineered, genetically encodable FP variants greatly enhances cellular imaging capabilities, which are dictated by excited-state structural dynamics of the embedded chromophore inside the protein pocket. Visualization of the molecular choreography of the photoexcited chromophore requires a spectroscopic technique capable of resolving atomic motions on the intrinsic timescale of femtosecond to picosecond. We use femtosecond stimulated Raman spectroscopy to study the excited-state conformational dynamics of a recently developed FP-calmodulin biosensor, GEM-GECO1, for calcium ion (Ca(2+)) sensing. This study reveals that, in the absence of Ca(2+), the dominant skeletal motion is a ∼ 170 cm(-1) phenol-ring in-plane rocking that facilitates excited-state proton transfer (ESPT) with a time constant of ∼ 30 ps (6 times slower than wild-type GFP) to reach the green fluorescent state. The functional relevance of the motion is corroborated by molecular dynamics simulations. Upon Ca(2+) binding, this in-plane rocking motion diminishes, and blue emission from a trapped photoexcited neutral chromophore dominates because ESPT is inhibited. Fluorescence properties of site-specific protein mutants lend further support to functional roles of key residues including proline 377 in modulating the H-bonding network and fluorescence outcome. These crucial structural dynamics insights will aid rational design in bioengineering to generate versatile, robust, and more sensitive optical sensors to detect Ca(2+) in physiologically relevant environments.

Vibrationally coherent crossing and coupling of electronic states during internal conversion in β-carotene

Coupling of nuclear and electronic degrees of freedom mediates energy flow in molecules after optical excitation. The associated coherent dynamics in polyatomic systems, however, remain experimentally unexplored. Here, we combined transient absorption spectroscopy with electronic population control to reveal nuclear wave packet dynamics during the S2 → S1 internal conversion in β-carotene. We show that passage through a conical intersection is vibrationally coherent and thereby provides direct feedback on the role of different vibrational coordinates in the breakdown of the Born-Oppenheimer approximation.

Anisotropic iron motion in nitrosyl iron porphyrinates: natural and synthetic hemes

The vibrational spectra of two five-coordinate nitrosyl iron porphyrinates, [Fe(OEP)(NO)] (OEP = dianion of 2,3,7,8,12,13,17,18-octaethylporphyrin) and [Fe(DPIX)(NO)] (DPIX = deuteroporphyrin IX), have been studied by oriented single-crystal nuclear resonance vibrational spectroscopy. Single crystals (both are in the triclinic crystal system) were oriented to give vibrational spectra perpendicular to the porphyrin plane. Additionally, two orthogonal in-plane measurements that were also either perpendicular or parallel to the projection of the FeNO plane onto the porphyrin plane yield the complete set of vibrations with iron motion. In addition to cleanly enabling the assignment of the FeNO bending and stretching modes, the measurements reveal that the two in-plane spectra from the parallel and perpendicular in-plane directions for both compounds have substantial differences. The assignment of these in-plane vibrations were aided by density functional theory predictions. The differences in the two in-plane directions result from the strongly bonded axial NO ligand. The direction of the in-plane iron motion is thus found to be largely parallel and perpendicular to the projection of the FeNO plane on the porphyrin plane. These axial ligand effects on the in-plane iron motion are related to the strength of the axial ligand-to-iron bond.

Long-range electron transfer with myoglobin immobilized at Au/mixed-SAM junctions: mechanistic impact of the strong protein confinement

Horse muscle myoglobin (Mb) was tightly immobilized at Au-deposited ~15-Å-thick mixed-type (1:1) alkanethiol SAMs, HS-(CH₂)₁₁-COOH/HS-(CH₂)₁₁-OH, and placed in contact with buffered H₂O or D₂O solutions. Fast-scan cyclic voltammetry (CV) and a Marcus-equation-based analysis were applied to determine unimolecular standard rate constants and reorganization free energies for electron transfer (ET), under variable-temperature (15-55 °C) and -pressure (0.01-150 MPa) conditions. The CV signal was surprisingly stable and reproducible even after multiple temperature and pressure cycles. The data analysis revealed the following values: standard rate constant, 33 s⁻¹ (25 °C, 0.01 MPa, H₂O); reorganization free energy, 0.5 ± 0.1 eV (throughout); activation enthalpy, 12 ± 3 kJ mol⁻¹; activation volume, -3.1 ± 0.2 cm³ mol⁻¹; and pH-dependent solvent kinetic isotope effect (k(H)⁰/k(D)⁰), 0.7-1.4. Furthermore, the values for the rate constant and reorganization free energy are very similar to those previously found for cytochrome c electrostatically immobilized at the monocomponent Au/HS-(CH₂)₁₁-COOH junction. In vivo, Mb apparently forms a natural electrostatic complex with cytochrome b₅ (cyt-b₅) through the "dynamic" (loose) docking pattern, allowing for a slow ET that is intrinsically coupled to the water's removal from the "defective" heme iron (altogether shaping the biological repair mechanism for Mb's "met" form). In contrary, our experiments rather mimic the case of a "simple" (tight) docking of the redesigned (mutant) Mb with cyt-b₅ (Nocek et al. J. Am. Chem. Soc. 2010, 132, 6165-6175). According to our analysis, in this configuration, Mb's distal pocket (linked to the "ligand channel") seems to be arrested within the restricted configuration, allowing the rate-determining reversible ET process to be coupled only to the inner-sphere reorganization (minimal elongation/shortening of an Fe-OH₂ bond) rather than the pronounced detachment (rebinding) of water and, hence, to be much faster.

Investigations of the low-frequency spectral density of cytochrome c upon equilibrium unfolding

The equilibrium unfolding process of ferric horse heart cytochrome c (cyt c), induced by guanidinium hydrochloride (GdHCl), was studied using UV-vis absorption spectroscopy, resonance Raman spectroscopy, and vibrational coherence spectroscopy (VCS). The unfolding process was successfully fit using a three-state model which included the fully folded (N) and unfolded (U) states, along with an intermediate (I) assigned to a Lys bound heme. The VCS spectra revealed for the first time several low-frequency heme modes that are sensitive to cyt c unfolding: γ(a) (~50 cm(-1)), γ(b) (~80 cm(-1)), γ(c) (~100 cm(-1)), and ν(s)(His-Fe-His) at 205 cm(-1). These out-of-plane modes have potential functional relevance and are activated by protein-induced heme distortions. The free energies for the N-I and the I-U transitions at pH 7.0 and 20 °C were found to be 4.6 kcal/M and 11.6 kcal/M, respectively. Imidazole was also introduced to replace the methionine ligand so the unfolding can be modeled as a two-state system. The intensity of the mode γ(b)~80 cm(-1) remains nearly constant during the unfolding process, while the amplitudes of the other low frequency modes track with spectral changes observed at higher frequency. This confirms that the heme deformation changes are coupled to the protein tertiary structural changes that take place upon unfolding. These studies also reveal that damping of the coherent oscillations depends sensitively on the coupling between heme and the surrounding water solvent.

Early time excited-state structural evolution of pyranine in methanol revealed by femtosecond stimulated Raman spectroscopy

To understand chemical reactivity of molecules in condensed phase in real time, a structural dynamics technique capable of monitoring molecular conformational motions on their intrinsic time scales, typically on femtoseconds to picoseconds, is needed. We have studied a strong photoacid pyranine (8-hydroxypyrene-1,3,6-trisulfonic acid, HPTS, pK(a)* ≈ 0) in pure methanol and observed that excited-state proton transfer (ESPT) is absent, in sharp contrast with our previous work on HPTS in aqueous solutions wherein ESPT prevails following photoexcitation. Two transient vibrational marker bands at ~1477 (1454) and 1532 (1528) cm(-1) appear in CH3OH (CD3OD), respectively, rising within the instrument response time of ~140 fs and decaying with 390-470 (490-1400) fs and ~200 ps time constants in CH3OH (CD3OD). We attribute the mode onset to small-scale coherent proton motion along the pre-existing H-bonding chain between HPTS and methanol, and the two decay stages to the low-frequency skeletal motion-modulated Franck-Condon relaxation within ~1 ps and subsequent rotational diffusion of H-bonding partners in solution before fluorescence. The early time kinetic isotope effect (KIE) of ~3 upon methanol deuteration argues active proton motions particularly within the first few picoseconds when coherent skeletal motions are underdamped. Pronounced quantum beats are observed for high-frequency modes consisting of strong phenolic COH rocking (1532 cm(-1)) or H-out-of-plane wagging motions (952 cm(-1)) due to anharmonic coupling to coherent low-frequency modes impulsively excited at ca. 96, 120, and 168 cm(-1). The vivid illustration of atomic motions of HPTS in varying H-bonding geometry with neighboring methanol molecules unravels the multidimensional energy relaxation pathways immediately following photoexcitation, and provides compelling evidence that, in lieu of ESPT, the photoacidity of HPTS promptly activates characteristic low-frequency skeletal motions to search phase space mainly concerning the phenolic end and to efficiently dissipate vibrational energy via skeletal deformation and proton shuttling motions within the intermediate, relatively confined excited-state HPTS-methanol complex on a solvent-dependent dynamic potential energy surface.

Fluorescence of tryptophan in designed hairpin and Trp-cage miniproteins: measurements of fluorescence yields and calculations by quantum mechanical molecular dynamics simulations

The quantum yield of tryptophan (Trp) fluorescence was measured in 30 designed miniproteins (17 β-hairpins and 13 Trp-cage peptides), each containing a single Trp residue. Measurements were made in D(2)O and H(2)O to distinguish between fluorescence quenching mechanisms involving electron and proton transfer in the hairpin peptides, and at two temperatures to check for effects of partial unfolding of the Trp-cage peptides. The extent of folding of all the peptides also was measured by NMR. The fluorescence yields ranged from 0.01 in some of the Trp-cage peptides to 0.27 in some hairpins. Fluorescence quenching was found to occur by electron transfer from the excited indole ring of the Trp to a backbone amide group or the protonated side chain of a nearby histidine, glutamate, aspartate, tyrosine, or cysteine residue. Ionized tyrosine side chains quenched strongly by resonance energy transfer or electron transfer to the excited indole ring. Hybrid classical/quantum mechanical molecular dynamics simulations were performed by a method that optimized induced electric dipoles separately for the ground and excited states in multiple π-π* and charge-transfer (CT) excitations. Twenty 0.5 ns trajectories in the tryptophan's lowest excited singlet π-π* state were run for each peptide, beginning by projections from trajectories in the ground state. Fluorescence quenching was correlated with the availability of a CT or exciton state that was strongly coupled to the π-π* state and that matched or fell below the π-π* state in energy. The fluorescence yields predicted by summing the calculated rates of charge and energy transfer are in good accord with the measured yields.

Ultrafast conformational dynamics of pyranine during excited state proton transfer in aqueous solution revealed by femtosecond stimulated Raman spectroscopy

Proton transfer reaction plays an essential role in a myriad of chemical and biological processes, and to reveal the choreography of the proton motion intra- and intermolecularly, a spectroscopic technique capable of capturing molecular structural snapshots on the intrinsic time scale of proton transfer motions is needed. The photoacid pyranine (8-hydroxypyrene-1,3,6-trisulfonic acid, HPTS) serves as a paradigm case to dissect excited state proton transfer (ESPT) events in aqueous solution, triggered precisely by photoexcitation. We have used femtosecond stimulated Raman spectroscopy (FSRS) to yield novel insights into the ultrafast conformational dynamics of photoexcited HPTS in complex with water and acetate molecules. Marker bands attributed to the deprotonated form of HPTS (1139 cm(-1), ∼220 fs rise) appear earlier and faster than the monomer acetic acid peak (864 cm(-1), ∼530 fs rise), indicating that water molecules actively participate in the ESPT chain. Several key low-frequency modes at 106, 150, 195, and 321 cm(-1) have been identified to facilitate ESPT at different stages from 300 fs, 1 ps, to 6 ps and beyond, having distinctive dynamics contributing through hydrogen bonds with 0, 1, and more intervening water molecules. The time-resolved FSRS spectroscopy renders a direct approach to observe the reactive coupling between the vibrational degrees of freedom of photoexcited HPTS in action, therefore revealing the anharmonicity matrix both within HPTS and between HPTS and the neighboring acceptor molecules. The observed excited state conformational dynamics are along the ESPT multidimensional reaction coordinate and are responsible for the photoacidity of HPTS in aqueous solution.

Nuclear resonance vibrational spectra of five-coordinate imidazole-ligated iron(II) porphyrinates

Nuclear resonance vibrational spectra have been obtained for six five-coordinate imidazole-ligated iron(II) porphyrinates, [Fe(Por)(L)] (Por = tetraphenylporphyrinate, octaethylporphyrinate, tetratolylporphyrinate, or protoporphyrinate IX and L = 2-methylimidazole or 1,2-dimethylimidazole). Measurements have been made on both powder and oriented crystal samples. The spectra are dominated by strong signals around 200-300 cm(-1). Although the in-plane and out-of-plane vibrations are seriously overlapped, oriented crystal spectra allow their deconvolution. Thus, oriented crystal experimental data, along with density functional theory (DFT) calculations, enable the assignment of key vibrations in the spectra. Molecular dynamics are also discussed. The nature of the Fe-N(Im) vibrations has been elaborated further than was possible from resonance Raman studies. Our study suggests that the Fe motions are coupled with the porphyrin core and peripheral groups motions. Both peripheral groups and their conformations have significant influence on the vibrational spectra (position and shape).

Spectroscopic identification of reactive porphyrin motions

Nuclear resonance vibrational spectroscopy (NRVS) reveals the vibrational dynamics of a Mössbauer probe nucleus. Here, (57)Fe NRVS measurements yield the complete spectrum of Fe vibrations in halide complexes of iron porphyrins. Iron porphine serves as a useful symmetric model for the more complex spectrum of asymmetric heme molecules that contribute to numerous essential biological processes. Quantitative comparison with the vibrational density of states (VDOS) predicted for the Fe atom by density functional theory calculations unambiguously identifies the correct sextet ground state in each case. These experimentally authenticated calculations then provide detailed normal mode descriptions for each observed vibration. All Fe-ligand vibrations are clearly identified despite the high symmetry of the Fe environment. Low frequency molecular distortions and acoustic lattice modes also contribute to the experimental signal. Correlation matrices compare vibrations between different molecules and yield a detailed picture of how heme vibrations evolve in response to (a) halide binding and (b) asymmetric placement of porphyrin side chains. The side chains strongly influence the energetics of heme doming motions that control Fe reactivity, which are easily observed in the experimental signal.