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.

Mechanical resistance in unstructured proteins.

Single-molecule pulling experiments on unstructured proteins linked to neurodegenerative diseases have measured rupture forces comparable to those for stable folded proteins. To investigate the structural mechanisms of this unexpected force resistance, we perform pulling simulations of the amyloid β-peptide (Aβ) and α-synuclein (αS), starting from simulated conformational ensembles for the free monomers. For both proteins, the simulations yield a set of rupture events that agree well with the experimental data.

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.

Unraveling hidden regulatory sites in structurally homologous metalloproteases.

Monitoring enzymatic activity in vivo of individual homologous enzymes such as the matrix metalloproteinases (MMPs) by antagonist molecules is highly desired for defining physiological and pathophysiological pathways. However, the rational design of antagonists targeting enzyme catalytic moieties specific to one of the homologous enzymes often appears to be an extremely difficult task. This is mainly due to the high structural homology at the enzyme active sites shared by members of the protein family.

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 as a molecular basis for allosteric regulation.

Allosteric regulation involves conformational transitions or fluctuations between a few closely related states, caused by the binding of effector molecules. We introduce a quantity called binding leverage that measures the ability of a binding site to couple to the intrinsic motions of a protein. We use Monte Carlo simulations to generate potential binding sites and either normal modes or pairs of crystal structures to describe relevant motions.

Monte Carlo study of the formation and conformational properties of dimers of Aβ42 variants.

Small soluble oligomers, and dimers in particular, of the amyloid β-peptide (Aβ) are believed to play an important pathological role in Alzheimer's disease. Here, we investigate the spontaneous dimerization of Aβ42, with 42 residues, by implicit solvent all-atom Monte Carlo simulations, for the wild-type peptide and the mutants F20E, E22G and E22G/I31E. The observed dimers of these variants share many overall conformational characteristics but differ in several aspects at a detailed level.

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.

Comparing the folding free-energy landscapes of Abeta42 variants with different aggregation properties.

The properties of the amyloid-beta peptide that lead to aggregation associated with Alzheimer's disease are not fully understood. This study aims at identifying conformational differences among four variants of full-length Abeta42 that are known to display very different aggregation properties. By extensive all-atom Monte Carlo simulations, we find that a variety of beta-sheet structures with distinct turns are readily accessible for full-length Abeta42.

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.

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.

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 (FnIII(10)) 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 three major intermediates lack two of seven native beta-strands, and share a quite similar extension.

Spontaneous beta-barrel formation: an all-atom Monte Carlo study of Abeta16-22 oligomerization.

Using all-atom Monte Carlo simulations with implicit water, combined with a cluster size analysis, we study the aggregation of Abeta(16) (-22), a peptide capable of forming amyloid fibrils. We consider a system of six initially randomly oriented Abeta(16) (-22) peptides, and investigate the thermodynamics and structural properties of aggregates formed by this system. The system is unaggregated without ordered secondary structure at high temperature, and forms beta-sheet rich aggregates at low temperature.

Differences in solution behavior among four semiconductor-binding peptides.

Recent experiments have identified peptides that adhere to GaAs and Si surfaces. Here, we use all-atom Monte Carlo simulations with implicit solvent to investigate the behavior in aqueous solution of four such peptides, all with 12 residues. At room temperature, we find that all four peptides are largely unstructured, which is consistent with experimental data.

Thermal versus mechanical unfolding of ubiquitin.

The authors studied the temperature-induced unfolding of ubiquitin by all-atom Monte Carlo simulations. The unfolding behavior is compared with that seen in previous simulations of the mechanical unfolding of this protein, based on the same model. In mechanical unfolding, secondary-structure elements were found to break in a quite well-defined order.

Dissecting the mechanical unfolding of ubiquitin.

The unfolding behavior of ubiquitin under the influence of a stretching force recently was investigated experimentally by single-molecule constant-force methods. Many observed unfolding traces had a simple two-state character, whereas others showed clear evidence of intermediate states. Here, we use Monte Carlo simulations to investigate the force-induced unfolding of ubiquitin at the atomic level.