Irina V. Gopich

Characterizing the unfolded states of proteins using single-molecule FRET spectroscopy and molecular simulations

To obtain quantitative information on the size and dynamics of unfolded proteins we combined single-molecule lifetime and intensity FRET measurements with molecular simulations. We compared the unfolded states of the 64-residue, α/β protein L and the 66-residue, all-β cold-shock protein CspTm. The average radius of gyration (Rg) calculated from FRET data on freely diffusing molecules was identical for the two unfolded proteins at guanidinium chloride concentrations >3 M, and the FRET-derived Rg of protein L agreed well with the Rg previously measured by equilibrium small-angle x-ray scattering. As the denaturant concentration was lowered, the mean FRET efficiency of the unfolded subpopulation increased, signaling collapse of the polypeptide chain, with protein L being slightly more compact than CspTm. A decrease in Rg with decreasing denaturant was also observed in all-atom molecular dynamics calculations in explicit water/urea solvent, and Langevin simulations of a simplified representation of the polypeptide suggest that collapse can result from either increased interresidue attraction or decreased excluded volume. In contrast to both the FRET and simulation results, previous time-resolved small-angle x-ray scattering experiments showed no collapse for protein L. Analysis of the donor fluorescence decay of the unfolded subpopulation of both proteins gives information about the end-to-end chain distribution and suggests that chain dynamics is slow compared with the donor life-time of ≈2 ns, whereas the bin-size independence of the small excess width above the shot noise for the FRET efficiency distributions may result from incomplete conformational averaging on even the 1-ms time scale.

Ultrafast dynamics of protein collapse from single-molecule photon statistics

We use the statistics of photon emission from single molecules to probe the ultrafast dynamics of an unfolded protein via Förster resonance energy transfer. Global reconfiguration of the chain occurs on a time scale of ≈50 ns and slows down concomitant with chain collapse under folding conditions. These diffusive dynamics provide a missing link between the phenomenological chemical kinetics commonly used in protein folding and a physical description in terms of quantitative free energy surfaces. The experiments demonstrate the potential of single-molecule methods in accessing the biologically important nanosecond time scales even in heterogeneous populations.

Effect of flexibility and cis residues in single-molecule FRET studies of polyproline

Polyproline has recently been used as a spacer between donor and acceptor chromophores to help establish the accuracy of distances determined from single-molecule Förster resonance energy transfer (FRET) measurements. This work showed that the FRET efficiency in water is higher than expected for a rigid spacer and was attributed to the flexibility of the polypeptide. Here, we investigate this issue further, using a combination of single-molecule fluorescence intensity and lifetime measurements, NMR, theory, and molecular dynamics simulations of polyproline-20 that include the dyes and their linkers to the polypeptide. NMR shows that in water ≈30% of the molecules contain internal cis prolines, whereas none are detectable in trifluoroethanol. Simulations suggest that the all-trans form of polyproline is relatively stiff, with persistence lengths of 9–13 nm using different established force fields, and that the kinks arising from internal cis prolines are primarily responsible for the higher mean FRET efficiency in water. We show that the observed efficiency histograms and distributions of donor fluorescence lifetimes are explained by the presence of multiple species with efficiencies consistent with the simulations and populations determined by NMR. In calculating FRET efficiencies from the simulation, we find that the fluctuations of the chromophores, attached to long flexible linkers, also play an important role. A similar simulation approach suggests that the flexibility of the chromophore linkers is largely responsible for the previously unexplained high value of R0 required to fit the data in the classic study of Stryer and Haugland.

Theory of photon statistics in single-molecule Förster resonance energy transfer

We present the theory for the distribution of the number of donor and acceptor photons detected in a time bin and the corresponding energy-transfer efficiency distribution obtained from single-molecule Forster resonance energy-transfer measurements. Photon counts from both immobilized and freely diffusing molecules are considered. Our starting point is the joint distribution for the donor and acceptor photons for a system described by an arbitrary kinetic scheme. This is simplified by exploiting the time scale separation between fast fluorescent transitions and slow processes which include conformational dynamics, intersystem conversion to a dark state, and translational diffusion in and out of the laser spot. The fast fluorescent transitions result in a Poisson distribution of the number of photons which is then averaged over slow fluctuations of the local transfer efficiency and the total number of photons. The contribution of various processes to the distribution and the variance of the energy-transfer efficiency are analyzed.

Theory of the statistics of kinetic transitions with application to single-molecule enzyme catalysis

Single-molecule spectroscopy can monitor transitions between two microscopic states when these transitions are associated with the emission of photons. A general formalism is developed for obtaining the statistics of such transitions from a microscopic model when the dynamics is described by master or rate equations or their continuum analog, multidimensional reaction-diffusion equations. The focus is on the distribution of the number of transitions during a fixed observation time, the distribution of times between transitions, and the corresponding correlation functions. It is shown how these quantities are related to each other and how they can be explicitly calculated in a straightforward way for both immobile and diffusing molecules. Our formalism reduces to renewal theory when the monitored transitions either go to or originate from a single state. The influence of dynamics slow compared with the time between monitored transitions is treated in a simple way, and the probability distributions are expressed in terms of Mandel-type formulas. The formalism is illustrated by a detailed analysis of the statistics of catalytic turnovers of enzymes. When the rates of conformational changes are slower than the catalytic rates which are in turn slower than the binding relaxation rate, (1) the mean number of turnovers is shown to have the classical Michaelis-Menten form, (2) the correlation function of the number of turnovers is a direct measure of the time scale of catalytic rate fluctuations, and (3) the distribution of the time between consecutive turnovers is determined by the steady-state distribution.

Decoding the pattern of photon colors in single-molecule FRET.

Conformational dynamics of a single molecule can be studied using Forster resonance energy transfer (FRET) by recording a sequence of photons emitted by a donor and an acceptor dye attached to the molecule. We describe a simple and robust method to estimate the rates of transitions between different conformational states and the FRET efficiencies associated with these states. For a photon trajectory with measured interphoton times, the pattern of colors is decoded by maximizing the appropriate likelihood function. This approach can be used to analyze bursts of photons from diffusing molecules as well as photon trajectories generated by immobilized molecules. The procedure is illustrated using simulated photon trajectories corresponding to two-state and three-state molecules. The method works even when the photon colors appear to be scrambled because of high background noise, the photophysical properties of the conformers are similar, or the conformational and photon count rates are comparable. The consistency of the model with the data can be checked by recoloring the photon trajectories and comparing the predicted and observed FRET efficiency histograms.

Kinetics of reversible diffusion influenced reactions: The self-consistent relaxation time approximation

The simplest general theory of the kinetics of reversible diffusion-influenced reactions that is exact both at short and long times for A+B⇌C and A+B⇌C+D is presented. The formalism is based on an approximate set of reaction-diffusion equations for the pair distribution functions which incorporate the influence of the chemical reaction by using effective rate constants that are determined self-consistently. For small deviations from equilibrium and contact reactivity, the relaxation function is given explicitly in the Laplace domain in terms of the Smoluchowski rate coefficient that describes the corresponding diffusion controlled irreversible reaction. Consequently, the kinetics can be easily obtained for arbitrary diffusion coefficients and equilibrium concentrations.

Extracting rate coefficients from single-molecule photon trajectories and FRET efficiency histograms for a fast-folding protein.

Recently developed statistical methods by Gopich and Szabo were used to extract folding and unfolding rate coefficients from single-molecule Förster resonance energy transfer (FRET) data for proteins with kinetics too fast to measure waiting time distributions. Two types of experiments and two different analyses were performed. In one experiment bursts of photons were collected from donor and acceptor fluorophores attached to a 73-residue protein, α(3)D, freely diffusing through the illuminated volume of a confocal microscope system. In the second, the protein was immobilized by linkage to a surface, and photons were collected until one of the fluorophores bleached. Folding and unfolding rate coefficients and mean FRET efficiencies for the folded and unfolded subpopulations were obtained from a photon by photon analysis of the trajectories using a maximum likelihood method. The ability of the method to describe the data in terms of a two-state model was checked by recoloring the photon trajectories with the extracted parameters and comparing the calculated FRET efficiency histograms with the measured histograms. The sum of the rate coefficients for the two-state model agreed to within 30% with the relaxation rate obtained from the decay of the donor-acceptor cross-correlation function, confirming the high accuracy of the method. Interestingly, apparently reliable rate coefficients could be extracted using the maximum likelihood method, even at low (<10%) population of the minor component where the cross-correlation function was too noisy to obtain any useful information. The rate coefficients and mean FRET efficiencies were also obtained in an approximate procedure by simply fitting the FRET efficiency histograms, calculated by binning the donor and acceptor photons, with a sum of three-Gaussian functions. The kinetics are exposed in these histograms by the growth of a FRET efficiency peak at values intermediate between the folded and unfolded peaks as the bin size increases, a phenomenon with similarities to NMR exchange broadening. When comparable populations of folded and unfolded molecules are present, this method yields rate coefficients in very good agreement with those obtained with the maximum likelihood method. As a first step toward characterizing transition paths, the Viterbi algorithm was used to locate the most probable transition points in the photon trajectories.

Theory of the energy transfer efficiency and fluorescence lifetime distribution in single-molecule FRET

In single-molecule FRET experiments with pulsed lasers, not only the colors of the photons but also the fluorescence lifetimes can be monitored. Although these quantities appear to be random, they are modulated by conformational dynamics. In order to extract information about such dynamics, we develop the theory of the joint distribution of FRET efficiencies and fluorescence lifetimes determined from bins (or bursts) of photons. Our starting point is a rigorous formal expression for the distribution of the numbers of donor and acceptor photons and donor lifetimes in a bin that treats the influence of conformational dynamics on all timescales. This formula leads to an analytic result for a two-state system interconverting on a timescale slower than the interphoton time and to an efficient simulation algorithm for multistate dynamics. The shape of the joint distribution contains more information about conformational dynamics than the FRET efficiency histogram alone. In favorable cases, the connectivity of the underlying conformational states can be determined directly by simple inspection of the projection of the joint distribution on the efficiency-lifetime plane.

Statistics of transitions in single molecule kinetics

The probability distribution of observing N state-to-state transitions in a finite time t is calculated for an arbitrary kinetic scheme. In the Laplace domain, this probability distribution can be expressed analytically in terms of the (sI−K)−1 where K is the rate matrix. For long times simple expressions are found for the mean and variance of the number of transitions.