Femtosecond studies of tryptophan fluorescence dynamics in proteins: local solvation and electronic quenching

We report our systematic examination of tryptophan fluorescence dynamics in proteins with femtosecond resolution. Distinct patterns of femtosecond-resolved fluorescence transients from the blue to the red side of emission have been characterized to distinguish local ultrafast solvation and electronic quenching. It is shown that tryptophan is an ideal local optical probe for hydration dynamics and protein-water interactions as well as an excellent local molecular reporter for ultrafast electron transfer in proteins, as demonstrated by a series of biological systems, here in melittin, human serum albumin, and human thioredoxin, and at lipid interfaces. These studies clarify the assignments in the literature about the ultrafast solvation or quenching dynamics of tryptophan in proteins. We also report a new observation of solvation dynamics at far red-side emission when the relaxation of the local environment is slower than 1 ps. These results provide a molecular basis for using tryptophan as a local molecular probe for ultrafast protein dynamics in general.     

Ultrafast hydration dynamics in melittin folding and aggregation: helix formation and tetramer self-assembly

Melittin, an amphipathic peptide from honeybee venom, consists of 26 amino acid residues and adopts different conformations from a random coil, to an alpha-helix, and to a self-assembled tetramer under certain aqueous environments. We report here our systematic studies of the hydration dynamics in these conformations using single intrinsic tryptophan (W19) as a molecular probe. With femtosecond resolution, we observed the solvation dynamics occurring in 0.62 and 14.7 ps in a random-coiled primary structure. The former represents bulklike water motion, and the latter reflects surface-type hydration dynamics of proteins. As a comparison, a model tripeptide (KWK) was also studied. At a membrane-water interface, melittin folds into a secondary alpha-helical structure, and the interfacial water motion was found to take as long as 114 ps, indicating a well-ordered water structure along the membrane surface. In high-salt aqueous solution, the dielectric screening and ionic solvation promote the hydrophobic core collapse in melittin aggregation and facilitate the tetramer formation. This self-assembled tertiary structure is also stabilized by the strong hydrophilic interactions of charged C-terminal residues and associated ions with water molecules in the two assembled regions. The hydration dynamics was observed to occur in 87 ps, significantly slower than typical water relaxation at protein surfaces but similar to water motion at membrane interfaces. Thus, the observed time scale of approximately 100 ps probably implies appropriate water mobility for mediating the formation of high-order structures of melittin in an alpha-helix and a self-assembled tetramer. These results elucidate the critical role of hydration dynamics in peptide conformational transitions and protein structural stability and integrity.     

Ultrafast dynamics of nonequilibrium resonance energy transfer and probing globular protein flexibility of myoglobin

Protein structural plasticity is critical to many biological activities and accurate determination of its temporal and spatial fluctuations is challenging and difficult. Here, we report our extensive characterization of global flexibility of a globular heme protein of myoglobin using resonance energy transfer as a molecular ruler. With site-directed mutagenesis, we use a tryptophan scan to examine local structural fluctuations from B to H helices utilizing 10 tryptophan-heme energy transfer pairs with femtosecond resolution. We observed ultrafast resonance energy transfer dynamics by following a nearly single exponential behavior in 10-100 ps, strongly indicating that the globular structure of myoglobin is relatively rigid, with no observable static or slow dynamic conformational heterogeneity. The observation is against our molecular dynamics simulations, which show large local fluctuations and give multiple exponential energy transfer behaviors, suggesting too flexible of the global structure and thus raising a serious issue of the force fields used in simulations. Finally, these ultrafast energy transfer dynamics all occur on the similar time scales of local environmental relaxations (solvation), leading to nonexponential processes caused by energy relaxations, not structural fluctuations. Our analyses of such processes reveal an intrinsic compressed- and/or stretched-exponential behaviors and elucidate the nature of inherent nonequilibrium of ultrafast resonance energy transfer in proteins. This new concept of compressed nonequilibrium transfer dynamics should be applied to all protein studies by time-resolved F?rster resonance energy transfer (FRET).     

Ultrafast Fluorescence Spectroscopy via Upconversion and Its Applications in Biophysics

In this review, the experimental set-up and functional characteristics of single-wavelength and broad-band femtosecond upconversion spectrophotofluorometers developed in our laboratory are described. We discuss applications of this technique to biophysical problems, such as ultrafast fluorescence quenching and solvation dynamics of tryptophan, peptides, proteins, reduced nicotinamide adenine dinucleotide (NADH), and nucleic acids. In the tryptophan dynamics field, especially for proteins, two types of solvation dynamics on different time scales have been well explored: ~1 ps for bulk water, and tens of picoseconds for “biological water”, a term that combines effects of water and macromolecule dynamics. In addition, some proteins also show quasi-static self-quenching (QSSQ) phenomena. Interestingly, in our more recent work, we also find that similar mixtures of quenching and solvation dynamics occur for the metabolic cofactor NADH. In this review, we add a brief overview of the emerging development of fluorescent RNA aptamers and their potential application to live cell imaging, while noting how ultrafast measurement may speed their optimization.     

Multiple time scales in solvation dynamics of DNA in aqueous solution: the role of water, counterions, and cross-correlations

Recent time domain experiments have explored solvation dynamics of a probe located inside a DNA duplex, in an effort to gain information, e.g., on the dynamics of water molecules in the DNA major and minor grooves and their environment. Multiple time constants in the range of a few picoseconds to several nanoseconds were obtained. We have carried out 15 ns long atomistic molecular dynamics simulations to study the solvation dynamics of bases of a 38 base-pair long DNA duplex in an aqueous solution containing counterions. We have computed the energy-energy time correlation function (TCF) of the four individual bases (A, T, G, and C) to characterize the solvation dynamics. All the TCFs display highly nonexponential decay with time. When the trajectories are analyzed with 100 fs time resolution, the TCF of each base shows initial ultrafast decay (with tau1 approximately equal 60-80 fs) followed by two intermediate components (tau2 approximately equal 1 ps, tau3 approximately equal 20-30 ps), in near complete agreement with a recent time domain experiment on DNA solvation. Interestingly, the solvation dynamics of each of the four different nucleotide bases exhibit rather similar time scales. To explore the existence of slow relaxation at longer times reported recently in a series of experiments, we also analyzed the solvation TCFs calculated with longer time trajectories and with a larger time resolution of 1 ps. In this case, an additional slow component with a time constant of the order of 250 ps is observed. Through an analysis of partial solvation TCFs, we find that the slow decay originates mainly from the interaction of the nucleotides with the dipolar water molecules and the counterions. An interesting negative cross-correlation between water and counterions is observed, which makes an important contribution to relaxation at intermediate to longer times.     

Temperature dependence of solvation dynamics and anisotropy decay in a protein: ANS in bovine serum albumin

Temperature dependence of solvation dynamics and fluorescence anisotropy decay of 8-anilino-1-naphthalenesulfonate (ANS) bound to a protein, bovine serum albumin (BSA), are studied. Solvation dynamics of ANS bound to BSA displays a component (300 ps) which is independent of temperature in the range of 278-318 K and a long component which decreases from 5800 ps at 278 K to 3600 ps at 318 K. The temperature independent part is ascribed to a dynamic exchange of bound to free water with a low barrier. The temperature variation of the long component of solvation dynamics corresponds to an activation energy of 2.1 kcal mol(-1). The activation energy is ascribed to local segmental motion of the protein along with the associated water molecules and polar residues. The time scale of solvation dynamics is found to be very different from the time scale of anisotropy decay. The anisotropy decays are analyzed in terms of the wobbling motion of the probe (ANS) and the overall tumbling of the protein.     

Probing polar solvation dynamics in proteins: a molecular dynamics simulation analysis

Measurements of time-resolved Stokes shifts on picosecond to nanosecond time scales have been used to probe the polar solvation dynamics of biological systems. Since it is difficult to decompose the measurements into protein and solvent contributions, computer simulations are useful to aid in understanding the details of the molecular behavior. Here we report the analysis of simulations of the electrostatic interactions of the rest of the protein and the solvent with 11 residues of the immunoglobulin binding domain B1 of protein G. It is shown that the polar solvation dynamics are position-dependent and highly heterogeneous. The contributions due to interactions with the protein and with the solvent are determined. The solvent contributions are found to vary from negligible after a few picoseconds to dominant on a scale of hundreds of picoseconds. The origin for the latter is found to involve coupled hydration and protein conformational dynamics. The resulting microscopic picture demonstrates that a wide range of possibilities have to be considered in the interpretation of time-resolved Stokes shift measurements.     

Hydration dynamics and time scales of coupled water-protein fluctuations

We report experimental and theoretical studies on water and protein dynamics following photoexcitation of apomyoglobin. Using site-directed mutation and with femtosecond resolution, we experimentally observed relaxation dynamics with a biphasic distribution of time scales, 5 and 87 ps, around the site Trp7. Theoretical studies using both linear response and direct nonequilibrium molecular dynamics (MD) calculations reproduced the biphasic behavior. Further constrained MD simulations with either frozen protein or frozen water revealed the molecular mechanism of slow hydration processes and elucidated the role of protein fluctuations. Observation of slow water dynamics in MD simulations requires protein flexibility, regardless of whether the slow Stokes shift component results from the water or protein contribution. The initial dynamics in a few picoseconds represents fast local motions such as reorientations and translations of hydrating water molecules, followed by slow relaxation involving strongly coupled water-protein motions. We observed a transition from one isomeric protein configuration to another after 10 ns during our 30 ns ground-state simulation. For one isomer, the surface hydration energy dominates the slow component of the total relaxation energy. For the other isomer, the slow component is dominated by protein interactions with the chromophore. In both cases, coupled water-protein motion is shown to be necessary for observation of the slow dynamics. Such biologically important water-protein motions occur on tens of picoseconds. One significant discrepancy exists between theory and experiment, the large inertial relaxation predicted by simulations but clearly absent in experiment. Further improvements required in the theoretical model are discussed.     

Classical trajectory simulations of photoionization dynamics of tryptophan: intramolecular energy flow, hydrogen-transfer processes and conformational transitions

One-photon and two-photon ionization dynamics of tryptophan is studied by classical trajectory simulations using the semiempirical parametric method number 3 (PM3) potential surface in "on the fly" calculations. The tryptophan conformer is assumed to be in the vibrational ground state prior to ionization. Initial conditions for the trajectories are weighted according to the Wigner distribution function computed for that state. Vertical ionization in the spirit of the classical Franck-Condon principle is assumed. For the two-photon ionization process the ionization is assumed to go resonantively through the first excited state. Most trajectories are computed, and the analysis is carried out for the first 10 ps. A range of interesting effects are observed. The main findings are as follows: (1) Multiple conformational transitions are observed in most of the trajectories within the ultrafast duration of 10 ps. (2) Hydrogen transfer from the carboxyl group to the amino group and back has been observed. A zwitterion is formed as a transient state. (3) Two new isomers are formed during the dynamics, which have apparently not been previously observed. (4) Fast energy flow between the ring modes and the amino acid backbone is observed for both one- and two-photon ionization. However, the effective vibrational temperatures only approach the same value after 90 ps. The conformation transition dynamics, the proton-transfer processes and the vibrational energy flow are discussed and analyzed.     

Ultrafast surface solvation dynamics and functionality of an enzyme alpha-chymotrypsin upon interfacial binding to a cationic micelle

In this contribution we report studies on enzymatic activity of alpha-chymotrypsin (CHT) upon complexation with cationic cetyltrimethylammonium bromide (CTAB) micelle. With picosecond time resolution, we examined solvation dynamics at the interface of CHT-micelle complex, and rigidity of the binding. We have used 5-(dimethyl amino) naphthalene-1-sulfonyl chloride (dansyl chloride; DC) that is covalently attached to the enzyme at the surface sites. The solvation processes at the surface of CHT in buffer solution are found to be mostly in the sub-50 ps time scale. However, at the interface the solvation correlation function decays with time constant 150 ps (65%) and 500 ps (35%), which is significantly different from those found at the enzyme and micellar surfaces. The binding structure of the enzyme-micelle complex was examined by local orientational motion of the probe DC and compared with the case without micelle. The orientational dynamics of the probe DC in the complex reveals a structural perturbation at the surface sites of CHT upon complexation, consistent with other reported structural studies. We also found possible entanglement of charge transfer dynamics of the probe DC on the measured solvation processes by using time-resolved area normalized emission spectroscopic technique. The interfacial solvation process and complex rigidity elucidate the strong recognition mechanism between CHT and the micelle, which is important to understand the biological function of CHT upon complexation with the micelle.     

Probing Deuterium Isotope Effect on Structure and Solvation Dynamics of Human Serum Albumin

The deuterium isotopic effect on the structure and solvation dynamics of the protein, human serum albumin (HSA), has been studied by using circular dichroism (CD), femtosecond up‐conversion, FRET, and single‐molecule spectroscopy. The CD spectra suggest that D₂O affects the structure of HSA, leading to a 20 % decrease in the helical structure. The FRET study indicates that the distance of C153 from the lone tryptophan residue of HSA is quite similar (≈21 Å) in H₂O and D₂O, and hence, the location of the probe in the protein remains the same in the two solvents. The single‐molecule study suggests that coumarin 153 (C153) binds almost exclusively (>96 %) to one site of HSA. Solvation dynamics of C153 in HSA is found to be markedly retarded in D₂O compared with H₂O. In H₂O, the solvation of C153 bound to HSA is found to be biexponential with one component of 7 ps (30 %) and a long component of 350 ps (70 %). In D₂O, we detected a short component of 4 ps (41 %) and a long component of 950 ps (59 %). Thus, the ultraslow component of the solvation dynamics of C153 bound to HSA in D₂O (950 ps) is 2.5‐fold slower than that in H₂O (350 ps). The marked deuterium isotope effect has been ascribed to water molecules confined in the protein environment and to a lesser extent to the structural modification of protein by D₂O.     

Temperature dependence of solvation dynamics of probe molecules in methanol-doped ice and in liquid ethanol

We have studied the solvation statics and dynamics of coumarin 343 and a strong photoacid (pK* approximately 0.7) 2-naphthol-6, 8-disulfonate (2N68DS) in methanol-doped ice (1% molar concentration of methanol) and in cold liquid ethanol in the temperature range of 160-270 K. Both probe molecules show a relatively fast solvation dynamics in ice, ranging from a few tens of picoseconds at about 240 K to nanoseconds at about 160 K. At about 160 K in doped ice, we observe a sharp decrease of the dynamic Stokes shift of both coumarin 343 and 2N68DS. Its value is approximately only 200 cm-1 at approximately 160 K compared to about 1100 cm-1 at T >/= 200 K (at times longer than t > 10 ps). We find a good correlation between the inefficient and slow excited-state proton-transfer rate at low-temperature ice, T < 180 K, and the dramatic decrease of the solvation energy, as measured by the dynamic band shift, at these low temperatures. We find that the average solvation rate in ice is similar to its value in liquid ethanol at all given temperatures in the range of 200-250 K. The surprisingly fast solvation rate in ice is explained by the relatively large freedom of the water hydrogen rotation in ice Ih.     

Understanding dynamics of myoglobin in heterogeneous aqueous environments using coupled water fractions

This work presents an analysis of near environment of myoglobin (Mb) in different aqueous solutions (in the presence of NaCl, sucrose, trehalose, urea, and glycerol) using the coupled water fractions measured using a quartz crystal microbalance (QCM). The secondary structural features of the protein from circular dichroic (CD) spectroscopy and the coupled water fractions give important clues to the overall dynamics of the protein. Using time resolved fluorescence, these leads have been applied to understand the observed lifetime relaxations of Mb. Though the time scales of observation of coupled water and the lifetimes are very different, our study suggests that the trends in coupled water fraction seem to be good indicators for regulation of the relaxation dynamics of the protein. The relaxations generally show a triphasic distribution of time scales. The initial relaxation in the picoseconds time scale represents the local motions of coupled water followed by a slightly slower decay in hundreds of picoseconds attributable to coupled water-‘quasi free’ water interactions. The third nanosecond lifetime is due to changes in transitions in isomers of hydrated protein. The dynamics of coupled water in Mb with NaCl is the fastest (around 21 ps) and is slowest in glycerol (250 ps). The results strongly indicate that it is the resident times of water molecules that play a dominant role in the overall stability of protein in a particular hydrated isomer and not just always the number of such water molecules in the hydrated protein.     

光谱法研究喜树碱与牛血清白蛋白的相互作用

采用多种光谱技术对喜树碱和牛血清白蛋白的相互作用进行了研究.结果表明喜树碱和牛血清白蛋白可形成基态复合物,引起牛血清白蛋白内源荧光猝灭.通过计算获得了二者在不同温度下的结合常数及结合位点数.根据喜树碱和牛血清白蛋白结合的热力学参数,确定了二者之间主要为疏水作用力.根据F(o)rster非辐射能量转移理论确定了喜树碱和牛血清白蛋白的作用距离.同步荧光光谱显示喜树碱主要与蛋白中色氨酸残基发生相互作用,改变其周围的局部构象.红外光谱提示喜树碱可引起蛋白的构象发生改变,α-螺旋二级结构减少.     

Ultrafast dynamics of resonance energy transfer in cryptochrome

In this communication, we report the ultrafast dynamics of resonance energy transfer in a blue-light photoreceptor, Vibrio cholerae cryptochrome. The transfer was observed to occur in 60 ps. We also studied the local rigidity and solvation around the binding site of the photoantenna molecule. The results for the first time show energy transfer in cryptochrome suggesting some mechanistic similarities between photolyase that repairs damaged DNA and cryptochrome that mediates blue-light signaling.     

Isothermal titration calorimetric studies of the pH induced conformational changes of bovine serum albumin

Bovine serum albumin (BSA) is a soft globular protein that undergoes conformational changes through several identified transition steps in the pH range 2–13.5. The ability to change conformation makes BSA capable of complexing different ligands from fatty acids to cations or drugs and carries them in the bloodstream. Microcalorimetric titration of BSA with NaOH solution was performed to measure the enthalpy of conformational changes. Two exothermic enthalpy changes were found in the course of the titration between pH 3 and 9.5, which can be identified with the E–F, and the F–N transitions. The enthalpy change at pH 3.5 (transition from the E to the F form of BSA, folding of intra-domain helices in domain I) is independent of the protein concentration. The second transition (F–N, folding of domain III) was observed at pH 4.8 for the 0.1% BSA solution, but it shifted to higher pH values as the protein concentration increased to 0.2% and 0.3%. The tightening of the protein structure with increasing pH was verified measuring intrinsic fluorescence of tryptophan residues. At even higher pH value, pH 10.5, fluorescence measurements revealed protein expansion. The BSA conformational changes were also measured by dynamic light scattering. The hydrodynamic diameter was smaller at the i.e.p. of BSA (5–7 nm at pH ~5) and larger at the two ends of the pH range (17.5 nm at pH 2 and 8.3 nm at pH 10).     

Ultrafast quenching of tryptophan fluorescence in proteins: Interresidue and intrahelical electron transfer

Quenching of tryptophan fluorescence in proteins has been critical to the understanding of protein dynamics and enzyme reactions using tryptophan as a molecular optical probe. We report here our systematic examinations of potential quenching residues with more than 40 proteins. With site-directed mutation, we placed tryptophan to desired positions or altered its neighboring residues to screen quenching groups among 20 amino acid residues and of peptide backbones. With femtosecond resolution, we observed the ultrafast quenching dynamics within 100 ps and identified two ultrafast quenching groups, the carbonyl- and sulfur-containing residues. The former is glutamine and glutamate residues and the later is disulfide bond and cysteine residue. The quenching by the peptide-bond carbonyl group as well as other potential residues mostly occurs in longer than 100 ps. These ultrafast quenching dynamics occur at van der Waals distances through intraprotein electron transfer with high directionality. Following optimal molecular orbital overlap, electron jumps from the benzene ring of the indole moiety in a vertical orientation to the LUMO of acceptor quenching residues. Molecular dynamics simulations were invoked to elucidate various correlations of quenching dynamics with separation distances, relative orientations, local fluctuations and reaction heterogeneity. These unique ultrafast quenching pairs, as recently found to extensively occur in high-resolution protein structures, may have significant biological implications.     

Fluorescence Study of the Conformational Properties of Recombinant Tick Anticoagulant Peptide (Ornithodorus moubata) Using Multifrequency Phase Fluorometry

Abstract— Steady-state and multifrequency phase fluorometry were used to characterize the conformational state and conformational dynamics of recombinant tick anticoagulant peptide ( Ornithodorus moubata ) (TAP). The TAP contains two tryptophan residues at positions 11 and 37. The fluorescence emission varies sigmoidally as a function of pH with a pK_a of 6.01 ± 0.07. This pH dependency suggests that tryptophan fluorescence is quenched by His43 at low pH. This is confirmed by modification of the his-tidine with diethylpyrocarbonate. At pH 9 the fluorescence decay is well described by a sum of three exponentials (0.52,1.9 and 5.4 ns), which decrease all three at pH 4 (0.25, 1.61 and 4.4 ns). From the reactivity of the fluorescence lifetimes toward N -bromosuccinimide and from the calculation of the accessibility we can attribute the long lifetime to Trpll, the short one to Trp37 and the middle one to both. The anisotropy decay was resolved into two components of 3.85 ns and 0.27 ns at pH 4 and 4.5 ns and 0.6 ns at pH 9. The long anisotropy decay time corresponds to the rotational correlation time of the protein, the short one to local mobility of the tryptophan residues.     

Effect of ionic liquid on the native and denatured state of a protein covalently attached to a probe: Solvation dynamics study

Effect of a room temperature ionic liquid (RTIL, [pmim][Br]) on the solvation dynamics of a probe covalently attached to a protein (human serum albumin (HSA)) has been studied using femtosecond up-conversion. For this study, a solvation probe, 7-diethylamino-3-(4-maleimidophenyl)-4-methylcoumarin (CPM) has been covalently attached to the lone cysteine group (cys-34) of the protein HSA. Addition of 1.5 M RTIL or 6 M GdnHCl causes a red shift of the emission maxima of CPM bound to HSA by 3 nm and 12 nm, respectively. The average solvation time ?τ(s)? decreases from 650 ps (in native HSA) to 260 ps (～2.5 times) in the presence of 1.5 M RTIL and to 60 ps (～11 times) in the presence of 6 M GdnHCl. This is ascribed to unfolding of the protein by RTIL or GdnHCl and therefore making the probe CPM more exposed. When 1.5 M RTIL is added to the protein denatured by 6 M GdnHCl in advance, a further ～5 nm red shift along with further ～2 fold faster solvent relaxation (?τ? ～30 ps) is observed. Our previous fluorescence correlation spectroscopy study [D. K. Sasmal, T. Mondal, S. Sen Mojumdar, A. Choudhury, R. Banerjee, and K. Bhattacharyya, J. Phys. Chem. B 115, 13075 (2011)] suggests that addition of RTIL to the protein denatured by 6 M GdnHCl causes a reduction in hydrodynamic radius (r(h)). It is demonstrated that in the presence of RTIL and GdnHCl, though the protein is structurally more compact, the local environment of CPM is very different from that in the native state.