Proteins attain their function only after folding into a highly organized three-dimensional structure. Much remains to be learned about the mechanisms of folding of large multidomain proteins, which may populate metastable intermediate states on their energy landscapes. Here we introduce a novel method, based on high-throughput single-molecule fluorescence experiments, which is specifically geared towards tracing the dynamics of folding in the presence of a plethora of intermediates. We employ this method to characterize the folding reaction of a three-domain protein, adenylate kinase. Using thousands of single-molecule trajectories and hidden Markov modelling, we identify six metastable states on adenylate kinase’s folding landscape. Remarkably, the connectivity of the intermediates depends on denaturant concentration; at low concentration, multiple intersecting folding pathways co-exist. We anticipate that the methodology introduced here will find broad applicability in the study of folding of large proteins, and will provide a more realistic scenario of their conformational dynamics.

The function of biological macromolecules involves large-scale conformational dynamics spanning multiple time scales, from microseconds to seconds. Such conformational motions, which may involve whole domains or subunits of a protein, play a key role in allosteric regulation. There is an urgent need for experimental methods to probe the fastest of these motions. Single-molecule fluorescence experiments can in principle be used for observing such dynamics, but there is a lack of analysis methods that can extract the maximum amount of information from the data, down to the microsecond time scale. To address this issue, we introduce HMM, a maximum likelihood estimation algorithm for photon-by-photon analysis of single-molecule fluorescence resonance energy transfer (FRET) experiments. HMM is based on analytical estimators for model parameters, derived using the Baum-Welch algorithm. An efficient and effective method for the calculation of these estimators is introduced. HMM is shown to accurately retrieve the reaction times from ∼1 s to ∼10 μs and even faster when applied to simulations of freely diffusing molecules. We further apply this algorithm to single-molecule FRET data collected from Holliday junction molecules and show that at low magnesium concentrations their kinetics are as fast as ∼10 s. The new algorithm is particularly suitable for experiments on freely diffusing individual molecules and is readily incorporated into existing analysis packages. It paves the way for the broad application of single-molecule fluorescence to study ultrafast functional dynamics of biomolecules.

SignificanceThe potential effect of conformational dynamics of enzymes on their chemical steps has been intensely debated recently. We use single-molecule FRET experiments on adenylate kinase (AK) to shed new light on this question. AK closes its domains to bring its two substrate close together for reaction. We show that domain closure takes only microseconds to complete, which is two orders of magnitude faster than the chemical reaction. Nevertheless, active-site mutants that reduce the rate of domain closure also reduce the reaction rate, suggesting a connection between the two phenomena. We propose that ultrafast domain closure is used by enzymes as a mechanism to optimize mutual orientation of substrates, a novel mode of coupling between conformational dynamics and catalysis.

No abstract

Experimental tools of increasing sophistication have been employed in recent years to study protein folding and misfolding. Folding is considered a complex process, and one way to address it is by studying small proteins, which seemingly possess a simple energy landscape with essentially only two stable states, either folded or unfolded. The B1-IgG binding domain of protein L (PL) is considered a model two-state folder, based on measurements using a wide range of experimental techniques. We applied single-molecule fluorescence resonance energy transfer (FRET) spectroscopy in conjunction with a hidden Markov model analysis to fully characterize the energy landscape of PL and to extract the kinetic properties of individual molecules of the protein. Surprisingly, our studies revealed the existence of a third state, hidden under the two-state behavior of PL due to its small population, ∼7%. We propose that this minority intermediate involves partial unfolding of the two C-terminal β strands of PL. Our work demonstrates that single-molecule FRET spectroscopy can be a powerful tool for a comprehensive description of the folding dynamics of proteins, capable of detecting and characterizing relatively rare metastable states that are difficult to observe in ensemble studies.