Probing Single-Molecule Protein Conformational Dynamics

Feb 10, 2005 - H. PETER LU*. Pacific Northwest National Laboratory, Fundamental Science. Division, P.O. Box 999, Richland, Washington 99352. Received ...
1 downloads 0 Views 386KB Size
Acc. Chem. Res. 2005, 38, 557-565

Probing Single-Molecule Protein Conformational Dynamics H. PETER LU* Pacific Northwest National Laboratory, Fundamental Science Division, P.O. Box 999, Richland, Washington 99352 Received June 15, 2004 ABSTRACT Protein conformational fluctuations and dynamics, often complex and associated with inhomogeneities, play a crucial role in biomolecular functions. It is extremely difficult to characterize such inhomogeneous dynamics in an ensemble-averaged measurement, especially when the proteins are involved in multiple-step, multipleconformation complex chemical interactions and transformations, such as in enzymatic reactions, protein-protein interactions, and ion-channel membrane protein processes. Alternatively, singlemolecule spectroscopy is a powerful approach to probing and analyzing protein conformational dynamics in real time.

Introduction Protein molecules typically engage in complex and dynamic interactions with other protein species or local components of their environments. Protein conformations are highly dynamic rather than static, and subtle conformational changes play a crucial role in protein function. It has been widely recognized that knowledge about static structure alone is not sufficient to understand real-time protein functions1-15 and that the conformational changes and ordered-disordered protein structure transitions6,7,16,17 play important roles in determining protein functions. Nevertheless, static structure analyses from ensembleaveraged measurements at equilibrium have been the norm in understanding of protein structure and functions. Recently, however, single-molecule spectroscopy has been demonstrated as a powerful complementary approach in characterizing dynamic conformations and understanding correlated protein functions in real time. Single-molecule studies allow characterization of temporal, spatial, and energetic properties of protein conformations, which involve complex interactions among protein molecules, enzyme-substrate interactions and enzymatic reactions, and ion channel processes in membranes. Single-molecule studies can thus result in a detailed, mechanistic understanding of some of the most important biological processes, including receptor sensing, cell redox respirations, cell signaling, and enzymatic reactions. Several elegant articles1-5 have provided extensive reviews of singlemolecule spectroscopy and its applications. Single-molecule fluorescence photon-stamping time trajectory detection6,9,18-21 and statistical analysis have H. Peter Lu received his B.S. from Peking University in 1982 and his Ph.D. degree in physical chemistry from Columbia University in 1991. After postdoctoral work at Northwestern University and the Pacific Northwest National Laboratory, he joined the Pacific Northwest National Laboratory as a senior research scientist in 1996. He is currently a senior research scientist II with primary research interests in single-molecule spectroscopy and microscopy, single-molecule protein conformational dynamics, and interfacial chemical and biophysical dynamics. 10.1021/ar0401451 CCC: $30.25 Published on Web 02/10/2005

 2005 American Chemical Society

revealed static and dynamic inhomogeneities in the dynamics of a number of protein systems.1-4,6-9,11-14 Static inhomogeneity describes inhomogeneous molecular dynamics from molecule to molecule, whereas dynamic inhomogeneity originates in the temporal fluctuations of the kinetic rates for individual molecules, which is beyond the scope of conventional kinetics.1,4,6-9,13,14,20-23 It is nearly impossible for an ensemble-averaged measurement to selectively identify and characterize both static and dynamic inhomogeneities. The following Account focuses on three areas of our recent work on single-molecule conformational dynamics: (1) the enzyme-substrate complex formation in T4 lysozyme hydrolization enzymatic reactions,8-10 (2) gramicidin-dimer ion channels at lipid bilayers,11,12 and (3) protein-protein interactions in a cell signaling system.6,18

Enzyme Conformational Dynamics in Catalytic Reactions Enzymatic reactions typically involve multiple kinetic steps, complex molecular interactions, inhomogeneous conformational changes, and confined local environments. These complexities, associated with intrinsic static and dynamic inhomogeneities, often make an ensembleaveraged measurement alone insufficient for characterizing enzymatic reaction dynamics. To obtain a molecularlevel understanding of the enzymatic reaction dynamics, conformational dynamics, intermediate enzyme-substrate complexes, and energy landscape of specific conformational changes along reaction coordinates,8 our approaches have combined single-molecule spectroscopy, molecular dynamics (MD) simulation, and theoretical modeling. One object of study has been the T4 lysozyme hydrolysis reaction of a polysaccharide.8-10 The wild-type T4 lysozyme, having two domains connected by an R-helix, adopts different conformations as the two domains undergo a hinge-bending motion.24 Previously, there have been no real-time observations of how much an active site opens to initiate the formation of an enzymesubstrate complex under physiological conditions, although the literature shows that the scale of the motion can be 4 to 8 Å.24-26 Furthermore, except for an assessment of the overall reaction rates, other features of the process, e.g., the multistep reaction dynamics, the enzymesubstrate intermediate states, and the correlated conformational dynamics, have also been unknown because there were no adequately selective and temporal characterization methods for such complex reaction systems. By covalently attaching a fluorescence resonance energy transfer (FRET)8,27 donor-acceptor TMR/Texas Red pair site-specifically to noninterfering sites, i.e., two cysteine residues on the enzyme (Figure 1A, inset), we were able to measure the hinge-bending motions of the enzyme by recording the time trajectories of the single-pair donor* E-mail: [email protected]. VOL. 38, NO. 7, 2005 / ACCOUNTS OF CHEMICAL RESEARCH

557

Single-Molecule Protein Conformational Dynamics Lu

actual donor-acceptor dipole-dipole distance change is estimated to be ∼5.5 Å. At the Fo¨rster distance R0 of 50 ( 5 Å, the distance change of 5.5 Å can cause a change of 2∼3 times in the intensity of donor fluorescence, although the change in the acceptor fluorescence intensity is less than 10%.8 The upper panel in Figure 1A shows an expanded portion of a trajectory (middle panel) recorded from donor fluorescence of a single-pair donor-acceptor labeled protein with substrate present. By comparison, the lower panel shows a portion of a donor-fluorescence trajectory recorded from a donor-only labeled T4 lysozyme, not sensitive to the enzyme open-close hinge-bending motions, under the same conditions. The large-amplitude, lower-frequency wiggling of the donor fluorescence intensity in the upper panel, with a bimodal fluorescence intensity distribution (inset), reflects the open and closed conformational states of a T4 lysozyme.8

FIGURE 1. (A) Probing of a single T4 lysozyme enzymatic reaction turnover trajectory. The upper panel shows an expanded portion of a trajectory from donor fluorescence of a donor-acceptor-labeled single T4 lysozyme under an enzymatic reaction. The lower panel shows a trajectory of a donor-alone-labeled enzyme. The fluorescence intensity distributions derived from the two trajectories are shown in the insets of the middle and lower panels. The solid lines are fit using bimodal and Gaussian functions, respectively. (Inset, top-right) Crystal structure of wild-type T4 lysozyme with donoracceptor labeling at Cys-54 and Cys-97. (B) Molecular dynamics simulation of T4 lysozyme hinge-bending motion. Three conformation trajectories of the distance between two -SH of Cys-54 and Cys97 of T4 lysozyme in solution at room temperature are presented: free enzyme without substrate (E) (left panel), nonspecific binding of polysaccharide (E + S f ES) (middle panel), and active enzymepolysaccharide complex (ES*) (right panel). Adapted with permission from ref 8. Copyright 2003 American Chemical Society. acceptor emission intensities. We characterized singlemolecule FRET and anisotropy dynamics over a wide range of time-scales, both the nanosecond fast dynamics and the sub-millisecond or longer slow dynamics.8,9 In our experiments, the dye-labeled T4 lysozyme proteins were covalently tethered on a hydrocarbon-modified glass surface through a bifunction-linker molecule under buffer solution.8,9 Single-molecule FRET fluorescence trajectories contain detailed information about the conformational motions associated with the enzymatic turnovers. Using MD simulation, we have calculated that the distance change between the two dye-tethered cysteine residues is from 30.5 to 35 Å.8 Considering the lengths of the two covalent linking groups of the donor and acceptor dipoles, the 558 ACCOUNTS OF CHEMICAL RESEARCH / VOL. 38, NO. 7, 2005

We have attributed each wiggle to a hinge-bending motion, i.e., an open-closed cycle, involved in either an enzymatic reaction turnover or nonproductive binding and releasing of the substrate, based on a series of ensemble-averaged and single-molecule control experiments.8 The enzyme and polysaccharide chains engage in a nonspecific interaction through diffusion to form a nonactivated complex state, ES. Then, the active-site cleft binds to the polysaccharide chain to form specific hydrogen bonds in the activated complex state, ES*. The subsequent release of product and the enzyme searching for the next reactive site in the substrate complete the enzymatic turnover cycle. The donor fluorescence intensity increases as the active site opens due to substrate inserting and decreases as the active site closes to form an active enzyme-substrate complex. This attribution is consistent with the results of the ensemble-averaged FRET measurements made with and without a substrate.8 Further evidence that we are measuring the hinge-bending motion comes from evaluating autocorrelation functions calculated from the single-molecule donor and acceptor trajectories when changing the laser excitation intensity, the pH, and the substrate concentration.8 Single-molecule FRET is sensitive to conformational change dynamics but not necessarily accurate in determining exact donor-acceptor distance changes of the conformations due to spectral fluctuation28 and dye photophysics.29 To evaluate quantitatively the amplitude of the protein-domain open-closed motions in forming an activated enzyme-substrate complex, we applied MD simulation to explore the conformations of the intermediate states of the enzymatic reaction.8 MD simulation is intrinsically a single-molecule computational experiment involving picosecond (ps) to nanosecond (ns) short-time trajectories as “snapshots,” as opposed to the millisecond time-resolved single-molecule experimental trajectories. Our MD simulation (Figure 1B) showed that when the substrate enters the active site, forming ES, the two domains of the enzyme exhibit a breathing-type hingebending motion along the R-helix, causing a 4.5 ( 0.2 Å distance increase (〈XN(t)〉/2) between the two cysteines.8

Single-Molecule Protein Conformational Dynamics Lu

on this model analysis, we were able to obtain the diffusion coefficient, D ) 3.8 × 10-14 cm2/s, of the T4 lysozyme substrate-enzyme complex formation conformational motions from measurable parameters, both the single-molecule spectroscopy data and MD simulation results:8

D)

FIGURE 2. (A) Active-complex, ES*, formation time (topen) distribution deduced from single T4 lysozyme fluorescence trajectories under enzymatic reaction. The topen is the duration time of each wiggling of the intensity trajectories above a threshold. (B) Energy potential surface of T4 lysozyme-substrate complex formation process. The conformational change dynamics involving multiple intermediate states is analyzed based on a one-dimensional random walk model coupled with the parameters obtained from our single-molecule experimental spectroscopy and MD simulation. Since the total energy change between E and ES* states is about -18 kcal/mol, assuming six intermediate states would result in an average energy difference of 3 kcal/mol for each associated conformational state along the R-helix coordinate during the hinge-bending motion. We postulate that the activation energy associated with the kf of the forward step is about 0∼5 kcal/mol, considering the slow forward rate and the entropy decrease in the complex formation process. It is reasonable to assume the energy potential surface along the conformational change nuclear coordinate to be parabolic for each intermediate state. Therefore, the averaged force constant of the potential surface for each intermediate state is calculated to be 3.4∼20 (kcal/mol)/Å2. Adapted with permission from ref 8. Copyright 2003 American Chemical Society. Experimentally, the formation time, topen, of ES and ES* is probed in real time by recording single-molecule FRET trajectories.8 Figure 2A shows a distribution of the opentime (topen) deduced from single-molecule trajectories. The first moment of the distribution, 〈topen〉 ) 19.5 ( 2 ms, corresponds to the mean formation time of the ES* state. The standard deviation of the distribution, x〈∆topen2〉 ) 8.3 ( 2 ms, reflects the distribution bandwidth. We model the hinge-bending motion associated with interactions between the enzyme and substrate as a classical particle one-dimensional multiple-step random walk in the presence of a force field (Figure 2B).8,30 Based

(x〈∆topen2〉〈XN(t)〉)2 2〈topen〉3

We also obtained the mean drifting velocity, 〈v〉 ) 〈XN(t)〉/〈topen〉, of the conformational change, 4.6 × 10-6 cm/s. We further characterized the energy landscape of the hinge-bending conformational change dynamics by estimating the minimum number of the intermediate conformational states involved in the complex formation and by calculating the averaged rates of forming these conformational states. From the random walk model (Figure 2B upper panel),8 we have D/〈v〉 ) L[(kf + kb)/(2(kf - kb))], and L ) 2D/〈v〉 ) 1.6 Å, when kb f 0; therefore, the minimum number of intermediate conformation states is m ) 〈XN(t)〉/L ) 5.6. The average rate for each step is m/〈topen〉 ) 280 s-1. This result (m > 2 at the limit of kb f 0) suggests that there are definitely more than two conformational intermediate states in addition to ES and ES*.8 Our model analysis, based on both experimental results and MD simulation, also identified the friction coefficient, the energy consumed by friction along the reaction coordinates, and the averaged force constant of the potential surface for each intermediate state (Figure 2B).8 The information about the existence of the multiple intermediate conformation states involving the enzymatic active complex formation, an experimental characterization of the energy landscape (Figure 2B), and dynamics of the complex formation process cannot be obtained solely by either an ensemble-averaged experiment, a single-molecule experiment, or a computational MD simulation approach. For studies of protein conformational dynamics, a combined approach, using both singlemolecule spectroscopy and MD simulations, is essential. Using the methodology described above, which has allowed a better understanding of hydrolysis enzymatic reaction in dissolving cell walls, we are now studying the spatially inhomogeneous and temporally intermittent single-molecule enzymatic reactions in living cell surfaces. Using single-molecule injection,10 we directly place a single-molecule T4 lysozyme on a cell wall surface and watch the highly effective enzyme proteins in action.

Probing Multiple-State Conformational Dynamics of Single Ion-Channel Membrane Proteins Ion channels are proteins that control ion flux across the cell membranes. Subtle conformational changes of ionchannel proteins are often critical in the functioning of ion channels. However, the conformational dynamics for VOL. 38, NO. 7, 2005 / ACCOUNTS OF CHEMICAL RESEARCH 559

Single-Molecule Protein Conformational Dynamics Lu

FIGURE 3. Experimental apparatus and measurement configuration of PCFM. Lipid bilayers were formed at the tip of a patch pipet by apposition of two monolayers of mixed lipids (4:1 dPhPE:DPhPC). Dye-labeled gramicidin monomers are introduced from both sides of the bilayer, and single-channel electric currents are recorded by using typical patch-clamp instrumentation. The wide-field or confocal excitation was focused on the pipet tip, and the fluorescence was imaged on a CCD or a pair of photostamping modules. Adapted with permission from ref 11. Copyright 2003 American Chemical Society. most ion-channel proteins are still largely unknown due to lack of direct measurements in real time. The patchclamp recording technique,31 which records ion currents of cations or anions such as sodium ions and protons flowing across biological membranes through ion channel proteins, has had a major impact on the understanding of ion-channel processes. However, the conformational dynamics of the ion channel proteins has been primarily extrapolated by model analyses of the electric current trajectories that are not particularly sensitive to subtle conformational intermediate states, including electrically undetectable “silent” conformational states. In 1999, Yanagida and co-workers32 reported singlechannel fluorescence images acquired before or after a single-channel current measurement on the same sample lipid bilayer, utilizing a micro-pinhole patch technique. During that time and since then, there have been significant efforts and advancements.11,12,33-41 These research efforts have begun to provide a new paradigm for studying ion-channel conformational dynamics and mechanisms and for obtaining information not obtainable by the conventional patch current recording measurements, a primary experimental approach for decades. The following discussion is focused on the new scientific information obtained in our single-molecule studies on the conformation-gated dynamics of the gramicidin ion channels in lipid bilayers by using a new approach,11,12,41 patch-clamp fluorescence microscopy (PCFM) (Figure 3), which combines simultaneous single-channel electric current patch recording with single-molecule fluorescence imaging of co-localization and FRET, nanosecond lifetime and anisotropy, and fluorescence selfquenching. Linear gramicidin, produced naturally from the bacterium Bacillus brevis, is probably one of the mostly studied ion channels.42-44 According to the literature, there 560 ACCOUNTS OF CHEMICAL RESEARCH / VOL. 38, NO. 7, 2005

are two conformational states involved in the channel open-closed activity. Its conductive and open state is a dimer of two 15-amino acid peptides forming β-helices, which meet head-to-head at the N-termini in the interior of a lipid bilayer; a dissociated dimer constitutes a closed state. However, even for this “simple” ion channel, the nature of the structural changes and the conformational dynamics that regulate the open and closed states are still being debated.44 Our experimental results, not obtainable by using ensemble-averaged analyses or traditional singlechannel patch electric recording, revealed that multiple intermediate conformational states are involved in the channel’s open-closed activities. Such results advanced our knowledge of ion channel dynamics and significantly modified the conventional two-state model of gramicidin. Using PCFM, we were able to record simultaneously the single-channel electric current and spectroscopic images and thus identify and characterize the single ionchannel protein conformational dynamics. Asymmetric incorporation of Cy5-gramicidin C and TMR-gramicidin C11,12 led to the formation of single channels tagged with pairs of FRET donors and acceptors. Co-localization and FRET within a single gramicidin heterodimer were probed by dual-color excitation and donor-acceptor two-channel fluorescence imaging (Figure 4A) and by single-molecule photon time-stamping.6,9,11,12,18-21 Dye molecule spin or wobbling motions were at a nanosecond time-scale8,12,18 and averaged out in the FRET imaging measurements at a millisecond time-scale.11,12 The mean spFRET efficiency of 0.59 ( 0.05 and the fwhm of 0.21 (Figure 4B) correspond to a mean donor-acceptor distance of ∼56 Å and a range of ∼52 Å e d e ∼63 Å. Dihedral angular distortion of the two monomers due to hydrogen bonding destabilization11,12 may also contributed to the broad distribution. Nevertheless, FRET results suggest the existence of the multiple intermediate conformational states, beyond the

Single-Molecule Protein Conformational Dynamics Lu

FIGURE 4. PCFM measurement of single TMR/Cy5-gramicidin heterodimer ion channels. The “green” and “red” arrows mark the timing of the toggled excitation at 514 and 632 nm, respectively. (A) Single-channel electric current trajectory (upper panel) with four arrows indicating the exposure of four consecutively acquired FRET imaging frames. Frames 1 and 3 were taken at 514-nm excitation, and Frames 2 and 4 were taken at 632-nm excitation. Frame 1 indicates that FRET occurred between the TMR/Cy5 spFRET pair when the gramicidin channel was at an open state. Frame 2 verifies that a single Cy5 molecule was present. Frame 3 indicates that no FRET occurred when the channel was at a closed state. Frame 4 verifies that the Cy5 was still present. (B) FRET efficiency distribution of “channel open” (shaded bars) and “channel closed” (open bars). It is experimentally difficult to differentiate FRET below the efficiency of 0.3 or no FRET (indicated by *). Adapted with permission from ref 11. Copyright 2003 American Chemical Society. conformational fluctuation within a conformational state or between two conformational states of open-closed states.11,12 The broad distributions of FRET efficiency, correlated with both open and closed states, suggest that the channels could be in closed states even when the two gramicidin monomers are still in an intermediate nondissociated state. Figure 5A shows a single-channel electric current trajectory simultaneously recorded with a synchronized fluorescence-quenching time trajectory. The fluorescence intensity and the electric current trajectories show a distinctive anti-correlated pattern when the channel is in an open state due to the close proximity of the pair of fluorophores across the lipid bilayer.12 A further analysis (Figure 5A-C) shows that there are statistically significant time mismatches between the two trajectories in the channel-opening and -closing events. We attributed the time mismatches to conformation-associated activity induction times or the silent times, ∆1 and ∆2, between the channel conformational changes and the channel opening and closing events, respectively (Figure 5A).12 The distributions of the stochastic ∆1 and ∆2 are essentially

FIGURE 5. Identification and distribution of the conformationassociated silent time. The fluorescence self-quenching of a TMRdye pair on a gramicidin dimer causes an intensity decrease. (A) Identification of the silent times, the time lags between the channel electric current changes, and fluorescence intensity changes. (B) The distribution of the silent times. (C) The collective distribution of the silent times. (D) Proposed multi-state model of gramicidin dynamics indicating a third “multi-state” or a set of intermediate conformation states of nonconducting or partially conducting channels with larger separation between the fluorophore probes (middle). Poisson and can be fitted with single-exponential decays (Figure 5B and 5C). The mean time of conformational change among the intermediate states and open-closed states was determined to be 200 ( 50 ms for ∆1 and 400 ( 50 ms for ∆2. The multiple intermediate states are involved in gramicidin channel dynamics, in addition to the channel open and closed states, and their interconversion time is at a subsecond time-scale. Based on our experimental results,11,12 we proposed a new model of multi-state gramicidin channel dynamics (Figure 5D). The open (upper left), closed (upper right), and intermediate (middle) states can presumably be the outcomes of intra- and intermolecular conformational changes. There is a subset of the intermediate states. We postulate that the intermediate conformers result predominantly from the six intermolecular H-bond fluctuaVOL. 38, NO. 7, 2005 / ACCOUNTS OF CHEMICAL RESEARCH 561

Single-Molecule Protein Conformational Dynamics Lu

tions that disrupt the pore structure of the gramicidin dimer.11 Lately, this postulated mechanism has been supported by MD simulations.45 Subtle conformational changes at the N-terminal linkage of the dimer may significantly disrupt the channel cavity structure and the electrostatic field within the channel, thereby perturbing the channel conductance. It is also noticeable that the gramicidin conformational dynamics are coupled with the peptide-lipid interactions.11,12 The unique real-time imaging and spectroscopy of the dynamics of individual gramicidin ion channels enabled us to reveal a hitherto undetected correlation of multiple conformational states of the channel consistent with transitions among closed, intermediate, and open states. The methodology we developed should be applicable for studying ion-channel protein conformational dynamics and mechanisms in more complex membrane protein systems, including the ion channel proteins that involve strong protein-lipid interacting modulations. It is now possible to explore the function-regulating and sitespecific conformational changes of single ion-channel membrane proteins in living cells and under physiological conditions.41

Protein-Protein Interaction Dynamics in Cell Signaling Proteins rarely act alone but often involve complex protein-protein interactions. Most of the important biological processes, including cell signaling, gene expression, enzymatic metabolism, and bacterial redox respiration, have crucial contributions from protein-protein interactions. Most often, protein-protein interactions are intrinsically single-molecule processes in living cells. For example, cell-signaling processes can be initiated through a few copies of interacting protein complexes, which determine pathway branching and are amplified along the signaling pathway. New information obtained from our single-molecule spectroscopy studies6,7 shows that protein-protein interactions in cell signaling6 and DNA damage recognition7 may involve highly flexible or disordered tertiary structures and slow conformational fluctuations at time-scales of sub-milliseconds to seconds. The large-amplitude conformational fluctuations underlie significant static and dynamic inhomogeneities in the interactions of protein complexes.6,7 To a certain extent, recent NMR static and ensemble-averaged structural analyses of protein complexes have also indicated that binding domains often undergo dramatic conformational changes, from disordered to ordered tertiary structures, upon induced-fit binding complex formation.16,17 It is reasonable to suggest that structural transitions of flexible conformational domains, which can be identified and characterized in detail by measuring single-molecule conformational fluctuation dynamics, are likely common in biomolecular recognition processes.6,7 We have applied single-molecule spectroscopy to study the interactions of an intracellular signaling protein Cdc42, 562 ACCOUNTS OF CHEMICAL RESEARCH / VOL. 38, NO. 7, 2005

a monomeric GTPase protein, with its downstream effector protein WASP (Wiskott-Aldrich Syndrome Protein).6 In this work, a WASP fragment, 13-kDa CBD (Cdc42 binding domain of WASP) that binds only the activated Cdc42, was used and labeled with a solvatochromic dye (I-SO)6 to probe the hydrophobic interactions at the protein-protein interface that are significant to Cdc42/WASP signaling recognition. Cdc42 acts as a molecular switch in signaling pathways to regulate various cellular responses.46 Cdc42 GTPase can be activated by binding to GTP, and the activated Cdc42 then binds and activates a series of effector proteins via direct and consecutive proteinprotein interactions.47,48 Our single-molecule spectroscopy study revealed both static and dynamic inhomogeneities on the conformational fluctuations of the protein complex involving bound and loosely bound states of the proteinprotein interaction complex.6 Collecting fluorescence emission images and photonstamping time trajectories (Figure 6A) from individual Cdc42-CBD complexes immobilized in an agarose gel (0.5%) thin film, we observed a more than 3-fold increase in fluorescence intensity fluctuations at a time-scale ranging over several orders of magnitude, from milliseconds to sub-seconds. To identify the origination of the fluorescence fluctuations, we conducted a series of control experiments and reached the following conclusions:6 (1) The measured conformational fluctuations were spontaneous, as the decay rate constants of the autocorrelation functions were independent of the laser excitation intensity at a level less than 250 W/cm2 under our experimental conditions. (2) Fluorescence fluctuation did not originate from dye interactions in the CBD biosensor because dye-labeled CBD alone did not give high enough fluorescence intensity to be detected at the single-molecule level (Figure 6B). Only the Cdc42-CBD complex was “visible” to our single-molecule measurements. (3) Fluorescence fluctuation was not a result of interaction between the dye-labeled CBD and the agarose gel because the fluorescence quantum yield of dye-labeled CBD remained as low in the agarose gel as in the buffer solution. (4) A triplet-state formation of the dye molecule in the Cdc42-CBD complex was unlikely to produce such fluorescence fluctuation because the fluctuation is spontaneous rather than photoinduced; besides, the triplet state dynamics would have been at a much faster timescale28,49 under ambient conditions than the 10-to-100ms time-scale observed in our fluorescence intensity fluctuations. The single-molecule control experiments6 and the ensemble-averaged measurements6,48 both support attributing dye fluorescence fluctuations to proteinprotein interactions and to conformational fluctuations at the protein-protein interface being probed.6 We attributed the dynamics of the conformational fluctuations to both bound (B) and loosely bound (LB) states of the protein complex.6 The LB states were presumably a subset of conformations that deviated from the bound equilibrium states, probed by the dye, without disrupting the sub-nanometer long-range interactions, so the overall protein complex is still associated. Compared

Single-Molecule Protein Conformational Dynamics Lu

FIGURE 6. Single-molecule spectroscopy and imaging study of Cdc42-CBD binding dynamics. (A) Single-molecule photon time-stamping trajectory. Each dot corresponds to a photon stamped with an arrival time (t) and a delay time (∆t). The delay time is the duration from the time of excitation of the single molecule to the emission by the same molecule, which is a one-event measurement of the fluorescence lifetime of this molecule. (B) Ensemble fluorescence assays of Cdc42 with CBD biosensor using fluorescence spectral measurements. Fluorescence spectra were obtained with excitation at 568 nm. The solid curve represents the active GTP-activated Cdc42 forming complexes with dye-labeled CBD; the dotted curve represents the GDP-loaded Cdc42 with dye-labeled CBD; and the dash-dotted curve represents the dye-labeled CBD alone. A non-hydrolyzable GTP analogue, GTP-γ-S, was used to lock Cdc42 in the active conformation. Thus, the GTP binding and unbinding process was eliminated from the single-molecule experiments that measured only the activated Cdc42 interacting with CBD. The effects of GTP-γ-S are the same as that of GTP, and the biological relevance and validity of using GTP-γ-S has been well established in the literature.48 (C) The autocorrelation function, C(t), calculated from an intensity trajectory of a single Cdc42-CBD complex, with biexponential decays of 250 ( 60 s-1 and 45 ( 10 s-1. (D) The occurrence histogram of the single-complex conformational fluctuation rates. Adapted with permission from ref 6. Copyright 2004 American Chemical Society.

with the B state, the LB state gives significantly lower fluorescence intensity, as the distorted protein-protein interaction interface becomes more solvent-accessible and hydrophilic. The spectroscopic characterization of the B and LB states showed them to be analogues of CDC42CBD and CBD-alone states measured in a control ensemble-averaged experiment (Figure 6B).6 Our MD simulation demonstrated that the solvent-accessible surface area of the dye molecule undergoes a decrease from ∼377 Å2 to ∼225 Å2 upon Cdc42-CBD binding.6 One of the most intriguing observations is that by using autocorrelation function analyses conformational fluctuation rates were found to be highly inhomogeneous (Figure 6C and 6D). Variations of more than 2 orders of magnitude occurred in the conformational fluctuation rates among individual protein complexes under the same solution condition. Although it is still difficult to identify exactly how many conformational states contribute to the inhomogeneous distribution, we postulated that there are at least two subgroups of states associated with conformational fluctuations at ∼10 s-1 and ∼40 s-1 (Figure 6D). About 25% of the autocorrelation functions of the singlecomplex fluorescence intensity trajectories demonstrated biexponential decays, indicating non-Poisson kinetics.6 The non-Poisson behavior suggests that the Cdc42-CBD interactions have both static and dynamic inhomogeneities.

The real-time observations of the large-amplitude and millisecond conformational fluctuations of the proteinprotein interactions in a cell signaling system can only be achieved by single-molecule spectroscopy, as the stochastic conformational fluctuations are typically averaged out in an ensemble-averaged measurement. The millisecond slow dynamic fluctuation of a highly flexible domain likely plays a crucial role in the protein-protein recognition complex, showing that the induced-fit proteinprotein interactions are highly dynamic rather than static.6,7 Is the slow conformational fluctuation essential for cell-signaling protein-protein interactions? A definitive answer has yet to be obtained and may be obtained by combined biological mutation, single-molecule spectroscopy, MD simulation, theoretical modeling, and structural analyses by NMR and X-ray crystallography.

Concluding Remarks Single-molecule spectroscopy is particularly powerful for identifying and characterizing spatially and temporally inhomogeneous protein conformational dynamics. For the last two decades, we have witnessed a rapid development of this field, and its future will definitely involve further advanced single-molecule approaches: single-molecule combined experimental multiple-parameter probing with high chemical selectivity and wide temporal and spatial VOL. 38, NO. 7, 2005 / ACCOUNTS OF CHEMICAL RESEARCH 563

Single-Molecule Protein Conformational Dynamics Lu

resolutions; computational molecular dynamics simulation analysis, based on static protein structures and correlated experimental single-molecule trajectory analyses; and theoretical modeling with experimentally accessible physical parameters. Already, single-molecule spectroscopy has substantially changed the landscape of the biological and biophysical sciences. We now have an open window for observing complex and inhomogeneous protein conformational fluctuations and motions in real time under physiological conditions. Ultimately, the singlemolecule conformational dynamics of protein-protein and protein-DNA interactions, enzymatic reactions, and ion-channel protein activities will be directly studied in living cells. The author gratefully acknowledges Dehong Hu, Yu Chen, Sofie Tan, Erich Vorpagel, Greg Harms, and Galya Orr for their crucial contributions to the work discussed here; Steve Colson, Greg Schenter, and Brian Thrall for helpful discussions. The author also thanks Prof. Brian Matthews of the University of Oregon, Prof. Klaus Hahn of the University of North Carolina School of Medicine, and Prof. Maurice Montal of University of California at San Diego for fruitful collaborations. The work was supported by the Office of Basic Energy Sciences within the Office of Science of the U.S. Department of Energy (DOE) and by the Pacific Northwest National Laboratory’s Laboratory Directed Research and Development Program. The Pacific Northwest National Laboratory is operated for the U.S. Department of Energy by Battelle Memorial Institute.

References (1) Xie, X.S.; Trautman, J. K. Optical studies of sinlge molecules at room temperature. Annu. Rev. Phys. Chem. 1998, 49, 441-480. (2) Moerner, W. E.; Orrit, M. Illuminating molecules in condensed matter. Science 1999, 283, 1670-1676. (3) Ishii, Y.; Yanagida, T. Single Molecule Detection in Life Sciences. Single Molecules 2000, 1, 5-16. (4) Rigler, R.; Orrit, M.; Basche, T. Single Molecule Spectroscopy: Nobel Conference Lectures; Springer-Verlag: Berlin, 2003. (5) Barkai, E.; Jung, Y. J.; Silbey, R. Theory of Single-Molecule Spectroscopy: Beyound the Ensemble Average. Annu. Rev. Phys. Chem. 2004, 55, 457-507. (6) Tan, X.; Nalbant, P.; Toutchkine, A.; Hu, D.; Vorpagel, E. R.; Hahn, K. M.; Lu, H. P. Single-molecule study of protein-protein interaction dynamics in a cell signaling system, J. Phys. Chem. B. 2004, 108, 737-744. (7) Lu, H. P.; Iakoucheva, L. M.; Ackerman, E. J. Single-molecule conformational dynamics of fluctuating noncovalent protein-DNA interactions in DNA damage recogntion, J. Am. Chem. Soc. 2001, 123, 9184-1985. (8) Chen, Y.; Hu, D.; Vorpegal, E. R.; Lu, H. P. Probing of singlemolecule T4 lysozyme conformational dynamics by intramolecular fluorescence energy transfer, J. Phys. Chem. B 2003, 107, 7947-7956. (9) Hu, D.; Lu, H. P. Single-molecule nanosecond anisotropy dynamics of tethered proteins. J. Phys. Chem. B 2003, 107, 618-626. (10) Hu, D.; Lu, H. P. Single Molecule Implanting of T4 Lysozyme on Bacterial Cell Surface: Towards Study Single Molecule Enzymatic Reaction in Living Cells. Biophys. J. 2004, 87, 656-661. (11) Harms, G. S.; Orr, G.; Montal, M.; Thrall, B. D.; Colson, S. D.; Lu, H. P. Probing Conformational Changes of Gramicidin Ion Channels by Single-Molecule Patch-Clamp Fluorescence Microscopy, Biophys. J. 2003, 85, 1826-1838. (12) Harms, G. S.; Orr, G.; Lu, H. P. Probing ion channel conformational dynamics using simultaneous single-molecule ultrafast spectroscopy and patch-clamp electric recording, Appl. Phys. Lett. 2004, 84, 1792-1794. (13) Lu, H. P.; Xun, L.; Xie, X. S. Single-molecule enzymatic dynamics. Science 1998, 282, 1877-1882. (14) Xie, X. S. Single-molecule approach to dispersed kinetics and dynamic disorder: Probing conformational fluctuation and enzymatic dynamics. J. Chem Phys. 2002, 117, 11024-11032. 564 ACCOUNTS OF CHEMICAL RESEARCH / VOL. 38, NO. 7, 2005

(15) Shoemaker, B. A.; Portman, J. J.; Wolynes, P. G. Speeding molecular recognition by using the folding funnel: The fly-casting mechanism. Proc. Natl. Acad. Sci. U.S.A. 2000, 97, 8868-8873. (16) For example, Wright, P. E.; Dyson, H. J. Intrinsically unstructured proteins: re-assessing the protein structure- function paradigm. J. Mol. Biol. 1999, 293, 321-331. (17) For example, Iakoucheva, L. M.; Dunker, A. K. Order, disorder, and flexibility: prediction from protein sequence. Structure 2003, 11, 1316-1317. (18) Tan, X.; Hu, D.; Squier, T. C.; Lu, H. P. Probing nanosecond protein motions of calmodulin by single-molecule fluorescence anisotropy. Appl. Phys. Lett. 2004, 85, 2420-2422. (19) Fries, J. R.; Brand, L.; Eggeling, C.; Kollner, M.; Seidel, C. A. M. Quantitative identification of different single molecules by selective time-resolved confocal fluorescence spectroscopy. J. Phys. Chem. A 1998, 102, 6601-6613. (20) Yang, H.; Xie, X. S. Probing single-molecule dynamics photon by photon. J. Chem. Phys. B 2002, 117, 10965-10979. (21) Tang, J.; Mei, E.; Green, C.; Kaplan, J.; DeGrado, W. F.; Smith, A. B., III; Hochstrasser, R. M., Probing structureal dynamics of individual calmodulin: Peptide complexes in hydrogels by singlemolecule confocal microscopy. J. Phys. Chem. B 2004, 108, 15910-15918. (22) Zwanzig, R. Rate-processes with dynamic disorder. Acc. Chem. Res. 1990, 23, 148-152. (23) Frauenfelder, H.; Sligar, S. G.; Wolynes, P. G. The energy landscapes and motions of proteins. Science 1991, 254, 15981603. (24) Matthews, B. W. Studies on protein stability with T4 lysozyme. Adv. Protein Chem. 1995, 46, 249-278. (25) Fang, C. J.; Hubbell, W. L. Conformation of T4 lysozyme in solution. Hinge-bending motion and the substrate-induced conformational transition studied by site-directed spin labeling. Biochemistry 1997, 36, 307-316. (26) Hubbell, W. L. Conformation of T4 lysozyme in solution. Hingebending motion and the substrate-induced conformational transition studied by site-directed spin labeling. Biochemistry 1997, 36, 307-316. (27) Weiss, S. Fluorescence spectroscopy of single biomolecules. Science 1999, 283, 1676-1683. (28) Lu, H. P.; Xie, X. S. Single-molecule spectral fluctuations at room temperature. Nature 1997, 385, 143-146. (29) Birks, J. B. Organic Molecular Photophysics; Wiley-VCH: New York, 1973. (30) Oppenheim, I.; Shuler, K. E.; Weiss, G. H. Stochastic Processes in Chemical Physics: The Master Equation; The MIT Press: Cambridge, MA, 1977. (31) Sakmann, B.; Neher, E. Single Channel Recordings, 2nd ed.; Kluwer: New York, 1995. Hille, B. Ion Channels of Excitable Membranes, 3rd ed.; Sinauer: Sunderland, MA, 2001. (32) Ide, T.; Yanagida, T. An artifical bilayer formed on an agarosecoated glass for simultaneoous electrical and optical measurement of single ion channels. Bioch. Biophys. Res. Commun. 1999, 265, 595-599. (33) Orr, G.; Montal, M.; Thrall, B. D.; Colson, S.; Lu, H. P. Singlechannel patch-clamp recording coupled with linear and nonlinear confocal scanning fluorescence spectroscopy: Towards the simultaneous probing of single-ion channel conformational changes and channel kinetics. Biophys. J. 2001, 80, 151a. (34) Ide, T.; Takeuchi, Y.; Yanagida, T. Development of an experimental apparatus for simultaneous observation of optical and electrical signals from single ion channels. Single Molecules 2002, 3, 3342. (35) Borisenko, V.; Lougheed, T.; Hesse, J.; Fureder-Kitzmuller, E.; Fertig, N.; Behrends, J. C.; Woolley, G. A.; Schutz, G. J. Simultaneous Optical and Electrical Recording of Single Gramicidin Channels. Biophys. J. 2003, 84, 612-622. (36) Sonnleitner, A.; Mannuzzu, L. M.; Terakawa, S.; Isacoff, E. Y. Structural rearrangements in single ion channels detected optically in living cells. Proc. Natl. Acad. Sci. U.S.A. 2002, 99, 1275912764. (37) Harms, G. S.; Cognet, L.; Lommerse, P. H. M.; Blab, G. A.; Kahr, H.; Gamsjager, R.; Spaink, H. P.; Soldator, N. M.; Romanin, C.; Schmidt, Th. Single-Molecule Imaging of L-Type Ca2+ Channels in Live Cells. Biophys. J. 2001, 81, 2639-2646. (38) Zheng, J.; Zagotta, W. N. Gating Rearrangements in Cyclic Neurotechnique Nucleotide-Gated Channels Revealed by PatchClamp Fluorometry. Neuron 2000, 28, 369-374. (39) Demuro, A.; Parker, I. Imaging the activity and localization of single voltage-gated Ca2+ channels by total internal reflection fluorescence microscopy. Biophys. J. 2004, 86, 3250-3259.

Single-Molecule Protein Conformational Dynamics Lu (40) Davis, R. W.; Patrick, E. L.; Meyer, L. A.; Ortiz, T. P.; Marshall, J. A.; Keller, D. J.; Brozik, S. M.; Brozik, J. A. Thermodynamic properties of single ion channel formation: Gramicidin, J. Phys. Chem. B 2004, 108, 15364-15369. (41) Orr, G.; Harms, G. S.; Thrall, B. D.; Montal, M.; Colson, S.; Lu, H. P. Probing Single-molecule ligand-channel interaction dynamics in a living cell. Biophys. J. 2002, 82, 255a. (42) Veatch, W. R.; Mathies, R.; Eisenberg, M.; Stryer, L. Simultaneous fluorescence and conductance studies of planar bilayer membranes containing a highly active fluorescent analog of gramicidin a. J. Mol. Biol. 1975, 99, 75-92. (43) Sawyer, D. B.; Koeppe, R. E.; Andersen, O. S. Induction of conductance heterogeneity in gramicidin channels. Biochemistry 1989, 28, 6571-6583. (44) Andersen, O. S. et al. Gramicidin channel controversy-the structure in a lopid environment. Nat. Struct. Biol. 1999, 6, 609.

(45) Milochevsky, G. V.; Gordon, P. C. Gating gramicidin channels in lipid bilayers: Reaction coordinates and the mechanism of dissociation. Biophys. J. 2004, 86, 92-104. (46) Johnson, D. I. Cdc42: An essential Rho-type GTPase controlling eukaryotic cell polarity. Microbiol. Mol. Biol. Rev. 1999, 63, 54105. (47) Abdul-Manan, N. et al. Structure of Cdc42 in complex with the GTPase-binding domain of the ‘Wiskott-Aldrich syndrome’ protein. Nature 1999, 399, 379-383. (48) Nalbant, P.; Hodgson, L.; Kraynov, V.; Toutchkine, A.; Hahn, K. M. Activation of endogenous Cdc42 visualized in living cells. Science 2004, 305, 1615-1619. (49) Biju, V.; Micic, M.; Hu, D.; Lu, H. P. Intermittent single-molecule interfacial electron-transfer dynamics. J. Am. Chem. Soc. 2004, 126, 9374-9381.

AR0401451

VOL. 38, NO. 7, 2005 / ACCOUNTS OF CHEMICAL RESEARCH 565