New force replica exchange method and protein folding pathways probed by force-clamp technique

TL;DR

This paper introduces a new force replica exchange method to study protein folding pathways under external force, revealing how terminal fixation influences refolding routes and providing insights into ubiquitin and titin domain behaviors.

Contribution

A novel extended replica exchange technique for analyzing thermodynamics of proteins under force, with findings on how terminal fixation affects folding pathways and transition states.

Findings

Refolding pathways depend on which terminus is fixed.

Anchoring the C-terminal does not alter folding pathways.

Force-extension maximum depends nonlinearly on temperature.

Abstract

We have developed a new extended replica exchange method to study thermodynamics of a system in the presence of external force. Our idea is based on the exchange between different force replicas to accelerate the equilibrium process. We have shown that the refolding pathways of single ubiquitin depend on which terminus is fixed. If the N-end is fixed then the folding pathways are different compared to the case when both termini are free, but fixing the C-terminal does not change them. Surprisingly, we have found that the anchoring terminal does not affect the pathways of individual secondary structures of three-domain ubiquitin, indicating the important role of the multi-domain construction. Therefore, force-clamp experiments, in which one end of a protein is kept fixed, can probe the refolding pathways of a single free-end ubiquitin if one uses either the poly-ubiquitin or a single domain with the C-terminus anchored. However, it is shown that anchoring one end does not affect refolding pathways of the titin domain I27, and the force-clamp spectroscopy is always capable to predict folding sequencing of this protein. We have obtained the reasonable estimate for unfolding barrier of ubiqutin. The linkage between residue Lys48 and the C-terminal of ubiquitin is found to have the dramatic effect on the location of the transition state along the end-to-end distance reaction coordinate, but the multi-domain construction leaves the transition state almost unchanged. We have found that the maximum force in the force-extension profile from constant velocity force pulling simulations depends on temperature nonlinearly. However, for some narrow temperature interval this dependence becomes linear, as have been observed in recent experiments.