-
Potentials of Mean Force for Protein Structure Prediction Vindicated, Formalized and Generalized
Authors:
Thomas Hamelryck,
Mikael Borg,
Martin Paluszewski,
Jonas Paulsen,
Jes Frellsen,
Christian Andreetta,
Wouter Boomsma,
Sandro Bottaro,
Jesper Ferkinghoff-Borg
Abstract:
Understanding protein structure is of crucial importance in science, medicine and biotechnology. For about two decades, knowledge based potentials based on pairwise distances -- so-called "potentials of mean force" (PMFs) -- have been center stage in the prediction and design of protein structure and the simulation of protein folding. However, the validity, scope and limitations of these potential…
▽ More
Understanding protein structure is of crucial importance in science, medicine and biotechnology. For about two decades, knowledge based potentials based on pairwise distances -- so-called "potentials of mean force" (PMFs) -- have been center stage in the prediction and design of protein structure and the simulation of protein folding. However, the validity, scope and limitations of these potentials are still vigorously debated and disputed, and the optimal choice of the reference state -- a necessary component of these potentials -- is an unsolved problem. PMFs are loosely justified by analogy to the reversible work theorem in statistical physics, or by a statistical argument based on a likelihood function. Both justifications are insightful but leave many questions unanswered. Here, we show for the first time that PMFs can be seen as approximations to quantities that do have a rigorous probabilistic justification: they naturally arise when probability distributions over different features of proteins need to be combined. We call these quantities reference ratio distributions deriving from the application of the reference ratio method. This new view is not only of theoretical relevance, but leads to many insights that are of direct practical use: the reference state is uniquely defined and does not require external physical insights; the approach can be generalized beyond pairwise distances to arbitrary features of protein structure; and it becomes clear for which purposes the use of these quantities is justified. We illustrate these insights with two applications, involving the radius of gyration and hydrogen bonding. In the latter case, we also show how the reference ratio method can be iteratively applied to sculpt an energy funnel. Our results considerably increase the understanding and scope of energy functions derived from known biomolecular structures.
△ Less
Submitted 23 November, 2010; v1 submitted 24 August, 2010;
originally announced August 2010.
-
Stop-and-go kinetics in amyloid fibrillation
Authors:
Jesper Fonslet,
Christian Beyschau Andersen,
Sandeep Krishna,
Simone Pigolotti,
Hisashi Yagi,
Yuji Goto,
Daniel Otzen,
Mogens H. Jensen,
Jesper Ferkinghoff-Borg
Abstract:
Many human diseases are associated with protein aggregation and fibrillation. We present experiments on in vitro glucagon fibrillation using total internal reflection fluorescence microscopy, providing real-time measurements of single-fibril growth. We find that amyloid fibrils grow in an intermittent fashion, with periods of growth followed by long pauses. The observed exponential distributions…
▽ More
Many human diseases are associated with protein aggregation and fibrillation. We present experiments on in vitro glucagon fibrillation using total internal reflection fluorescence microscopy, providing real-time measurements of single-fibril growth. We find that amyloid fibrils grow in an intermittent fashion, with periods of growth followed by long pauses. The observed exponential distributions of stop and growth times support a Markovian model, in which fibrils shift between the two states with specific rates. Remarkably, the probability of being in the growing (stopping) state is very close to 1/4 (3/4) in all experiments, even if the rates vary considerably. This finding suggests the presence of 4 independent conformations of the fibril tip; we discuss this possibility in terms of the existing structural knowledge.
△ Less
Submitted 21 October, 2009;
originally announced October 2009.
-
What thermodynamic features characterize good and bad folders? Results from a simplified off-lattice protein model
Authors:
A. Amatori,
J. Ferkinghoff-Borg,
G. Tiana,
R. A. Broglia
Abstract:
The thermodynamics of the small SH3 protein domain is studied by means of a simplified model where each bead-like amino acid interacts with the others through a contact potential controlled by a 20x20 random matrix. Good folding sequences, characterized by a low native energy, display three main thermodynamical phases, namely a coil-like phase, an unfolded globule and a folded phase (plus other…
▽ More
The thermodynamics of the small SH3 protein domain is studied by means of a simplified model where each bead-like amino acid interacts with the others through a contact potential controlled by a 20x20 random matrix. Good folding sequences, characterized by a low native energy, display three main thermodynamical phases, namely a coil-like phase, an unfolded globule and a folded phase (plus other two phases, namely frozen and random coil, populated only at extremes temperatures). Interestingly, the unfolded globule has some regions already structured. Poorly designed sequences, on the other hand, display a wide transition from the random coil to a frozen state. The comparison with the analytic theory of heteropolymers is discussed.
△ Less
Submitted 12 December, 2005;
originally announced December 2005.
-
Design of amino acid sequences to fold into C_alpha-model proteins
Authors:
A. Amatori,
G. Tiana,
L. Sutto,
J. Ferkinghoff-Borg,
A. Trovato,
R. A. Broglia
Abstract:
In order to extend the results obtained with minimal lattice models to more realistic systems, we study a model where proteins are described as a chain of 20 kinds of structureless amino acids moving in a continuum space and interacting through a contact potential controlled by a 20x20 quenched random matrix. The goal of the present work is to design and characterize amino acid sequences folding…
▽ More
In order to extend the results obtained with minimal lattice models to more realistic systems, we study a model where proteins are described as a chain of 20 kinds of structureless amino acids moving in a continuum space and interacting through a contact potential controlled by a 20x20 quenched random matrix. The goal of the present work is to design and characterize amino acid sequences folding to the SH3 conformation, a 60-residues recognition domain common to many regulatory proteins. We show that a number of sequences can fold, starting from a random conformation, to within a distance root mean square deviation (dRMSD) of 2.6A from the native state. Good folders are those sequences displaying in the native conformation an energy lower than a sequence--independent threshold energy.
△ Less
Submitted 16 May, 2005;
originally announced May 2005.
-
Competition between Diffusion and Fragmentation: An Important Evolutionary Process of Nature
Authors:
Jesper Ferkinghoff-Borg,
Mogens H. Jensen,
Joachim Mathiesen,
Poul Olesen,
Kim Sneppen
Abstract:
We investigate systems of nature where the common physical processes diffusion and fragmentation compete. We derive a rate equation for the size distribution of fragments. The equation leads to a third order differential equation which we solve exactly in terms of Bessel functions. The stationary state is a universal Bessel distribution described by one parameter, which fits perfectly experiment…
▽ More
We investigate systems of nature where the common physical processes diffusion and fragmentation compete. We derive a rate equation for the size distribution of fragments. The equation leads to a third order differential equation which we solve exactly in terms of Bessel functions. The stationary state is a universal Bessel distribution described by one parameter, which fits perfectly experimental data from two very different system of nature, namely, the distribution of ice crystal sizes from the Greenland ice sheet and the length distribution of alpha-helices in proteins.
△ Less
Submitted 29 December, 2003; v1 submitted 17 October, 2003;
originally announced October 2003.
-
Mean first passage time analysis reveals rate-limiting steps, parallel pathways and dead ends in a simple model of protein folding
Authors:
M. A. Micheelsen,
C. Rischel,
J. Ferkinghoff-Borg,
R. Guerois,
L. Serrano
Abstract:
We have analyzed dynamics on the complex free energy landscape of protein folding in the FOLD-X model, by calculating for each state of the system the mean first passage time to the folded state. The resulting kinetic map of the folding process shows that it proceeds in jumps between well-defined, local free energy minima. Closer analysis of the different local minima allows us to reveal seconda…
▽ More
We have analyzed dynamics on the complex free energy landscape of protein folding in the FOLD-X model, by calculating for each state of the system the mean first passage time to the folded state. The resulting kinetic map of the folding process shows that it proceeds in jumps between well-defined, local free energy minima. Closer analysis of the different local minima allows us to reveal secondary, parallel pathways as well as dead ends.
△ Less
Submitted 28 March, 2003;
originally announced March 2003.
-
Compact phases of polymers with hydrogen bonding
Authors:
Antonio Trovato,
Jesper Ferkinghoff-Borg,
Mogens H. Jensen
Abstract:
We propose an off-lattice model for a self-avoiding homopolymer chain with two different competing attractive interactions, mimicking the hydrophobic effect and the hydrogen bond formation respectively. By means of Monte Carlo simulations, we are able to trace out the complete phase diagram for different values of the relative strength of the two competing interactions. For strong enough hydroge…
▽ More
We propose an off-lattice model for a self-avoiding homopolymer chain with two different competing attractive interactions, mimicking the hydrophobic effect and the hydrogen bond formation respectively. By means of Monte Carlo simulations, we are able to trace out the complete phase diagram for different values of the relative strength of the two competing interactions. For strong enough hydrogen bonding, the ground state is a helical conformation, whereas with decreasing hydrogen bonding strength, helices get eventually destabilized at low temperature in favor of more compact conformations resembling $β$-sheets appearing in native structures of proteins. For weaker hydrogen bonding helices are not thermodynamically relevant anymore.
△ Less
Submitted 23 April, 2003; v1 submitted 5 December, 2001;
originally announced December 2001.