Geometry of the energy landscape and folding transition in a simple model of a protein

A geometric analysis of the global properties of the energy landscape of a minimalistic model of a polypeptide is presented, which is based on the relation between dynamical trajectories and geodesics of a suitable manifold, whose metric is completely determined by the potential energy. We consider different sequences, some with a definite protein-like behavior, a unique native state and a folding transition, and the others undergoing a hydrophobic collapse with no tendency to a unique native state. The global geometry of the energy landscape appears to contain relevant information on the behavior of the various sequences: in particular, the fluctuations of the curvature of the energy landscape, measured by means of numerical simulations, clearly mark the folding transition and allow to distinguish the protein-like sequences from the others.

💡 Research Summary

The paper presents a geometric framework for analyzing the global properties of the energy landscape of a minimalistic polypeptide model, with the aim of identifying signatures of protein‑like folding behavior. Starting from the observation that the dynamics of a Hamiltonian system can be recast as geodesic motion on a suitably chosen Riemannian manifold, the authors adopt the Eisenhart metric—a pseudo‑Riemannian metric defined on an extended configuration‑time space (M × ℝ²). This metric has the crucial property that its geodesics, when projected onto the physical configuration space, reproduce the Newtonian trajectories of the original system.

Within this geometric setting, the Ricci curvature K_R in the direction of motion reduces to the Laplacian of the potential energy, ΔV. Because ΔV is simply the sum of second derivatives of the potential with respect to all coordinates, it provides a scalar measure of the average curvature “felt” by the system along its trajectory. Positive values of K_R correspond to locally convex regions (energy minima), while negative values indicate saddle‑like regions that can destabilize the motion.

To test whether this curvature carries information about folding, the authors study a three‑dimensional off‑lattice model originally introduced by Thirumalai and co‑workers. The model represents a polymer chain of N beads, each bead being one of three residue types: polar (P), hydrophobic (H), or neutral (N). The total potential V consists of bonded terms (stretching, bending, dihedral) and a non‑bonded term that depends on the residue types and their distances. By varying the sequence of residues, the model can display two qualitatively different behaviors: (i) sequences that possess a unique native state and undergo a thermally driven folding transition, and (ii) sequences that merely collapse hydrophobically without forming a well‑defined native structure.

Molecular dynamics simulations are performed over a range of temperatures. For each temperature the authors compute the time‑average ⟨K_R⟩ and its standard deviation σ_K_R along the trajectory. The results reveal a striking contrast between the two classes of sequences. In the protein‑like sequences, σ_K_R exhibits a pronounced peak at the folding temperature T_f. This peak signals a large fluctuation of the curvature as the system moves from a landscape dominated by minima (low‑temperature folded basin) to one where saddles become increasingly accessible (near the transition). In the hydrophobic‑collapse sequences, σ_K_R remains essentially flat across the whole temperature range, indicating that the curvature does not experience any abrupt change.

Traditional thermodynamic observables such as the specific heat C_v show similar smooth peaks for both types of sequences, making them insufficient to discriminate a true folding transition from a generic collapse. By contrast, the curvature fluctuation provides a clear, geometry‑based order parameter that directly reflects the presence of a “folding funnel” in the energy landscape.

The authors argue that because K_R is analytically simple (just ΔV) it can be efficiently evaluated even for more realistic potentials, opening the possibility of using curvature‑based diagnostics in protein design and in the study of larger, atomistic models. They also note that the curvature approach captures global features of the landscape that are otherwise difficult to access through local minima‑saddle analyses, which are computationally prohibitive for high‑dimensional systems.

In summary, the work demonstrates that the Ricci curvature of the Eisenhart‑induced manifold encodes essential information about the topology of the energy landscape. Its fluctuations serve as a sensitive indicator of folding transitions, allowing one to distinguish protein‑like sequences that possess a well‑defined native basin from generic polymers that only undergo a hydrophobic collapse. This geometric perspective offers a promising new tool for probing the folding problem and for guiding the design of sequences with desired folding properties.