Unfolding times for proteins in a force clamp

The escape process from the native valley for proteins subjected to a constant stretching force is examined using a model for a Beta-barrel. For a wide range of forces, the unfolding dynamics can be treated as one-dimensional diffusion, parametrized in terms of the end-to-end distance. In particular, the escape times can be evaluated as first passage times for a Brownian particle moving on the protein free-energy landscape, using the Smoluchowski equation. At strong forces, the unfolding process can be viewed as a diffusive drift away from the native state, while at weak forces thermal activation is the relevant mechanism. An escape-time analysis within this approach reveals a crossover from an exponential to an inverse Gaussian escape-time distribution upon passing from weak to strong forces. Moreover, a single expression valid at weak and strong forces can be devised both for the average unfolding time as well as for the corresponding variance. The analysis offers a possible explanation of recent experimental findings for ddFLN4 and ubiquitin.Comment: 6 pages, 4 figures, submitted for pubblication to Physical Review Letter

A geometry-based generic predictor for catalytic and allosteric sites

An important aspect of understanding protein allostery, and of artificial effector design, is the characterization and prediction of substrate- and effector-binding sites. To find binding sites in allosteric enzymes, many of which are oligomeric with allosteric sites at domain interfaces, we devise a local centrality measure for residue interaction graphs, which behaves well for both small/monomeric and large/multimeric proteins. The measure is purely structure based and has a clear geometrical interpretation and no free parameters. It is not biased towards typically catalytic residues, a property that is crucial when looking for non-catalytic effector sites, which are potent drug targets.Norwegian Research Council/FUGE IIacceptedVersio

Coherent Conformational Degrees of Freedom as a Structural Basis for Allosteric Communication

Conformational changes in allosteric regulation can to a large extent be described as motion along one or a few coherent degrees of freedom. The states involved are inherent to the protein, in the sense that they are visited by the protein also in the absence of effector ligands. Previously, we developed the measure binding leverage to find sites where ligand binding can shift the conformational equilibrium of a protein. Binding leverage is calculated for a set of motion vectors representing independent conformational degrees of freedom. In this paper, to analyze allosteric communication between binding sites, we introduce the concept of leverage coupling, based on the assumption that only pairs of sites that couple to the same conformational degrees of freedom can be allosterically connected. We demonstrate how leverage coupling can be used to analyze allosteric communication in a range of enzymes (regulated by both ligand binding and post-translational modifications) and huge molecular machines such as chaperones. Leverage coupling can be calculated for any protein structure to analyze both biological and latent catalytic and regulatory sites

SPACER: server for predicting allosteric communication and effects of regulation

The SPACER server provides an interactive framework for exploring allosteric communication in proteins with different sizes, degrees of oligomerization and function. SPACER uses recently developed theoretical concepts based on the thermodynamic view of allostery. It proposes easily tractable and meaningful measures that allow users to analyze the effect of ligand binding on the intrinsic protein dynamics. The server shows potential allosteric sites and allows users to explore communication between the regulatory and functional sites. It is possible to explore, for instance, potential effector binding sites in a given structure as targets for allosteric drugs. As input, the server only requires a single structure. The server is freely available at http://allostery.bii.a-star.edu.sgpublishedVersio

Changing the mechanical unfolding pathway of FnIII10 by tuning the pulling strength

We investigate the mechanical unfolding of the tenth type III domain from fibronectin, FnIII10, both at constant force and at constant pulling velocity, by all-atom Monte Carlo simulations. We observe both apparent two-state unfolding and several unfolding pathways involving one of three major, mutually exclusive intermediate states. All the three major intermediates lack two of seven native beta-strands, and share a quite similar extension. The unfolding behavior is found to depend strongly on the pulling conditions. In particular, we observe large variations in the relative frequencies of occurrence for the intermediates. At low constant force or low constant velocity, all the three major intermediates occur with a significant frequency. At high constant force or high constant velocity, one of them, with the N- and C-terminal beta-strands detached, dominates over the other two. Using the extended Jarzynski equality, we also estimate the equilibrium free-energy landscape, calculated as a function of chain extension. The application of a constant pulling force leads to a free-energy profile with three major local minima. Two of these correspond to the native and fully unfolded states, respectively, whereas the third one can be associated with the major unfolding intermediates.Comment: 15 pages, 9 figure

An effective all-atom potential for proteins

We describe and test an implicit solvent all-atom potential for simulations of protein folding and aggregation. The potential is developed through studies of structural and thermodynamic properties of 17 peptides with diverse secondary structure. Results obtained using the final form of the potential are presented for all these peptides. The same model, with unchanged parameters, is furthermore applied to a heterodimeric coiled-coil system, a mixed alpha/beta protein and a three-helix-bundle protein, with very good results. The computational efficiency of the potential makes it possible to investigate the free-energy landscape of these 49--67-residue systems with high statistical accuracy, using only modest computational resources by today's standards

FreeSASA: An open source C library for solvent accessible surface area calculations

Calculating solvent accessible surface areas (SASA) is a run-of-the-mill calculation in structural biology. Although there are many programs available for this calculation, there are no free-standing, open-source tools designed for easy tool-chain integration. FreeSASA is an open source C library for SASA calculations that provides both command-line and Python interfaces in addition to its C API. The library implements both Lee and Richards’ and Shrake and Rupley’s approximations, and is highly configurable to allow the user to control molecular parameters, accuracy and output granularity. It only depends on standard C libraries and should therefore be easy to compile and install on any platform. The library is well-documented, stable and efficient. The command-line interface can easily replace closed source legacy programs, with comparable or better accuracy and speed, and with some added functionality

Protein dynamics: aggregation and mechanical unfolding

The subject of this thesis is protein dynamics. Papers I--IV and VI study either of two different processes: mechanical unfolding and aggregation. Paper V presents a computationally efficient all-atom model for proteins, variants of which are used to perform Monte Carlo simulations in the other papers. Mechanical unfolding experiments probe properties of proteins at the single molecule level. The only information obtained the experiments is the extension and resisting force of the molecule. We perform all-atom simulations to generate a detailed description of the unfolding process. Papers I and II discuss the mechanical and thermal unfolding of the protein ubiquitin. The principal finding of Paper I is that ubiquitin unfolds through a well-defined pathway and that the experimentally observed non-obligatory unfolding intermediate lies on this pathway. Paper II compares mechanical unfolding pathways with thermal unfolding pathways. In Paper IV we study the mechanical unfolding of the protein FNIII-10 and find that it has three important, mutually exclusive, unfolding pathways and that the balance between the three can be shifted by changing the pulling strength. Paper III describes oligomerization of six-chain systems of the disease-related peptide Abeta(16-22). We find that disordered oligomers of different sizes dominate at high temperatures and as temperature is lowered, larger, more structured, oligomers form. In particular a very stable beta-barrel structure forms. Paper VI is an investigation of the effect of mutations on the folding properties of the peptide Abeta42 from Alzheimer's disease. Small aggregates of this peptide are believedto be important toxic agents. We find that a disease-related mutant peptide, with an elevated aggregation propensity, has a larger conformational diversity than the wild-type peptide, whereas a mutation that is known to inhibit aggregation has the opposite effect