Signatures of beta-sheet secondary structures in linear and two-dimensional infrared spectroscopy

Using idealized models for parallel and antiparallel beta sheets, we calculate the linear and two-dimensional infrared spectra of the amide I vibration as a function of size and secondary structure. The model assumes transition-dipole coupling between the amide I oscillators in the sheet and accounts for the anharmonic nature of these oscillators. Using analytical and numerical methods, we show that the nature of the one-quantum vibrational eigenstates, which govern the linear spectrum, is, to a large extent, determined by the symmetry of the system and the relative magnitude of interstrand interactions. We also find that the eigenstates, in particular their trends with system size, depend sensitively on the secondary structure of the sheet. While in practice these differences may be difficult to distinguish in congested linear spectra, we demonstrate that they give rise to promising markers for secondary structure in the two-dimensional spectra. In particular, distinct differences occur between the spectra of parallel and antiparallel beta sheets and between beta hairpins and extended beta sheets.     

Signatures of beta-peptide unfolding in two-dimensional vibrational echo spectroscopy: a simulation study.

An ensemble of exciton Hamiltonians for the amide-I band of the folded and unfolded states of a helical beta-heptapeptide is generated using a molecular dynamics (MD) simulation. The correlated fluctuations of its parameters and their signatures in two-dimensional (2D) vibrational echo spectroscopy are computed. This technique uses infrared pulse sequences to provide ultrafast snapshots of molecular structural fluctuations, in analogy with multidimensional NMR. The present study demonstrates that, by combining a method of calculating the vibrational Hamiltonian from MD snapshots and the nonlinear exciton equations (NEE), it may be possible to simulate realistic multidimensional IR spectra of chemically and biologically interesting systems.     

Heterogeneous exciton dynamics revealed by two-dimensional optical spectroscopy

We show that optical two-dimensional (2D) spectroscopy can recover ultrafast heterogeneous dynamics of closely spaced delocalized exciton states from a molecular exciton manifold characterized by a single absorption band. The complete experimental third-order nonlinear optical response from room-temperature J-aggregates in liquid phase is reproduced for the first time with self-consistent Frenkel exciton theory combined with modified Redfield theory. We show that exciton relaxation between the exciton states and nuclear-motion-induced exchange-narrowed energy fluctuations of individual delocalized exciton states can be distinguished because these two processes lead to a distinctively different evolution of the absolute 2D spectrum. Our technique also allows recovery of the variation of the exciton relaxation rates as well as the degree of exciton delocalization across the absorption band.     

Backbone dynamics of proteins studied by two-dimensional heteronuclear NMR spectroscopy and molecular dynamics simulations

To investigate the backbone dynamics of proteins ¹⁵N longitudinal and transverse relaxation experiments combined with {¹H, ¹⁵N{ NOE measurements together with molecular dynamics simulations were carried out using ribonuclease T₁ and the complex of ribonuclease T₁ with 2′GMP as a model protein. The intensity decay of individual amide cross peaks in a series of (¹H, ¹⁵N)HSQC spectra with appropriate relaxation periods was fitted to a single exponential by using a simplex algorithm in order to obtain ¹⁵N T₁ and T₂ relaxation times. The relaxation times were analyzed in terms of the “model-free” approach introduced by Lipari and Szabo. In addition, a nanosecond molecular dynamics (MD ) simulation of ribonuclease T₁ and its 2′GMP complex in water was carried out. The angular reorientations of the backbone amide groups were classified with several coordinate frames following a transformation of NH vector trajectories. In this study, NH librations and backbone dihedral angle fluctuations were distinguished. The NH bond librations were found to be similar for all amides as characterized by correlation times of librational motions in a subpicosecond scale. The angular amplitudes of these motions were found to be about 10°–12° for out-of-plane displacements and 3°–5° for the in-plane displacement. The contributions from the much slower backbone dihedral angle fluctuations strongly depend on the secondary structure. The dependence of the amplitude of local motion on the residue location in the backbone is in good agreement with the results of NMR relaxation measurements and the X-ray data. The protein dynamics is characterized by a highly restricted local motion of those parts of the backbone with defined secondary structure as well as by a high flexibility in loop regions. Comparison of the MD and NMR data of the free liganded enzyme ribonuclease T₁ clearly indicates a restriction of the mobility within certain regions of the backbone upon inhibitor binding. © 1996 John Wiley & Sons, Inc.     

Excited state dynamics in photosynthetic reaction center and light harvesting complex 1

Key to efficient harvesting of sunlight in photosynthesis is the first energy conversion process in which electronic excitation establishes a trans-membrane charge gradient. This conversion is accomplished by the photosynthetic reaction center (RC) that is, in case of the purple photosynthetic bacterium Rhodobacter sphaeroides studied here, surrounded by light harvesting complex 1 (LH1). The RC employs six pigment molecules to initiate the conversion: four bacteriochlorophylls and two bacteriopheophytins. The excited states of these pigments interact very strongly and are simultaneously influenced by the surrounding thermal protein environment. Likewise, LH1 employs 32 bacteriochlorophylls influenced in their excited state dynamics by strong interaction between the pigments and by interaction with the protein environment. Modeling the excited state dynamics in the RC as well as in LH1 requires theoretical methods, which account for both pigment-pigment interaction and pigment-environment interaction. In the present study we describe the excitation dynamics within a RC and excitation transfer between light harvesting complex 1 (LH1) and RC, employing the hierarchical equation of motion method. For this purpose a set of model parameters that reproduce RC as well as LH1 spectra and observed oscillatory excitation dynamics in the RC is suggested. We find that the environment has a significant effect on LH1-RC excitation transfer and that excitation transfers incoherently between LH1 and RC.     

A phenomenological approach to modeling chemical dynamics in nonlinear and two-dimensional spectroscopy

We present an approach for calculating nonlinear spectroscopic observables, which overcomes the approximations inherent to current phenomenological models without requiring the computational cost of performing molecular dynamics simulations. The trajectory mapping method uses the semi-classical approximation to linear and nonlinear response functions, and calculates spectra from trajectories of the system's transition frequencies and transition dipole moments. It rests on identifying dynamical variables important to the problem, treating the dynamics of these variables stochastically, and then generating correlated trajectories of spectroscopic quantities by mapping from the dynamical variables. This approach allows one to describe non-Gaussian dynamics, correlated dynamics between variables of the system, and nonlinear relationships between spectroscopic variables of the system and the bath such as non-Condon effects. We illustrate the approach by applying it to three examples that are often not adequately treated by existing analytical models--the non-Condon effect in the nonlinear infrared spectra of water, non-Gaussian dynamics inherent to strongly hydrogen bonded systems, and chemical exchange processes in barrier crossing reactions. The methods described are generally applicable to nonlinear spectroscopy throughout the optical, infrared and terahertz regions.     

Excited and ground state vibrational dynamics revealed by two-dimensional electronic spectroscopy

Broadband two-dimensional electronic spectroscopy (2DES) can assist in understanding complex electronic and vibrational signatures. In this paper, we use 2DES to examine the electronic structure and dynamics of a long chain cyanine dye (1,1-diethyl-4,4-dicarbocyanine iodide, or DDCI-4), a system with a vibrational progression. Using broadband pulses that span the resonant electronic transition, we measure two-dimensional spectra that show a characteristic six peak pattern from coherently excited ground and excited state vibrational modes. We model these features using a spectral density formalism and the vibronic features are assigned to Feynman pathways. We also examine the dynamics of a particular set of peaks demonstrating anticorrelated peak motion, a signature of oscillatory wavepacket dynamics on the ground and excited states. These dynamics, in concert with the general structure of vibronic two-dimensional spectra, can be used to distinguish between pure electronic and vibrational quantum coherences.     

Reaction dynamics of a molecular switch unveiled by coherent two-dimensional electronic spectroscopy

Coherent two-dimensional electronic spectroscopy is usually employed on molecular species with fixed geometric configuration. Here we present two-dimensional Fourier-transform electronic spectra of dissolved 6,8-dinitro-1',3',3'-trimethylspiro[2H-1-benzopyran-2,2'-indoline] (6,8-dinitro-BIPS), a photochromic system present in two ring-open forms differing in the cis/trans configuration of a double bond, which both undergo a photoinduced ring closure. The two-dimensional spectra, recorded with 20 fs pump pulses centered at 605 nm and a supercontinuum probe covering the complete visible spectral range, allow for a detailed analysis of the photophysics and photochemistry of the two isomers and directly reveal that cis/trans isomerization among them does not play a major role. This experiment demonstrates the potential of two-dimensional electronic spectroscopy for reactive processes.     

Lattice theory of ultrafast excitonic and charge-transfer dynamics in DNA

We propose a lattice fermion model suitable for studying the ultrafast photoexcitation dynamics of ordered chains of deoxyribonucleic acid (DNA) polymers. The model includes both parallel (intrachain) and perpendicular (cross-chain) terms as well as diagonal cross-chain terms coupling neighboring bases. The general form of our Hamiltonian is borrowed from lattice fermion models of quantum chromodynamics. The band structure for this model can be determined analytically, and we use this as a basis for computing the singly excited states of the poly(dA)poly(dT) DNA duplex using configuration interaction singles. Parameters for the model are taken from various literature sources and our own ab initio calculations. Results indicate that the excited states consist of a low energy band of dark charge-separated states followed by separate bands of delocalized excitonic states which have weak mixing between the thymidine and adenosine sides of the DNA chain. We then propose a lattice exciton model based upon the transition dipole-dipole couplings between bases and compare the analytical results for the survival probability of an initially localized exciton to exact numerical results. The results herein underscore the competing role of excitonic and charge-transfer dynamics in these systems.     

Influence of environment induced correlated fluctuations in electronic coupling on coherent excitation energy transfer dynamics in model photosynthetic systems

Two-dimensional photon-echo experiments indicate that excitation energy transfer between chromophores near the reaction center of the photosynthetic purple bacterium Rhodobacter sphaeroides occurs coherently with decoherence times of hundreds of femtoseconds, comparable to the energy transfer time scale in these systems. The original explanation of this observation suggested that correlated fluctuations in chromophore excitation energies, driven by large scale protein motions could result in long lived coherent energy transfer dynamics. However, no significant site energy correlation has been found in recent molecular dynamics simulations of several model light harvesting systems. Instead, there is evidence of correlated fluctuations in site energy-electronic coupling and electronic coupling-electronic coupling. The roles of these different types of correlations in excitation energy transfer dynamics are not yet thoroughly understood, though the effects of site energy correlations have been well studied. In this paper, we introduce several general models that can realistically describe the effects of various types of correlated fluctuations in chromophore properties and systematically study the behavior of these models using general methods for treating dissipative quantum dynamics in complex multi-chromophore systems. The effects of correlation between site energy and inter-site electronic couplings are explored in a two state model of excitation energy transfer between the accessory bacteriochlorophyll and bacteriopheophytin in a reaction center system and we find that these types of correlated fluctuations can enhance or suppress coherence and transfer rate simultaneously. In contrast, models for correlated fluctuations in chromophore excitation energies show enhanced coherent dynamics but necessarily show decrease in excitation energy transfer rate accompanying such coherence enhancement. Finally, for a three state model of the Fenna-Matthews-Olsen light harvesting complex, we explore the influence of including correlations in inter-chromophore couplings between different chromophore dimers that share a common chromophore. We find that the relative sign of the different correlations can have profound influence on decoherence time and energy transfer rate and can provide sensitive control of relaxation in these complex quantum dynamical open systems.     

Multidimensional infrared spectroscopy of water. I. Vibrational dynamics in two-dimensional IR line shapes

In this and the following paper, we describe the ultrafast structural fluctuations and rearrangements of the hydrogen bonding network of water using two-dimensional (2D) infrared spectroscopy. 2D IR spectra covering all the relevant time scales of molecular dynamics of the hydrogen bonding network of water were studied for the OH stretching absorption of HOD in D2O. Time-dependent evolution of the 2D IR line shape serves as a spectroscopic observable that tracks how different hydrogen bonding environments interconvert while changes in spectral intensity result from vibrational relaxation and molecular reorientation of the OH dipole. For waiting times up to the vibrational lifetime of 700 fs, changes in the 2D line shape reflect the spectral evolution of OH oscillators induced by hydrogen bond dynamics. These dynamics, characterized through a set of 2D line shape analysis metrics, show a rapid 60 fs decay, an underdamped oscillation on a 130 fs time scale induced by hydrogen bond stretching, and a long time decay constant of 1.4 ps. 2D surfaces for waiting times larger than 700 fs are dominated by the effects of vibrational relaxation and the thermalization of this excess energy by the solvent bath. Our modeling based on fluctuations with Gaussian statistics is able to reproduce the changes in dispersed pump-probe and 2D IR spectra induced by these relaxation processes, but misses the asymmetry resulting from frequency-dependent spectral diffusion. The dynamical origin of this asymmetry is discussed in the companion paper.     

Ultrafast intermolecular dynamics of liquid water: a theoretical study on two-dimensional infrared spectroscopy

Physical and chemical properties of liquid water are dominated by hydrogen bond structure and dynamics. Recent studies on nonlinear vibrational spectroscopy of intramolecular motion provided new insight into ultrafast hydrogen bond dynamics. However, our understanding of intermolecular dynamics of water is still limited. We theoretically investigated the intermolecular dynamics of liquid water in terms of two-dimensional infrared (2D IR) spectroscopy. The 2D IR spectrum of intermolecular frequency region (<1000 cm(-1)) is calculated by using the equilibrium and nonequilibrium hybrid molecular dynamics method. We find the ultrafast loss of the correlation of the libration motion with the time scale of approximately 110 fs. It is also found that the energy relaxation from the libration motion to the low frequency motion takes place with the time scale of about 180 fs. We analyze the effect of the hindered translation motion on these ultrafast dynamics. It is shown that both the frequency modulation of libration motion and the energy relaxation from the libration to the low frequency motion significantly slow down in the absence of the hindered translation motion. The present result reveals that the anharmonic coupling between the hindered translation and libration motions is essential for the ultrafast relaxation dynamics in liquid water.     

Two-dimensional optical three-pulse photon echo spectroscopy. II. Signatures of coherent electronic motion and exciton population transfer in dimer two-dimensional spectra

Using the nonperturbative approach to the calculation of nonlinear optical spectra developed in a foregoing paper [Mancal et al., J. Chem. Phys. 124, 234504 (2006), preceding paper], calculations of two-dimensional electronic spectra of an excitonically coupled dimer model system are presented. The dissipative exciton transfer dynamics is treated within the Redfield theory and energetic disorder within the molecular ensemble is taken into account. The manner in which the two-dimensional spectra reveal electronic couplings in the aggregate system and the evolution of the spectra in time is studied in detail. Changes in the intensity and shape of the peaks in the two-dimensional relaxation spectra are related to the coherent and dissipative dynamics of the system. It is shown that coherent electronic motion, an electronic analog of a vibrational wave packet, can manifest itself in two-dimensional optical spectra of molecular aggregate systems as a periodic modulation of both the diagonal and off-diagonal peaks.     

Study of water dynamics and distances in paramagnetic solids by variable-temperature two-dimensional 2H NMR spectroscopy

A recently proposed two-dimensional (2)H NMR experiment is used to measure the (2)H (spin I=1) quadrupolar and paramagnetic shift anisotropy interactions in powdered CuCl(2).2D(2)O as a function of temperature. The principal components of the quadrupolar and paramagnetic shift anisotropy tensors and the Euler angles describing the orientations of the tensors in the molecular frame are determined at each temperature. For this purpose an analytical approach is introduced to extract desired parameters from motionally averaged two-dimensional line shapes where the averaging is introduced by rapid 180 degrees flips around C(2) axes of D(2)O molecules. This approach can be readily applied to study various materials containing water of crystallization. It is also clearly shown that the rapid continuous rotation of D(2)O molecules around their C(2) axes is not taking place in the studied solid in the range of temperatures between 209 and 344 K. Once the paramagnetic shift anisotropy of a deuterium atom is measured accurately it is used to estimate the distance between deuterium and the nearest copper atom bearing an unpaired electron. Excellent agreement is found between structural parameters obtained in this study and those provided by neutron and x-ray diffraction, showing that the paramagnetic shift anisotropy is a sensitive probe of distances in paramagnetic solids.     

Ion-pairing dynamics of Li+ and SCN- in dimethylformamide solution: chemical exchange two-dimensional infrared spectroscopy

Ultrafast two-dimensional infrared (2DIR) spectroscopy has been proven to be an exceptionally useful method to study chemical exchange processes between different vibrational chromophores under thermal equilibria. Here, we present experimental results on the thermal equilibrium ion pairing dynamics of Li(+) and SCN(-) ions in N,N-dimethylformamide. Li(+) and SCN(-) ions can form a contact ion pair (CIP). Varying the relative concentration of Li(+) in solution, we could control the equilibrium CIP and free SCN(-) concentrations. Since the CN stretch frequency of Li-SCN CIP is blue-shifted by about 16 cm(-1) from that of free SCN(-) ion, the CN stretch IR spectrum is a doublet. The temperature-dependent IR absorption spectra reveal that the CIP formation is an endothermic (0.57 kJ∕mol) process and the CIP state has larger entropy by 3.12 J∕(K?mol) than the free ion states. Since the two ionic configurations are spectrally distinguishable, this salt solution is ideally suited for nonlinear IR spectroscopic investigations to study ion pair association and dissociation dynamics. Using polarization-controlled IR pump-probe methods, we first measured the lifetimes and orientational relaxation times of these two forms of ionic configurations. The vibrational population relaxation times of both the free ion and CIP are about 32 ps. However, the orientational relaxation time of the CIP, which is ～47 ps, is significantly longer than that of the free SCN(-), which is ～7.7 ps. This clearly indicates that the effective moment of inertia of the CIP is much larger than that of the free SCN(-). Then, using chemical exchange 2DIR spectroscopy and analyzing the diagonal peak and cross-peak amplitude changes with increasing the waiting time, we determined the contact ion pair association and dissociation time constants that are found to be 165 and 190 ps, respectively. The results presented and discussed in this paper are believed to be important, not only because the ion-pairing dynamics is one of the most fundamental physical chemistry problems but also because such molecular ion-ion interactions are of critical importance in understanding Hofmeister effects on protein stability.     

Electronic spectroscopy and excited state dynamics of the Al-H2/D2 complex

The electronic spectra of the Al-H2 and Al-D2 complexes are investigated in a collaborative experimental and theoretical study. The complexes were prepared in a pulsed supersonic beam and detected with laser fluorescence excitation spectroscopy. Transitions to bound vibrational levels in electronic states correlating with the excited-state Al(3d, 4p, 4d) + H2/D2 asymptotes were observed by monitoring emission from lower excited Al atomic levels, formed in the non-radiative decay of the excited complex. Fluorescence depletion has also been used to verify that the observed Al-H2 bands all involve the same molecular carrier. The bands have been assigned to the more strongly bound Al-oH2 and Al-pD2 nuclear spin modifications. In contrast to our previous observations for Al(5s)-H2 [X. Yang and P. J. Dagdigian, J. Chem. Phys., 1998, 109, 8920], for which only one potential energy surface (PES) emanates from the dissociation asymptote, the Lorentzian widths of the different vibrational bands in the 3d, 4p, 4d<--3p transitions vary widely, in some cases allowing resolution of the rotational structure of the bands. With the help of the calculated Al(3p)-oH2/pD2 dissociation energies, binding energies of the observed excited vibronic levels are reported. The mechanism of predissociation is investigated theoretically through ab initio calculation of C2 nu cuts of the excited PESs. It is concluded that predissociation occurs through coupling with the repulsive Al(4s)-H2 PES. With these calculations, a qualitative interpretation of the observed bands could be made.     

Site-specific vibrational dynamics of the CD3zeta membrane peptide using heterodyned two-dimensional infrared photon echo spectroscopy

Heterodyned two-dimensional infrared (2D IR) spectroscopy has been used to study the amide I vibrational dynamics of a 27-residue peptide in lipid vesicles that encompasses the transmembrane domain of the T-cell receptor CD3zeta. Using 1-(13)C[Double Bond](18)O isotope labeling, the amide I mode of the 49-Leucine residue was spectroscopically isolated and the homogeneous and inhomogeneous linewidths of this mode were measured by fitting the 2D IR spectrum collected with a photon echo pulse sequence. The pure dephasing and inhomogeneous linewidths are 2 and 32 cm(-1), respectively. The population relaxation time of the amide I band was measured with a transient grating, and it contributes 9 cm(-1) to the linewidth. Comparison of the 49-Leucine amide I mode and the amide I band of the entire CD3zeta peptide reveals that the vibrational dynamics are not uniform along the length of the peptide. Possible origins for the large amount of inhomogeneity present at the 49-Leucine site are discussed.     

Submolecular dynamics of amorphous polymers probed by two-dimensional infrared (2D IR) spectroscopy

Two-dimensional infrared (2D IR) spectroscopy was used to probe the submolecular dynamics of atactic polystyrene. 2D IR is a powerful analytical technique especially suited for the elucidation of localized motions of polymer segments. In 2D IR, a polymer sample is excited by a small-amplitude oscillatory strain at a frequency in the acoustic range. The fluctuation of IR dichroism signals resulting from the strain-induced reorientation of electric dipole-transition moments is monitored with a time-resolved spectrometer. Spectra defined by two independent wavenumber axes are constructed by applying a correlation analysis to such signals. The 2D spectra provide detailed information about the local dynamics of submolecular constituents of the system. From the sign of cross peaks in the synchronous 2D IR spectrum of glassy polystyrene, it is shown that the main-chain backbone of polystyrene reorients in the direction of applied strain. Cross peaks in an asynchronous 2D IR spectrum reveal highly localized reorientational motions of phenyl side groups occurring more or less independently of the main-chain realignment. In the glassy state, the phenyl ring tends to fold back along the main-chain, indicating that there exists a highly constrained local distortion of side groups during deformation.