Mapping the energy landscape of biomolecules using single molecule force correlation spectroscopy(FCS): Theory and applications

¹ University of Maryland
² George Mason University

^* To whom correspondence should be addressed. E-mail: thirum{at}glue.umd.edu.

Submitted on October 12, 2005
Revised on November 14, 2005
Accepted on 12 January 2006

	Abstract

We present a new theory that takes internal dynamics of proteins into account to describe forced-unfolding and force-quench refolding in single molecule experiments. In the current experimental setup (Atomic Force Microscopy or Laser Optical Tweezers) the distribution of unfolding times, P(t), is measured by applying a constant stretching force f_S from which the apparent f_S-dependent unfolding rate is obtained. To describe the complexity of the underlying energy landscape requires additional probes that can incorporate the dynamics of tension propagation and relaxation of the polypeptide chain upon force quench. We introduce a theory of force correlation spectroscopy (FCS) to map the parameters of the energy landscape of proteins. In the FCS the joint distribution P(T,t) of folding and unfolding times is constructed by repeated application of cycles of stretching at constant f_S separated by release periods T during which the force is quenched to f_Q < f_S. During the release period, the protein can collapse to a manifold of compact states or refold. We show that P(T,t) at various f_S and f_Q values can be used to resolve the kinetics of unfolding as well as formation of native contacts. We also present methods to extract the parameters of the energy landscape using chain extension as the reaction coordinate and P(T,t). The theory and a worm-like chain model for the unfolded states allows us to obtain the persistence length l_p and the f_Q-dependent relaxation time, that gives an estimate of collapse timescale at the single molecular level, in the coil states of the polypeptide chain. Thus, a more complete description of landscape of protein native interactions can be maped out if unfolding time data are collected at several values of f_S and f_Q. We illustrate the utility of the proposed formalism by analyzing simulations of unfolding-refolding trajectories of a coarse-grained protein (S1) with -sheet architecture for several values of f_S, T and f_Q=0. The simulations of stretch-relax trajectories are used to map many of the parameters that characterize the energy landscape of S1.

Key Words: beta-sheet protein, coil dynamics, force correlation spectroscopy (FCS), force-induced unfolding, force-quench refolding