Energy Landscape of Ubiquitin is Weakly Multidimensional

Abstract

AbstractSingle molecule pulling experiments report time-dependent changes in the extension (X) of a biomolecule as a function of the applied force (f). By fitting the data to one-dimensional analytical models of the energy landscape, the hopping rates between the folded and unfolded states in two-state folders, the height and the location of the transition state (TS) can be extracted. Although this approach is remarkably insightful, there are cases for which the energy landscape is multidimensional (catch bonds being the most prominent). To assess if the unfolding energy landscape in small single domain proteins could be one dimensional, we simulated force-induced unfolding of Ubiquitin (Ub) using the coarse-grained Self-Organized Polymer-Side Chain (SOP-SC) model. Brownian dynamics simulations using the SOP-SC model reveal that the Ub energy landscape is weakly multidimensional (WMD) governed predominantly by a single barrier. The unfolding pathway is confined to a narrow reaction pathway that could be described as diffusion in a quasi 1D X-dependent free energy profile. However, a granular analysis using the Pfold analysis, which does not assume any form for the reaction coordinate, shows that X alone does not account for the height, and more importantly, the location of the TS. The f-dependent TS location moves towards the folded state as f increases, in accord with the Hammond postulate. Our study shows that, in addition to analyzing the f-dependent hopping rates, the transition state ensemble must also be determined without resorting to X as a reaction coordinate in order to describe the unfolding energy landscapes of single domain proteins, especially if they are only WMD.TOC Graphic