Vibrational Echo Studies of Protein Dynamics

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

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

Echo dephasing and heat capacity from constrained and unconstrained dynamics of triiodothyronine nuclear receptor protein

The objective of this study is to observe the echo feature curves, vibrational dephasing, and heat capacity of a protein-hormone system taking thyroid hormone receptor-beta (THR-β) as an example. Constrained and unconstrained molecular dynamics simulations are performed by implementing the theory of velocity reassignments to probe the phase coherent state in terms of echo pulses. The constrained vibrations are incorporated by adjusting rigid bonds to all hydrogen atoms with an integrator parameter of 2 fs/step in order to reduce the degrees of freedom whereas 1 fs/step is used in the free vibrations of the atomic cluster. The nature of temperature auto-correlation functions changes so that echo feature curves also show a distinct nature in the cases of constrained and unconstrained vibrations. There is a large variation in kinetic temperature and internal potential energy in the echo time zone. The temperature rate of change of internal potential energy is the main contributor to the heat capacity of the native state protein-hormone system. The heat capacity of proteins estimated from this technique is in good agreement with the values from experiments. This study shows that triiodothyronine (T3) hormone makes some differences in heat capacity upon binding to the THR-β ligand binding domain (LBD). The physical properties of unliganded THR-β and T3-bound THR-β LBD in the cases of constrained and unconstrained dynamics are observed distinctly under the effect of anharmonicity on the phase coherent state of normal modes and the dephasing time lies in a range of 0.6-0.8 ps when the systems are perturbed suddenly.

Ribonuclease S dynamics measured using a nitrile label with 2D IR vibrational echo spectroscopy

A nitrile-labeled amino acid, p-cyanophenylalanine, is introduced near the active site of the semisynthetic enzyme ribonuclease S to serve as a probe of protein dynamics and fluctuations. Ribonuclease S is the limited proteolysis product of subtilisin acting on ribonuclease A, and consists of a small fragment including amino acids 1-20, the S-peptide, and a larger fragment including residues 21-124, the S-protein. A series of two-dimensional vibrational echo experiments performed on the nitrile-labeled S-peptide and the RNase S are described. The time-dependent changes in the two-dimensional infrared vibrational echo line shapes are analyzed using the center line slope method to obtain the frequency-frequency correlation function (FFCF). The observations show that the nitrile probe in the S-peptide has dynamics that are similar to, but faster than, those of the single amino acid p-cyanophenylalanine in water. In contrast, the dynamics of the nitrile label when the peptide is bound to form ribonuclease S are dominated by homogeneous dephasing (motionally narrowed) contributions with only a small contribution from very fast inhomogeneous structural dynamics. The results provide insights into the nature of the structural dynamics of the ribonuclease S complex. The equilibrium dynamics of the nitrile labeled S-peptide and the ribonuclease S complex are also investigated by molecular dynamics simulations. The experimentally determined FFCFs are compared to the FFCFs obtained from the molecular dynamics simulations, thereby testing the capacity of simulations to determine the amplitudes and time scales of protein structural fluctuations on fast time scales under thermal equilibrium conditions.

Dynamics of the folded and unfolded villin headpiece (HP35) measured with ultrafast 2D IR vibrational echo spectroscopy

A series of two-dimensional infrared vibrational echo experiments performed on nitrile-labeled villin headpiece [HP35-(CN)(2)] is described. HP35 is a small peptide composed of three alpha helices in the folded configuration. The dynamics of the folded HP35-(CN)(2) are compared to that of the guanidine-induced unfolded peptide, as well as the nitrile-functionalized phenylalanine (PheCN), which is used to differentiate the peptide dynamic contributions to the observables from those of the water solvent. Because the viscosity of solvent has a significant effect on fast dynamics, the viscosity of the solvent is held constant by adding glycerol. For the folded peptide, the addition of glycerol to the water solvent causes observable slowing of the peptide's dynamics. Holding the viscosity constant as GuHCl is added, the dynamics of unfolded peptide are much faster than those of the folded peptide, and they are very similar to that of PheCN. These observations indicate that the local environment of the nitrile in the unfolded peptide resembles that of PheCN, and the dynamics probed by the CN are dominated by the fluctuations of the solvent molecules, in contrast to the observations on the folded peptide.

Dynamics of a myoglobin mutant enzyme: 2D IR vibrational echo experiments and simulations

Myoglobin (Mb) double mutant T67R/S92D displays peroxidase enzymatic activity in contrast to the wild type protein. The CO adduct of T67R/S92D shows two CO absorption bands corresponding to the A(1) and A(3) substates. The equilibrium protein dynamics for the two distinct substates of the Mb double mutant are investigated by using two-dimensional infrared (2D IR) vibrational echo spectroscopy and molecular dynamics (MD) simulations. The time-dependent changes in the 2D IR vibrational echo line shapes for both of the substates are analyzed using the center line slope (CLS) method to obtain the frequency-frequency correlation function (FFCF). The results for the double mutant are compared to those from the wild type Mb. The experimentally determined FFCF is compared to the FFCF obtained from molecular dynamics simulations, thereby testing the capacity of a force field to determine the amplitudes and time scales of protein structural fluctuations on fast time scales. The results provide insights into the nature of the energy landscape around the free energy minimum of the folded protein structure.

Communication: Decoherence in a nonequilibrium environment: an analytically solvable model

We describe an analytically solvable model of quantum decoherence in a nonequilibrium environment. The model considers the effect of a bath driven from equilibrium by, for example, an ultrafast excitation of a quantum chromophore. The nonequilibrium response of the environment is represented by a nonstationary random function corresponding to the fluctuating transition frequency between two quantum states coupled to the surroundings. The nonstationary random function is characterized by a Fourier series with the phase of each term starting initially with a definite value across the ensemble but undergoing random diffusion with time. The decay of the off-diagonal density matrix element is shown to depend significantly on the particular pattern of initial phases of the terms in the Fourier series, or equivalently, the initial phases of bath modes coupled to the quantum subsystem. This suggests the possibility of control of quantum decoherence by the detailed properties of an environment that is driven from thermal equilibrium.

Residue-specific structural kinetics of proteins through the union of isotope labeling, mid-IR pulse shaping, and coherent 2D IR spectroscopy

We describe a methodology for studying protein kinetics using a rapid-scan technology for collecting 2D IR spectra. In conjunction with isotope labeling, 2D IR spectroscopy is able to probe the secondary structure and environment of individual residues in polypeptides and proteins. It is particularly useful for membrane and aggregate proteins. Our rapid-scan technology relies on a mid-IR pulse shaper that computer generates the pulse shapes, much like in an NMR spectrometer. With this device, data collection is faster, easier, and more accurate. We describe our 2D IR spectrometer, as well as protocols for (13)C(18)O isotope labeling, and then illustrate the technique with an application to the aggregation of the human islet amyloid polypeptide implicated in type 2 diabetes.

Dynamical transition in a small helical peptide and its implication for vibrational energy transport

The two-dimensional infrared spectrum of an octameric helical peptide in chloroform was measured as a function of temperature. Isotope labeling of the carbonyl group of one of the amino acids was used to obtain information for an isolated vibration. The antidiagonal width of the 2D-IR signal, which is a measure of the homogeneous dephasing time T(2), is constant from 220 to 260 K (within experimental error), and increases steeply above. The homogeneous dephasing time of the carbonyl vibration is attributed to the flexibility of the system and/or its immediate surrounding. The system undergoes a dynamical transition at about 270 K, with similarities to the protein dynamical transition. Furthermore, the temperature dependence of the antidiagonal width strongly resembles that of the efficiency of vibrational energy transport along the helix, which has been studied in a recent paper (J. Phys. Chem. B 2008, 112, 15487). The connection between the two processes, structural flexibility and energy transport mechanism, is discussed.

Examination of enzymatic H-tunneling through kinetics and dynamics

In recent years, kinetic measurements of isotope effects of enzyme-catalyzed reactions and their temperature dependence led to the development of theoretical models that were used to rationalize the findings. These models suggested that motions at the femto- to picosecond (fs to ps) time scale modulate the environment of the catalyzed reaction. Due to the fast nature of motions that directly affect the cleavage of a covalent bond, it is challenging to correlate the enzyme kinetics and dynamics related to that step. We report a study of formate dehydrogenase (FDH) that compares the temperature dependence of intrinsic kinetic isotope effects (KIEs) to measurements of the environmental dynamics at the fs-ps time scale (Bandaria et al. J. Am. Chem. Soc. 2008, 130, 22-23). The findings from this comparison of experimental kinetics and dynamics are consistent with models of environmentally coupled H-tunneling models, also known as Marcus-like models. Apparently, at tunneling ready conformations, the donor-acceptor distance, orientation, and fluctuations seems to be well tuned for H-transfer and are not affected by thermal fluctuations slower than 10 ps. This phenomenon has been suggested in the past to be quite general in enzymatic reactions. Here, the kinetics and the dynamics measurements on a single chemical step and on fs-ps time scale, respectively, provide new insight and support for the relevant theoretical models. Furthermore, this methodology could be applied to other systems and be used to examine mutants for which the organization of the donor and acceptor is not ideal, or enzymes with different rigidity and different temperature optimum.

Viscosity-dependent protein dynamics

Spectrally resolved stimulated vibrational echo spectroscopy is used to investigate the dependence of fast protein dynamics on bulk solution viscosity at room temperature in four heme proteins: hemoglobin, myoglobin, a myoglobin mutant with the distal histidine replaced by a valine (H64V), and a cytochrome c552 mutant with the distal methionine replaced by an alanine (M61A). Fructose is added to increase the viscosity of the aqueous protein solutions over many orders of magnitude. The fast dynamics of the four globular proteins were found to be sensitive to solution viscosity and asymptotically approached the dynamical behavior that was previously observed in room temperature sugar glasses. The viscosity-dependent protein dynamics are analyzed in the context of a viscoelastic relaxation model that treats the protein as a deformable breathing sphere. The viscoelastic model is in qualitative agreement with the experimental data but does not capture sufficient system detail to offer a quantitative description of the underlying fluctuation amplitudes and relaxation rates. A calibration method based on the near-infrared spectrum of water overtones was constructed to accurately determine the viscosity of small volumes of protein solutions.

Probing dynamics of complex molecular systems with ultrafast 2D IR vibrational echo spectroscopy

Ultrafast 2D IR vibrational echo spectroscopy is described and a number of experimental examples are given. Details of the experimental method including the pulse sequence, heterodyne detection, and determination of the absorptive component of the 2D spectrum are outlined. As an initial example, the 2D spectrum of the stretching mode of CO bound to the protein myoglobin (MbCO) is presented. The time dependence of the 2D spectrum of MbCO, which is caused by protein structural evolution, is presented and its relationship to the frequency-frequency correlation function is described and used to make protein structural assignments based on comparisons to molecular dynamics simulations. The 2D vibrational echo experiments on the protein horseradish peroxidase are presented. The time dependence of the 2D spectra of the enzyme in the free form and with a substrate bound at the active site are compared and used to examine the influence of substrate binding on the protein's structural dynamics. The application of 2D vibrational echo spectroscopy to the study of chemical exchange under thermal equilibrium conditions is described. 2D vibrational echo chemical exchange spectroscopy is applied to the study of formation and dissociation of organic solute-solvent complexes and to the isomerization around a carbon-carbon single bond of an ethane derivative.

Cytochrome c552 Mutants: Structure and Dynamics at the Active Site Probed by Multidimensional NMR and Vibration Echo Spectroscopy

Spectrally resolved infrared stimulated vibrational echo experiments are used to measure the vibrational dephasing of a CO ligand bound to the heme cofactor in two mutated forms of the cytochrome c552 from Hydrogenobacter thermophilus. The first mutant (Ht-M61A) is characterized by a single mutation of Met61 to an Ala (Ht-M61A), while the second variant is doubly modified to have Gln64 replaced by an Asn in addition to the M61A mutation (Ht-M61A/Q64N). Multidimensional NMR experiments determined that the geometry of residue 64 in the two mutants is consistent with a non-hydrogen-bonding and hydrogen-bonding interaction with the CO ligand for Ht-M61A and Ht-M61A/Q64N, respectively. The vibrational echo experiments reveal that the shortest time scale vibrational dephasing of the CO is faster in the Ht-M61A/Q64N mutant than that in Ht-M61A. Longer time scale dynamics, measured as spectral diffusion, are unchanged by the Q64N modification. Frequency-frequency correlation functions (FFCFs) of the CO are extracted from the vibrational echo data to confirm that the dynamical difference induced by the Q64N mutation is primarily an increase in the fast (hundreds of femtoseconds) frequency fluctuations, while the slower (tens of picoseconds) dynamics are nearly unaffected. We conclude that the faster dynamics in Ht-M61A/Q64N are due to the location of Asn64, which is a hydrogen bond donor, above the heme-bound CO. A similar difference in CO ligand dynamics has been observed in the comparison of the CO derivative of myoglobin (MbCO) and its H64V variant, which is caused by the difference in axial residue interactions with the CO ligand. The results suggest a general trend for rapid ligand vibrational dynamics in the presence of a hydrogen bond donor.

Dynamics of proteins encapsulated in silica sol-gel glasses studied with IR vibrational echo spectroscopy

Spectrally resolved infrared stimulated vibrational echo spectroscopy is used to measure the fast dynamics of heme-bound CO in carbonmonoxy-myoglobin (MbCO) and -hemoglobin (HbCO) embedded in silica sol-gel glasses. On the time scale of approximately 100 fs to several picoseconds, the vibrational dephasing of the heme-bound CO is measurably slower for both MbCO and HbCO relative to that of aqueous protein solutions. The fast structural dynamics of MbCO, as sensed by the heme-bound CO, are influenced more by the sol-gel environment than those of HbCO. Longer time scale structural dynamics (tens of picoseconds), as measured by the extent of spectral diffusion, are the same for both proteins encapsulated in sol-gel glasses compared to that in aqueous solutions. A comparison of the sol-gel experimental results to viscosity-dependent vibrational echo data taken on various mixtures of water and fructose shows that the sol-gel-encapsulated MbCO exhibits dynamics that are the equivalent of the protein in a solution that is nearly 20 times more viscous than bulk water. In contrast, the HbCO dephasing in the sol-gel reflects only a 2-fold increase in viscosity. Attempts to alter the encapsulating pore size by varying the molar ratio of silane precursor to water (R value) used to prepare the sol-gel glasses were found to have no effect on the fast or steady-state spectroscopic results. The vibrational echo data are discussed in the context of solvent confinement and protein-pore wall interactions to provide insights into the influence of a confined environment on the fast structural dynamics experienced by a biomolecule.

FT-IR spectroscopic study on the variations of molecular structures of some carboxyl acids induced by free electron laser

Free electron laser has been developed as tunable lasers over a wide range of wavelengths. Devices irradiating in the region of 6-16 microm (1666-625 cm-1) are operable in the Beijing free electron lasers facilities (BFEL). For understanding the interactions between FEL and biological tissues, in this study wavelength-selective infrared-induced structure changes of substances under irradiation by FEL were measured using FT-IR spectroscopy. The carboxyl acids and carboxylates samples investigated include salicylic acid, sulfosalicylic acid, cholic acid, deoxycholic acid, sodium cholate and sodium deoxycholate. The changes of the FT-IR spectra of the molecules prove that the spectral variations of the samples induced by FEL are closely related to their hydrogen bond networks.

Probing electric fields in protein cavities by using the vibrational stark effect of carbon monoxide

To determine the magnitude and direction of the internal electric field in the Xe4 cavity of myoglobin mutant L29W-S108L, we have studied the vibrational Stark effect of carbon monoxide (CO) using infrared spectroscopy at cryogenic temperatures. CO was photodissociated from the heme iron and deposited selectively in Xe4. Its infrared spectrum exhibits Stark splitting into two bands associated with CO in opposite orientations. Two different photoproduct states can be distinguished, C' and C'', with markedly different properties. For C', characteristic temperature-dependent changes of the area, shift, and width were analyzed, based on a dynamic model in which the CO performs fast librations within a double-well model potential. For the barrier between the wells, a height of approximately 1.8 kJ/mol was obtained, in which the CO performs oscillations at an angular frequency of approximately 25 cm(-1). The magnitude of the electric field in the C' conformation was determined as 11.1 MV/cm; it is tilted by an angle of 29 degrees to the symmetry axis of the potential. Above 140 K, a protein relaxation leads to a significantly altered photoproduct, C'', with a smaller Stark splitting and a more confining potential (barrier >4 kJ/mol) governing the CO librations.

Flexibility and molecular recognition in the immune system

Photon echo spectroscopy has been used to measure the response of three antibody-binding sites to perturbation from electronic excitation of a bound antigen, fluorescein. The three antibodies show motions that range in time scale from tens of femtoseconds to nanoseconds. Relative to the others, one antibody, 4-4-20, possesses a rigid binding site that likely results from a short and inflexible heavy chain complementarity-determining region 3 (HCDR3) loop and a critical Tyr that acts as a "molecular splint," rigidifying the antigen across its most flexible internal degree of freedom. The remaining two antibodies, 34F10 and 40G4, despite being generated against the same antigen, possess binding sites that are considerably more flexible. The more flexible combining sites likely result from longer HCDR3 loops and a deletion in the light chain complementarity-determining region 1 (LCDR1) that removes the critical Tyr residue. The binding site flexibilities may result in varying mechanisms of antigen recognition including lock-and-key, induced-fit, and conformational selection.

Higher-order optical correlation spectroscopy in liquids

Linear optical spectroscopies have long been used to study the behavior of liquids. Laser technology has progressed to the point that it has become possible to perform nonlinear optical experiments that probe higher-order correlation functions in liquids, opening a new window into our understanding of the microscopic details of solution-phase processes. Here we review advances that have been made in recent years in employing higher-order electronic and vibrational spectroscopies to study liquid-state dynamics and structure.

Disorder and funneling effects on exciton migration in treelike dendrimers

The center-bound excitonic diffusion on dendrimers subjected to several types of nonhomogeneous funneling potentials is considered. We first study the mean first passage time (MFPT) for diffusion in a linear potential with different types of correlated and uncorrelated random perturbations. Increasing the funneling force, there is a transition from a phase in which the MFPT grows exponentially with the number of generations g to one in which it does so linearly. Overall the disorder slows down the diffusion, but the effect is much more pronounced in the exponential compared to the linear phase. When the disorder gives rise to uncorrelated random forces there is, in addition, a transition as the temperature T is lowered. This is a transition from a high-T regime in which all paths contribute to the MFPT to a low-T regime in which only a few of them do. We further explore the funneling within a realistic nonlinear potential for extended dendrimers in which the dependence of the lowest excitonic energy level on the segment length was derived using the time-dependent Hatree-Fock approximation. Under this potential the MFPT grows initially linearly with g but crosses over, beyond a molecular-specific and T-dependent optimal size, to an exponential increase. Finally we consider geometrical disorder in the form of a small concentration of long connections as in the small world model. Beyond a critical concentration of connections the MFPT decreases significantly and it changes to a power law or to a logarithmic scaling with g, depending on the strength of the funneling force.

Conformational dynamics of the transcriptional regulator CooA protein studied by subpicosecond mid-infrared vibrational spectroscopy

CooA, which is a transcriptional regulator heme protein allosterically triggered by CO, is studied by femtosecond visible-pump mid-IR-probe spectroscopy. Transient bleaching upon excitation of the heme in the Soret band is detected at approximately 1979 cm(-1), which is the absorption region of the CO bound to the heme. The bleach signal shows a nonexponential decay with time constants of 56 and 290 ps, caused by the rebinding of the CO to the heme. About 98% of dissociated CO recombines geminately. The geminate recombination rate in CooA is significantly faster than those in myoglobin and hemoglobin. The angle of the bound CO with respect to the porphyrin plane is calculated to be about 78 degrees on the basis of the anisotropy measurements. A shift of the bleached mid-IR spectrum of the bound CO is detected and has a characteristic time of 160 ps. It is suggested that the spectral shift is caused by a difference in the frequency of the bound CO in different protein conformations, particularly in an active conformation and in an intermediate one, which is on the way toward an inactive conformation. Thus, the biologically relevant conformation change in CooA was traced. Possible assignment of the observed conformation change is discussed.

Fast protein dynamics probed with infrared vibrational echo experiments

IR vibrational echo experiments are used to study dynamics in myoglobin (Mb) by investigating the dephasing of the CO-stretching mode of CO bound at the active site of the protein (Mb-CO). The temperature dependence and the viscosity dependence of Mb-CO pure dephasing have been measured in several solvents. In low-temperature, glassy solvents, the pure dephasing has a power law temperature dependence, T(1.3), that reflects glasslike protein dynamics. In liquids, the temperature dependence is much steeper and arises from a combination of pure temperature dependence and the influence of decreasing solvent viscosity with increasing temperature. As the solvent viscosity decreases, the ability of the protein's surface to undergo topological fluctuations increases, which in turn increases the internal protein-structural fluctuations. The protein-structural motions are coupled to the CO bound at the active site by electric field fluctuations that accompany movements of polar residues. The dynamic electric field-coupling mechanism is tested by observing differences in the temperature dependence of the pure dephasing of Mb-CO mutations.

Multidimensional femtosecond correlation spectroscopies of electronic and vibrational excitations

Femtosecond visible and infrared analogues of multiple-pulse nuclear magnetic resonance techniques provide novel snapshot probes into the structure and electronic and vibrational dynamics of complex molecular assemblies such as photosynthetic antennae, proteins, and hydrogen-bonded liquids. A classical-oscillator description of these spectroscopies in terms of interacting quasiparticles (rather than transitions among global eigenstates) is developed and sets the stage for designing new pulse sequences and inverting the multidimensional signals to yield molecular structures. Considerable computational advantages and a clear physical insight into the origin of the response and the relevant coherence sizes are provided by a real-space analysis of the underlying coherence-transfer pathways in Liouville space.

Excitonic funneling in extended dendrimers with nonlinear and random potentials

The mean first passage time (MFPT) for photoexcitations diffusion in a funneling potential of artificial treelike light-harvesting antennas (phenylacetylene dendrimers with generation-dependent segment lengths) is computed. Effects of the nonlinearity of the realistic funneling potential and slow random solvent fluctuations considerably slow down the center-bound diffusion beyond a temperature-dependent optimal size. Diffusion on a disordered Cayley tree with a linear potential is investigated analytically. At low temperatures we predict a phase in which the MFPT is dominated by a few paths.

Universal low-frequency vibrations of proteins from a simple interaction potential

A pairwise Born potential connecting the heavy atom sites within a prescribed cutoff, and the equation of motion method (EOM), reproduce the existence of a universal singularity in the low-frequency vibrational density of states of typical globular proteins. This is due to quasilocalization of acoustic waves and an analogy with a similar feature found in glasses is stressed. We explain the dependence of this anomaly with the effective dimensionality of the protein. The EOM method allows for the study of even the largest proteins with a simple personal computer.