Monte Carlo simulations of proteins in cages: influence of confinement on the stability of intermediate states

We present a theoretical study of the folding of small proteins inside confining potentials. Proteins are described in the framework of an effective potential model which contains the Ramachandran angles as degrees of freedom and does not need any {\it a priori} information about the native state. Hydrogen bonds, dipole-dipole- and hydrophobic interactions are taken explicitly into account. An interesting feature displayed by this potential is the presence of some intermediates between the unfolded and native states. We consider different types of confining potentials in order to study the structural properties of proteins folding inside cages with repulsive or attractive walls. Using the Wang-Landau algorithm we determine the density of states (DOS) and analyze in detail the thermodynamical properties of the confined proteins for different sizes of the cages. We show that confinement dramatically reduces the phase space available to the protein and that the presence of intermediate states can be controlled by varying the properties of the confining potential. Cages with strongly attractive walls lead to the disappearance of the intermediate states and to a two-state folding into a less stable configuration. However, cages with slightly attractive walls make the native structure more stable than in the case of pure repulsive potentials, and the folding process occurs through intermediate configurations. In order to test the metastable states we analyze the free energy landscapes as a function of the configurational energy and of the end-to-end distance as an order parameter.

💡 Research Summary

This paper presents a comprehensive theoretical investigation of how spatial confinement influences the folding behavior of small proteins, with a particular focus on the emergence and stability of intermediate states. The authors adopt an effective potential model in which the only degrees of freedom are the Ramachandran dihedral angles (ϕ, ψ) of the peptide backbone. By doing so, they avoid any a‑priori knowledge of the native structure and retain the full rotational freedom of the backbone. The potential explicitly includes three physically motivated interaction terms: (i) hydrogen‑bonding described by angle‑ and distance‑dependent functions, (ii) dipole‑dipole interactions derived from the partial charges on backbone atoms, and (iii) a non‑linear hydrophobic term that depends on solvent‑exposed surface area and pairwise distances. This combination naturally generates metastable intermediate conformations that lie between the fully unfolded coil and the native compact state.

To explore the thermodynamics of the system, the authors employ the Wang‑Landau (WL) algorithm, which iteratively estimates the density of states (DOS) across the entire energy spectrum. Unlike conventional Metropolis Monte Carlo, WL does not suffer from trapping in local minima because it forces a flat histogram in energy space, thereby providing accurate estimates of the partition function at any temperature. From the DOS they compute temperature‑dependent quantities such as the free energy, enthalpy, entropy, and heat capacity, and they construct two‑dimensional free‑energy landscapes as functions of total configurational energy and the end‑to‑end distance (R).

The confinement is modeled as a spherical cage with a hard‑core repulsive wall or a Lennard‑Jones‑type wall that can be tuned from purely repulsive to strongly attractive. Four cage radii (10 Å, 15 Å, 20 Å, 30 Å) are examined, and the wall‑protein interaction strength ε is varied from 0 (purely repulsive) to 1 kcal mol⁻¹ (strongly attractive). This systematic variation allows the authors to disentangle pure steric confinement from specific wall‑protein affinity.

Key findings are as follows:

1. Entropy Reduction and Transition Temperature Shift – As the cage radius decreases, the available conformational space shrinks dramatically, leading to a lower DOS at all energies. Consequently, the folding‑unfolding transition temperature (Tm) shifts upward and the heat‑capacity peak becomes sharper, indicating a more cooperative transition.

2. Presence of Intermediate States under Pure Repulsion – In the hard‑wall cages, the free‑energy surface exhibits a shallow local minimum separated from the native basin by two modest barriers. This corresponds to an intermediate ensemble characterized by a partially collapsed core and an end‑to‑end distance of roughly 15–20 Å. The intermediate is thermodynamically metastable and contributes a secondary shoulder in the heat‑capacity curve.

3. Stabilization of Intermediates by Weak Attractive Walls – Introducing a modest attractive component (ε≈0.2 kcal mol⁻¹) lowers the free‑energy of the intermediate relative to the unfolded state, deepening its basin and reducing the barrier to the native state. The folding pathway becomes multi‑step: unfolded → intermediate → native. The native basin is also slightly stabilized, as reflected by a modest decrease in its free‑energy relative to the unfolded ensemble.

4. Suppression of Intermediates by Strong Attraction – When the wall is strongly attractive (ε≈1 kcal mol⁻¹), the protein adheres to the cage surface, forming a non‑native compact configuration that is energetically favored over the native fold. The intermediate basin disappears; the free‑energy landscape collapses into a two‑state picture (unfolded ↔ adsorbed non‑native). The native structure becomes less stable, and the heat‑capacity curve shows a single, broader peak.

5. Thermodynamic Signatures – Heat‑capacity (Cv) profiles corroborate the landscape analysis. Purely repulsive cages display a primary peak at Tm and a smaller secondary peak associated with the intermediate. Weakly attractive cages show a broadened primary peak and an enhanced secondary feature, while strongly attractive cages exhibit only one dominant peak.

6. Two‑Dimensional Free‑Energy Maps – By plotting F(E,R) = –kBT ln P(E,R), the authors demonstrate that the intermediate basin aligns with a specific range of R, confirming that the intermediate is defined more by a geometric constraint (partial extension) than by a unique energy value. This visualization also highlights the metastable nature of the adsorbed non‑native state in the strongly attractive case.

The authors discuss the practical implications of these results for nanotechnological applications such as protein encapsulation, drug delivery, and biosensor design. For instance, a mildly hydrophilic cage coating could be used to preserve native-like folding pathways while still providing confinement, whereas a strongly hydrophobic coating would promote rapid adsorption and potentially misfolded states—an effect that could be exploited to trap proteins in a desired conformation.

In summary, the study demonstrates that confinement is a powerful lever for modulating protein folding landscapes. By combining a Ramachandran‑angle‑based effective potential with Wang‑Landau sampling, the authors provide a robust computational framework that captures both the thermodynamic and structural aspects of folding under spatial restriction. The work opens avenues for future research, including experimental validation with engineered nanocages, extension to larger or multi‑domain proteins, and integration with kinetic Monte Carlo or molecular dynamics to explore folding pathways in real time.